datatypes/int/sc_nbutils.cpp

644   sc_digit carry = 0;   // Also used as sum to save space.
648     carry += (*u++) + (*v++);
649     (*w++) = carry & DIGIT_MASK;
650     carry >>= BITS_PER_DIGIT;
653   // Propagate the carry.
654   while (carry && (u < uend)) {
655     carry = (*u++) + 1;
656     (*w++) = carry & DIGIT_MASK;
657     carry >>= BITS_PER_DIGIT;
664   // Propagate the carry if it is still 1.
665   if (carry)
688   sc_digit carry = 0;   // Also used as sum to save space.
692     carry += (*u) + (*v++);
693     (*u++) = carry & DIGIT_MASK;
694     carry >>= BITS_PER_DIGIT;
697   // Propagate the carry.
698   while (carry && (u < uend)) {
699     carry = (*u) + 1;
700     (*u++) = carry & DIGIT_MASK;
701     carry >>= BITS_PER_DIGIT;
705   if( carry != 0 ) {
736   sc_digit carry = 0;   // Also used as sum to save space.
740     carry += (*u) + (*v++);
741     (*u++) = carry & DIGIT_MASK;
742     carry >>= BITS_PER_DIGIT;
746   if( carry != 0 ) {
771   sc_digit carry = (*u++) + v;
772   (*w++) = carry & DIGIT_MASK;
773   carry >>= BITS_PER_DIGIT;
775   // Propagate the carry.
776   while (carry && (u < uend)) {
777     carry = (*u++) + 1;
778     (*w++) = carry & DIGIT_MASK;
779     carry >>= BITS_PER_DIGIT;
786   // Propagate the carry if it is still 1.
787   if (carry)
1040      carry = A * C + high_half(prod_h)
1053      NOTE: sc_max(prod_l, prod_h, carry) <= 2 * x^2 - 1, so any
1067 #define prod_h carry
1083     sc_digit carry = 0;
1101       sc_digit prod_l = (*w) + u_l * v_l + low_half(carry);
1103       prod_h = u_h * v_l + u_l * v_h + high_half(prod_l) + high_half(carry);
1107       carry = u_h * v_h + high_half(prod_h);
1111     (*w) = carry;
1132 #define prod_h carry
1136   sc_digit carry = 0;
1147     sc_digit prod_l = v * low_half(u_AB) + low_half(carry);
1149     prod_h = v * high_half(u_AB) + high_half(prod_l) + high_half(carry);
1153     carry = high_half(prod_h);
1157   (*w) = carry;
1174 #define prod_h carry
1176   sc_digit carry = 0;
1185     sc_digit prod_l = v * low_half(u[i]) + low_half(carry);
1187     prod_h = v * high_half(u[i]) + high_half(prod_l) + high_half(carry);
1191     carry = high_half(prod_h);
1198   if( carry != 0 ) {
1295       sc_digit carry = 0;
1298         carry += y[i] * qk;
1299         sc_digit diff = (xk[i] + BYTE_RADIX) - (carry & BYTE_MASK);
1301         carry = (carry >> BITS_PER_BYTE) + (1 - (diff >> BITS_PER_BYTE));
1304       // If carry, qk may be one too large.
1305       if (carry) {
1308         carry = (xk[ylen] + BYTE_RADIX) - carry;
1309         xk[ylen] = (uchar)(carry & BYTE_MASK);
1310         carry = 1 - (carry >> BITS_PER_BYTE);
1312         if (carry) {
1319           carry = 0;
1322             carry += xk[i] + y[i];
1323             xk[i] = (uchar)(carry & BYTE_MASK);
1324             carry >>= BITS_PER_BYTE;
1327           if (carry)
1330         }  // second if carry
1331       }  // first if carry
1496       sc_digit carry = 0;
1499         carry += y[i] * qk;
1500         sc_digit diff = (xk[i] + BYTE_RADIX) - (carry & BYTE_MASK);
1502         carry = (carry >> BITS_PER_BYTE) + (1 - (diff >> BITS_PER_BYTE));
1505       if (carry) {
1508         carry = (xk[ylen] + BYTE_RADIX) - carry;
1509         xk[ylen] = (uchar)(carry & BYTE_MASK);
1510         carry = 1 - (carry >> BITS_PER_BYTE);
1512         if (carry) {
1518           carry = 0;
1521             carry += xk[i] + y[i];
1522             xk[i] = (uchar)(carry & BYTE_MASK);
1523             carry >>= BITS_PER_BYTE;
1526           if (carry)
1529         }  // second if carry
1530       } // first if carry
1691     sc_digit carry = 0;
1695       (*viter++) = (((vval & mask) << BITS_PER_BYTE) | carry);
1696       carry = vval >> nsr;
1700       (*viter) = carry;
1757   sc_digit carry = 0;
1761     (*uiter++) = (((uval & mask) << nsl) | carry);
1762     carry = uval >> nsr;
1766     (*uiter) = carry;
1826   sc_digit carry = (fill & mask) << nsl;
1830     (*uiter) = (uval >> nsr) | carry;
1831     carry = (uval & mask) << nsl;