dt/int/sc_nbutils.cc

619 // - ulen >= vlen
620 // - wlen >= sc_max(ulen, vlen) + 1
622 vec_add(int ulen, const sc_digit *u, int vlen, const sc_digit *v, sc_digit *w)
625     sc_assert((ulen > 0) && (u != NULL));
628     sc_assert(ulen >= vlen);
631     const sc_digit *uend = (u + ulen);
661 // - ulen >= vlen
663 vec_add_on(int ulen, sc_digit *ubegin, int vlen, const sc_digit *v)
666     sc_assert((ulen > 0) && (ubegin != NULL));
668     sc_assert(ulen >= vlen);
672     const sc_digit *uend = (u + ulen);
703 // - ulen < vlen
705 vec_add_on2(int ulen, sc_digit *ubegin, int,
712     sc_assert((ulen > 0) && (ubegin != NULL));
714     sc_assert(ulen < vlen);
718     const sc_digit *uend = (u + ulen);
742 vec_add_small(int ulen, const sc_digit *u, sc_digit v, sc_digit *w)
745     sc_assert((ulen > 0) && (u != NULL));
749     const sc_digit *uend = (u + ulen);
774 vec_add_small_on(int ulen, sc_digit *u, sc_digit v)
777     sc_assert((ulen > 0) && (u != NULL));
782     while (v && (i < ulen)) {
799 // - ulen >= vlen
800 // - wlen >= sc_max(ulen, vlen)
802 vec_sub(int ulen, const sc_digit *u, int vlen, const sc_digit *v, sc_digit *w)
805     sc_assert((ulen > 0) && (u != NULL));
808     sc_assert(ulen >= vlen);
811     const sc_digit *uend = (u + ulen);
841 // - ulen >= vlen
843 vec_sub_on(int ulen, sc_digit *ubegin, int vlen, const sc_digit *v)
846     sc_assert((ulen > 0) && (ubegin != NULL));
848     sc_assert(ulen >= vlen);
852     const sc_digit *uend = (u + ulen);
878 // - ulen <= vlen or ulen > ulen
880 vec_sub_on2(int ulen, sc_digit *ubegin, int vlen, const sc_digit *v)
883     sc_assert((ulen > 0) && (ubegin != NULL));
888     const sc_digit *uend = (u + sc_min(ulen, vlen));
911 vec_sub_small(int ulen, const sc_digit *u, sc_digit v, sc_digit *w)
914     sc_assert(ulen > 0);
918     const sc_digit *uend = (u + ulen);
944 vec_sub_small_on(int ulen, sc_digit *u, sc_digit v)
947     sc_assert((ulen > 0) && (u != NULL));
950     for (int i = 0; i < ulen; ++i) {
963 vec_mul(int ulen, const sc_digit *u, int vlen, const sc_digit *vbegin,
1024     sc_assert((ulen > 0) && (u != NULL));
1030     const sc_digit *uend = (u + ulen);
1071 vec_mul_small(int ulen, const sc_digit *u, sc_digit v, sc_digit *w)
1074     sc_assert((ulen > 0) && (u != NULL));
1081     const sc_digit *uend = (u + ulen);
1101 vec_mul_small_on(int ulen, sc_digit *u, sc_digit v)
1104     sc_assert((ulen > 0) && (u != NULL));
1110     for (int i = 0; i < ulen; ++i) {
1135 vec_div_large(int ulen, const sc_digit *u, int vlen, const sc_digit *v,
1139     sc_assert((ulen > 0) && (u != NULL));
1151     int xlen = BYTES_PER_DIGIT * ulen + 1;
1168     xlen = vec_to_char(ulen, u, xlen, x);
1260     vec_from_char(xlen - ylen + 1, q, ulen, w);
1273 vec_div_small(int ulen, const sc_digit *u, sc_digit v, sc_digit *q)
1300     sc_assert((ulen > 0) && (u != NULL));
1309     u += ulen;
1310     q += ulen;
1332 vec_rem_large(int ulen, const sc_digit *u, int vlen, const sc_digit *v,
1336     sc_assert((ulen > 0) && (u != NULL));
1343     int xlen = BYTES_PER_DIGIT * ulen + 1;
1357     xlen = vec_to_char(ulen, u, xlen, x);
1444     vec_from_char(ylen, x, ulen, w);
1457 vec_rem_small(int ulen, const sc_digit *u, sc_digit v)
1461     sc_assert((ulen > 0) && (u != NULL));
1470     u += ulen;
1487 vec_rem_on_small(int ulen, sc_digit *u, sc_digit v)
1490     sc_assert((ulen > 0) && (u != NULL));
1498     u += ulen;
1517 vec_to_char(int ulen, const sc_digit *u, int vlen, uchar *v)
1520     sc_assert((ulen > 0) && (u != NULL));
1524     int nbits = ulen * BITS_PER_DIGIT;
1536             if (left_digit < ulen)
1552 vec_from_char(int ulen, const uchar *u, int vlen, sc_digit *v)
1555     sc_assert((ulen > 0) && (u != NULL));
1565     (*v) = (sc_digit) u[ulen - 1];
1567     for (int i = ulen - 2; i >= 0; --i) {
1587 vec_shift_left(int ulen, sc_digit *u, int nsl)
1590     sc_assert((ulen > 0) && (u != NULL));
1609             for (int j = ulen - 1; j >= nd; --j)
1612             vec_zero(sc_min(nd, ulen), u);
1620     sc_digit *uend = uiter + ulen;
1640 vec_shift_right(int ulen, sc_digit *u, int nsr, sc_digit fill)
1643     sc_assert((ulen > 0) && (u != NULL));
1663             for (int j = 0; j < (ulen - nd); ++j)
1667                 for (int j = ulen - sc_min( nd, ulen ); j < ulen; ++j)
1670                 vec_zero(ulen - sc_min( nd, ulen ), ulen, u);
1679     sc_digit *uiter = (ubegin + ulen);