| /dragonfly/contrib/gmp/mpf/ |
| HD | sub.c | 29 mp_size_t usize, vsize, rsize; in mpf_sub() local 36 usize = u->_mp_size; in mpf_sub() 40 if (usize == 0) in mpf_sub() 53 if ((usize ^ vsize) < 0) in mpf_sub() 66 negate = usize < 0; in mpf_sub() 74 usize = u->_mp_size; in mpf_sub() 78 usize = ABS (usize); in mpf_sub() 95 if (up[usize - 1] == vp[vsize - 1]) in mpf_sub() 100 usize--; in mpf_sub() 104 if (usize == 0) in mpf_sub() [all …]
|
| HD | ui_sub.c | 29 mp_size_t usize, vsize, rsize; in mpf_ui_sub() local 71 usize = ABS (vsize); in mpf_ui_sub() 83 usize = 1; in mpf_ui_sub() 100 usize--; in mpf_ui_sub() 102 if (up[usize] != vp[vsize]) in mpf_ui_sub() 105 if (usize == 0) in mpf_ui_sub() 110 usize++; in mpf_ui_sub() 118 while (usize != 0 && up[usize - 1] == 0) in mpf_ui_sub() 119 usize--, uexp--; in mpf_ui_sub() 120 if (usize == 0) in mpf_ui_sub() [all …]
|
| HD | add_ui.c | 28 mp_size_t usize, sumsize; in mpf_add_ui() local 32 usize = u->_mp_size; in mpf_add_ui() 33 if (usize <= 0) in mpf_add_ui() 35 if (usize == 0) in mpf_add_ui() 43 u_negated._mp_size = -usize; in mpf_add_ui() 57 sumsize = MIN (usize, prec + 1); in mpf_add_ui() 58 MPN_COPY (sum->_mp_d, up + usize - sumsize, sumsize); in mpf_add_ui() 78 if (uexp > usize) in mpf_add_ui() 84 MPN_COPY_DECR (sump + uexp - usize, up, usize); in mpf_add_ui() 86 MPN_ZERO (sump + 1, uexp - usize - 1); in mpf_add_ui() [all …]
|
| HD | add.c | 28 mp_size_t usize, vsize, rsize; in mpf_add() local 36 usize = u->_mp_size; in mpf_add() 40 if (usize == 0) in mpf_add() 54 if ((usize ^ vsize) < 0) in mpf_add() 67 negate = usize < 0; in mpf_add() 74 usize = u->_mp_size; in mpf_add() 78 usize = ABS (usize); in mpf_add() 88 if (usize > prec) in mpf_add() 90 up += usize - prec; in mpf_add() 91 usize = prec; in mpf_add() [all …]
|
| HD | cmp_si.c | 28 mp_size_t usize; in mpf_cmp_si() local 35 usize = u->_mp_size; in mpf_cmp_si() 38 if ((usize < 0) == (vval < 0)) /* don't use xor, type size may differ */ in mpf_cmp_si() 41 if (usize == 0) in mpf_cmp_si() 46 return usize != 0; in mpf_cmp_si() 52 return usize >= 0 ? 1 : -1; in mpf_cmp_si() 57 usign = usize >= 0 ? 1 : -1; in mpf_cmp_si() 58 usize = ABS (usize); in mpf_cmp_si() 76 ulimb = up[usize - 1]; in mpf_cmp_si() 78 if (usize >= 2 && uexp == 2) in mpf_cmp_si() [all …]
|
| HD | cmp.c | 27 mp_size_t usize, vsize; in mpf_cmp() local 35 usize = u->_mp_size; in mpf_cmp() 39 if ((usize ^ vsize) >= 0) in mpf_cmp() 42 if (usize == 0) in mpf_cmp() 47 return usize != 0; in mpf_cmp() 53 return usize >= 0 ? 1 : -1; in mpf_cmp() 58 usign = usize >= 0 ? 1 : -1; in mpf_cmp() 66 usize = ABS (usize); in mpf_cmp() 78 usize--; in mpf_cmp() 87 if (usize > vsize) in mpf_cmp() [all …]
|
| HD | div.c | 50 mp_size_t usize, vsize, rsize, prospective_rsize, tsize, zeros; in mpf_div() local 56 usize = SIZ(u); in mpf_div() 58 sign_quotient = usize ^ vsize; in mpf_div() 59 usize = ABS (usize); in mpf_div() 66 if (usize == 0) in mpf_div() 80 prospective_rsize = usize - vsize + 1; /* quot from using given u,v sizes */ in mpf_div() 88 usize -= chop; in mpf_div() 91 tsize = usize + zeros; /* size for possible copy of u */ in mpf_div() 98 MPN_COPY (tp+zeros, up, usize); in mpf_div() 100 usize = tsize; in mpf_div() [all …]
|
| HD | mul.c | 27 mp_size_t usize, vsize; in mpf_mul() local 33 usize = u->_mp_size; in mpf_mul() 35 sign_product = usize ^ vsize; in mpf_mul() 37 usize = ABS (usize); in mpf_mul() 42 if (usize > prec) in mpf_mul() 44 up += usize - prec; in mpf_mul() 45 usize = prec; in mpf_mul() 53 if (usize == 0 || vsize == 0) in mpf_mul() 65 rsize = usize + vsize; in mpf_mul() 67 cy_limb = (usize >= vsize in mpf_mul() [all …]
|
| HD | cmp_ui.c | 27 mp_size_t usize; in mpf_cmp_ui() local 32 usize = u->_mp_size; in mpf_cmp_ui() 35 if (usize < 0) in mpf_cmp_ui() 40 return usize != 0; in mpf_cmp_ui() 57 ulimb = up[usize - 1]; in mpf_cmp_ui() 59 if (usize >= 2 && uexp == 2) in mpf_cmp_ui() 63 ulimb = (ulimb << GMP_NUMB_BITS) | up[usize - 2]; in mpf_cmp_ui() 64 usize--; in mpf_cmp_ui() 67 usize--; in mpf_cmp_ui() 79 usize--; in mpf_cmp_ui() [all …]
|
| HD | eq.c | 29 mp_size_t usize, vsize, minsize, maxsize, n_limbs, i, size; in mpf_eq() local 37 usize = u->_mp_size; in mpf_eq() 41 if ((usize ^ vsize) >= 0) in mpf_eq() 44 if (usize == 0) in mpf_eq() 63 usize = ABS (usize); in mpf_eq() 69 up += usize; /* point just above most significant limb */ in mpf_eq() 79 usize = MIN (usize, n_limbs); in mpf_eq() 85 up++, usize--; in mpf_eq() 90 minsize = MIN (usize, vsize); in mpf_eq() 91 maxsize = usize + vsize - minsize; in mpf_eq() [all …]
|
| HD | div_ui.c | 30 mp_size_t usize; in mpf_div_ui() local 53 usize = u->_mp_size; in mpf_div_ui() 54 sign_quotient = usize; in mpf_div_ui() 55 usize = ABS (usize); in mpf_div_ui() 61 if (usize == 0) in mpf_div_ui() 76 if (usize > tsize) in mpf_div_ui() 78 up += usize - tsize; in mpf_div_ui() 79 usize = tsize; in mpf_div_ui() 84 MPN_ZERO (tp, tsize - usize); in mpf_div_ui() 85 rtp = tp + (tsize - usize); in mpf_div_ui() [all …]
|
| HD | sqrt.c | 55 mp_size_t usize; in mpf_sqrt() local 61 usize = u->_mp_size; in mpf_sqrt() 62 if (usize <= 0) in mpf_sqrt() 64 if (usize < 0) in mpf_sqrt() 87 if (usize > tsize) in mpf_sqrt() 89 up += usize - tsize; in mpf_sqrt() 90 usize = tsize; in mpf_sqrt() 95 MPN_ZERO (tp, tsize - usize); in mpf_sqrt() 96 MPN_COPY (tp + (tsize - usize), up, usize); in mpf_sqrt()
|
| /dragonfly/contrib/gmp/mpz/ |
| HD | mul.c | 36 mp_size_t usize; in mpz_mul() local 47 usize = SIZ (u); in mpz_mul() 49 sign_product = usize ^ vsize; in mpz_mul() 50 usize = ABS (usize); in mpz_mul() 53 if (usize < vsize) in mpz_mul() 56 MP_SIZE_T_SWAP (usize, vsize); in mpz_mul() 68 MPZ_REALLOC (w, usize+vsize); in mpz_mul() 71 cy_limb = mpn_mul_1 (wp, PTR(u), usize, PTR(v)[0]); in mpz_mul() 74 cy_limb = mpn_mul_2 (wp, PTR(u), usize, PTR(v)); in mpz_mul() 75 usize++; in mpz_mul() [all …]
|
| HD | hamdist.c | 29 mp_size_t usize, vsize; in mpz_hamdist() local 32 usize = SIZ(u); in mpz_hamdist() 38 if (usize >= 0) in mpz_hamdist() 45 if (usize < vsize) in mpz_hamdist() 46 MPN_SRCPTR_SWAP (up,usize, vp,vsize); in mpz_hamdist() 52 usize -= vsize; in mpz_hamdist() 53 if (usize != 0) in mpz_hamdist() 54 count += mpn_popcount (up + vsize, usize); in mpz_hamdist() 68 usize = -usize; in mpz_hamdist() 74 ASSERT (usize > 0); in mpz_hamdist() [all …]
|
| HD | gcd.c | 40 mp_size_t usize = ABS (u->_mp_size); in mpz_gcd() local 47 if (usize == 0) in mpz_gcd() 61 g->_mp_size = usize; in mpz_gcd() 64 if (g->_mp_alloc < usize) in mpz_gcd() 65 _mpz_realloc (g, usize); in mpz_gcd() 66 MPN_COPY (g->_mp_d, up, usize); in mpz_gcd() 70 if (usize == 1) in mpz_gcd() 80 g->_mp_d[0] = mpn_gcd_1 (up, usize, vp[0]); in mpz_gcd() 90 usize -= u_zero_limbs; in mpz_gcd() 93 up = TMP_ALLOC_LIMBS (usize); in mpz_gcd() [all …]
|
| HD | lcm.c | 29 mp_size_t usize, vsize, size; in mpz_lcm() local 32 usize = SIZ (u); in mpz_lcm() 34 if (usize == 0 || vsize == 0) in mpz_lcm() 39 usize = ABS (usize); in mpz_lcm() 49 MPZ_REALLOC (r, usize+1); in mpz_lcm() 53 gl = mpn_gcd_1 (up, usize, vl); in mpz_lcm() 57 c = mpn_mul_1 (rp, up, usize, vl); in mpz_lcm() 58 rp[usize] = c; in mpz_lcm() 59 usize += (c != 0); in mpz_lcm() 60 SIZ(r) = usize; in mpz_lcm() [all …]
|
| HD | lcm_ui.c | 28 mp_size_t usize; in mpz_lcm_ui() local 49 usize = SIZ(u); in mpz_lcm_ui() 50 if (usize == 0 || v == 0) in mpz_lcm_ui() 55 usize = ABS(usize); in mpz_lcm_ui() 57 MPZ_REALLOC (r, usize+1); in mpz_lcm_ui() 60 g = (unsigned long) mpn_gcd_1 (up, usize, (mp_limb_t) v); in mpz_lcm_ui() 64 c = mpn_mul_1 (rp, up, usize, (mp_limb_t) v); in mpz_lcm_ui() 65 rp[usize] = c; in mpz_lcm_ui() 66 usize += (c != 0); in mpz_lcm_ui() 67 SIZ(r) = usize; in mpz_lcm_ui()
|
| HD | gcdext.c | 29 mp_size_t asize, bsize, usize, vsize; in mpz_gcdext() local 51 usize = asize; in mpz_gcdext() 53 up = TMP_ALLOC_LIMBS (usize + 1); in mpz_gcdext() 55 MPN_COPY (up, ap, usize); in mpz_gcdext() 64 usize = bsize; in mpz_gcdext() 66 up = TMP_ALLOC_LIMBS (usize + 1); in mpz_gcdext() 68 MPN_COPY (up, bp, usize); in mpz_gcdext() 76 tmp_gp = TMP_ALLOC_LIMBS (usize + 1); in mpz_gcdext() 77 tmp_sp = TMP_ALLOC_LIMBS (usize + 1); in mpz_gcdext() 83 MPN_COPY (tmp_gp, up, usize); in mpz_gcdext() [all …]
|
| HD | aorsmul.c | 101 mp_size_t usize = wsize; in mpz_aorsmul() local 103 if (usize < tsize) in mpz_aorsmul() 106 usize = tsize; in mpz_aorsmul() 110 wsize = usize; in mpz_aorsmul() 113 c = mpn_add (wp, up,usize, tp,tsize); in mpz_aorsmul() 120 mp_size_t usize = wsize; in mpz_aorsmul() local 122 if (mpn_cmp_twosizes_lt (up,usize, tp,tsize)) in mpz_aorsmul() 125 usize = tsize; in mpz_aorsmul() 129 wsize = usize; in mpz_aorsmul() 133 ASSERT_NOCARRY (mpn_sub (wp, up,usize, tp,tsize)); in mpz_aorsmul() [all …]
|
| HD | aors.h | 61 mp_size_t usize, vsize, wsize; variable 65 usize = u->_mp_size; 67 abs_usize = ABS (usize); 74 MP_SIZE_T_SWAP (usize, vsize); 90 if ((usize ^ vsize) < 0) 101 if (usize < 0) 109 if (usize >= 0) 117 if (usize < 0) 127 if (usize < 0)
|
| HD | cfdiv_r_2exp.c | 36 mp_size_t usize, abs_usize, limb_cnt, i; in cfdiv_r_2exp() local 41 usize = SIZ(u); in cfdiv_r_2exp() 42 if (usize == 0) in cfdiv_r_2exp() 50 abs_usize = ABS (usize); in cfdiv_r_2exp() 56 if ((usize ^ dir) < 0) in cfdiv_r_2exp() 77 SIZ(w) = usize; in cfdiv_r_2exp() 120 usize = -usize; in cfdiv_r_2exp() 141 SIZ(w) = (usize >= 0 ? limb_cnt : -limb_cnt); in cfdiv_r_2exp()
|
| HD | cmp_si.c | 28 mp_size_t usize = u->_mp_size; in _mpz_cmp_si() local 52 if (usize != vsize) in _mpz_cmp_si() 53 return usize - vsize; in _mpz_cmp_si() 55 if (usize == 0) in _mpz_cmp_si() 64 return usize; in _mpz_cmp_si() 66 return -usize; in _mpz_cmp_si()
|
| HD | cfdiv_q_2exp.c | 33 mp_size_t wsize, usize, abs_usize, limb_cnt, i; in cfdiv_q_2exp() local 38 usize = SIZ (u); in cfdiv_q_2exp() 39 abs_usize = ABS (usize); in cfdiv_q_2exp() 46 SIZ(w) = (usize == 0 || (usize ^ dir) < 0 ? 0 : dir); in cfdiv_q_2exp() 57 rmask = ((usize ^ dir) >= 0 ? MP_LIMB_T_MAX : 0); in cfdiv_q_2exp() 88 SIZ(w) = (usize >= 0 ? wsize : -wsize); in cfdiv_q_2exp()
|
| /dragonfly/contrib/gmp/mpn/generic/ |
| HD | perfsqr.c | 96 #define PERFSQR_MOD_34(r, up, usize) \ argument 98 (r) = mpn_mod_34lsub1 (up, usize); \ 107 #define PERFSQR_MOD_PP(r, up, usize) \ argument 109 if (BELOW_THRESHOLD (usize, PREINV_MOD_1_TO_MOD_1_THRESHOLD)) \ 111 (r) = mpn_preinv_mod_1 (up, usize, PERFSQR_PP_NORM, \ 117 (r) = mpn_mod_1 (up, usize, PERFSQR_PP); \ 169 mpn_perfect_square_p (mp_srcptr up, mp_size_t usize) in mpn_perfect_square_p() argument 171 ASSERT (usize >= 1); in mpn_perfect_square_p() 173 TRACE (gmp_printf ("mpn_perfect_square_p %Nd\n", up, usize)); in mpn_perfect_square_p() 195 up++, lo = up[0], usize--; in mpn_perfect_square_p() [all …]
|
| /dragonfly/contrib/mpfr/src/ |
| HD | eq.c | 32 mp_size_t usize, vsize, size, i; in mpfr_eq() local 59 usize = MPFR_LIMB_SIZE (u); in mpfr_eq() 62 if (vsize > usize) /* exchange u and v */ in mpfr_eq() 67 vsize = usize; in mpfr_eq() 68 usize = size; in mpfr_eq() 77 MPFR_ASSERTD(usize >= vsize); in mpfr_eq() 79 if (usize > vsize) in mpfr_eq() 86 k = usize - vsize - 1; in mpfr_eq() 104 size = usize; in mpfr_eq() 119 up += usize - size; in mpfr_eq()
|