/dragonfly/contrib/gmp/mpn/generic/ |
H A D | hgcd2.c | 101 dh = (dh << 1) | (dl >> (GMP_LIMB_BITS - 1)); in div2() 113 dl = (dh << (GMP_LIMB_BITS - 1)) | (dl >> 1); in div2() 123 dh = (dh << 1) | (dl >> (GMP_LIMB_BITS - 1)); in div2() 129 dl = (dh << (GMP_LIMB_BITS - 1)) | (dl >> 1); in div2() 179 dl = (dh << (GMP_LIMB_BITS - 1)) | (dl >> 1); 252 if (ah < (CNST_LIMB(1) << (GMP_LIMB_BITS / 2))) in mpn_hgcd2() 254 ah = (ah << (GMP_LIMB_BITS / 2) ) + (al >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2() 255 bh = (bh << (GMP_LIMB_BITS / 2) ) + (bl >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2() 295 if (bh < (CNST_LIMB(1) << (GMP_LIMB_BITS / 2))) in mpn_hgcd2() 297 ah = (ah << (GMP_LIMB_BITS / 2) ) + (al >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2() [all …]
|
H A D | perfpow.c | 85 y = 2 + size / GMP_LIMB_BITS; in pow_equals() 115 mp_bitcnt_t order[GMP_LIMB_BITS * 2]; in binv_root() 140 if ((b % GMP_LIMB_BITS) != 0) in binv_root() 141 rp[(b - 1) / GMP_LIMB_BITS] &= (((mp_limb_t) 1) << (b % GMP_LIMB_BITS)) - 1; in binv_root() 179 bn = 1 + b / GMP_LIMB_BITS; in binv_sqroot() 192 rp[b / GMP_LIMB_BITS] &= (((mp_limb_t) 1) << (b % GMP_LIMB_BITS)) - 1; in binv_sqroot() 218 rn = 1 + b / GMP_LIMB_BITS; in is_kth_power() 278 if (b % GMP_LIMB_BITS) in perfpow() 279 yp[(b - 1) / GMP_LIMB_BITS] &= (((mp_limb_t) 1) << (b % GMP_LIMB_BITS)) - 1; in perfpow() 354 s = twos / GMP_LIMB_BITS; in mpn_perfect_power_p() [all …]
|
H A D | get_d.c | 166 #if GMP_LIMB_BITS == 64 in mpn_get_d() 175 nbits = GMP_LIMB_BITS - lshift; in mpn_get_d() 197 #if GMP_LIMB_BITS == 32 in mpn_get_d() 211 nbits = GMP_LIMB_BITS - lshift - 21; in mpn_get_d() 217 nbits = GMP_LIMB_BITS - lshift; in mpn_get_d() 223 mlo = x << GMP_LIMB_BITS - nbits - 11; in mpn_get_d() 281 #if GMP_LIMB_BITS > 53 in mpn_get_d() 291 lshift = GMP_LIMB_BITS - rshift; in mpn_get_d() 326 rshift = GMP_LIMB_BITS - lshift; in mpn_get_d() 345 rshift = GMP_LIMB_BITS - lshift; in mpn_get_d() [all …]
|
H A D | dump.c | 39 if ((ptr[n] >> GMP_LIMB_BITS / 2) != 0) in mpn_dump() 41 printf ("%lX", (unsigned long) (ptr[n] >> GMP_LIMB_BITS / 2)); in mpn_dump() 42 printf ("%0*lX", (GMP_LIMB_BITS / 2 / 4), (unsigned long) ptr[n]); in mpn_dump() 52 printf ("%0*lX", (GMP_NUMB_BITS - GMP_LIMB_BITS / 2) / 4, in mpn_dump() 53 (unsigned long) (ptr[n] >> GMP_LIMB_BITS / 2)); in mpn_dump() 54 printf ("%0*lX", GMP_LIMB_BITS / 2 / 4, (unsigned long) ptr[n]); in mpn_dump()
|
H A D | perfsqr.c | 136 ASSERT ((d) <= GMP_LIMB_BITS); \ 153 ASSERT ((d) <= 2*GMP_LIMB_BITS); \ 158 m = ((int) idx - GMP_LIMB_BITS < 0 ? (mlo) : (mhi)); \ 159 idx %= GMP_LIMB_BITS; \ 179 if (((sq_res_0x100[idx / GMP_LIMB_BITS] in mpn_perfect_square_p() 180 >> (idx % GMP_LIMB_BITS)) & 1) == 0) in mpn_perfect_square_p()
|
H A D | mod_1_2.c | 46 B1modb = -b * ((bi >> (GMP_LIMB_BITS-cnt)) | (CNST_LIMB(1) << cnt)); in mpn_mod_1s_2p_cps() 91 udiv_qrnnd_preinv (q, r, rl >> (GMP_LIMB_BITS - cnt), in mpn_mod_1s_2p() 131 r = (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)); in mpn_mod_1s_2p() 133 udiv_qrnnd_preinv (q, r, rh >> (GMP_LIMB_BITS - cnt), in mpn_mod_1s_2p() 134 (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)), b, bi); in mpn_mod_1s_2p()
|
H A D | udiv_w_sdiv.c | 44 if (a1 < d - a1 - (a0 >> (GMP_LIMB_BITS - 1))) in mpn_udiv_w_sdiv() 52 sub_ddmmss (c1, c0, a1, a0, d >> 1, d << (GMP_LIMB_BITS - 1)); in mpn_udiv_w_sdiv() 56 q += (mp_limb_t) 1 << (GMP_LIMB_BITS - 1); in mpn_udiv_w_sdiv() 63 c0 = (a1 << (GMP_LIMB_BITS - 1)) + (a0 >> 1); in mpn_udiv_w_sdiv()
|
H A D | powm_sec.c | 161 ((p[(bi - 1) / GMP_LIMB_BITS] >> (bi - 1) % GMP_LIMB_BITS) & 1) 177 i = bi / GMP_LIMB_BITS; /* word index of low bit to extract */ in getbits() 178 bi %= GMP_LIMB_BITS; /* bit index in low word */ in getbits() 180 nbits_in_r = GMP_LIMB_BITS - bi; /* number of bits now in r */ in getbits() 232 ebi = (mp_bitcnt_t) en * GMP_LIMB_BITS - cnt; in mpn_powm_sec()
|
H A D | mod_1_3.c | 46 B1modb = -b * ((bi >> (GMP_LIMB_BITS-cnt)) | (CNST_LIMB(1) << cnt)); in mpn_mod_1s_3p_cps() 138 r = (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)); in mpn_mod_1s_3p() 140 udiv_qrnnd_preinv (q, r, rh >> (GMP_LIMB_BITS - cnt), in mpn_mod_1s_3p() 141 (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)), b, bi); in mpn_mod_1s_3p()
|
H A D | pow_1.c | 59 cnt = GMP_LIMB_BITS; in mpn_pow_1() 77 for (i = GMP_LIMB_BITS - cnt - 1;;) in mpn_pow_1() 102 for (i = GMP_LIMB_BITS - cnt - 1;;) in mpn_pow_1()
|
H A D | mod_1_4.c | 46 B1modb = -b * ((bi >> (GMP_LIMB_BITS-cnt)) | (CNST_LIMB(1) << cnt)); in mpn_mod_1s_4p_cps() 151 r = (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)); in mpn_mod_1s_4p() 153 udiv_qrnnd_preinv (q, r, rh >> (GMP_LIMB_BITS - cnt), in mpn_mod_1s_4p() 154 (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)), b, bi); in mpn_mod_1s_4p()
|
H A D | divrem_2.c | 93 n1 = n1 - d1 - (n0 >> GMP_LIMB_BITS - 1); in mpn_divrem_2() 158 n1 = n1 - (n0 >> GMP_LIMB_BITS - 1); in mpn_divrem_2() 171 n1 = r - n1 - (n0 >> GMP_LIMB_BITS - 1); in mpn_divrem_2()
|
H A D | powlo.c | 27 ((p[(bi - 1) / GMP_LIMB_BITS] >> (bi - 1) % GMP_LIMB_BITS) & 1) 85 ebi = (mp_bitcnt_t) en * GMP_LIMB_BITS - cnt; in mpn_powlo()
|
H A D | get_str.c | 149 #define BUF_ALLOC (GET_STR_THRESHOLD_LIMIT * GMP_LIMB_BITS * 7 / 11) in mpn_sb_get_str() 151 #define BUF_ALLOC (GET_STR_PRECOMPUTE_THRESHOLD * GMP_LIMB_BITS * 7 / 11) in mpn_sb_get_str() 220 digit = frac >> (GMP_LIMB_BITS - 4); in mpn_sb_get_str() 361 powers_t powtab[GMP_LIMB_BITS]; in mpn_get_str() 439 mp_size_t n_pows, xn, pn, exptab[GMP_LIMB_BITS], bexp; in mpn_get_str()
|
H A D | sqrtrem.c | 110 ASSERT_ALWAYS (GMP_LIMB_BITS == 32 || GMP_LIMB_BITS == 64); in mpn_sqrtrem1() 117 abits = a0 >> (GMP_LIMB_BITS - 1 - 8); /* extract bits for table lookup */ in mpn_sqrtrem1() 123 a1 = a0 >> (GMP_LIMB_BITS - 1 - 32); in mpn_sqrtrem1()
|
H A D | mod_1_1.c | 47 B1modb = -b * ((bi >> (GMP_LIMB_BITS-cnt)) | (CNST_LIMB(1) << cnt)); in mpn_mod_1_1p_cps() 91 rh = (rh << cnt) | (rl >> (GMP_LIMB_BITS - cnt)); in mpn_mod_1_1p()
|
H A D | popham.c | 74 #if GMP_LIMB_BITS > 32 in FNAME() 107 #if GMP_LIMB_BITS > 32 in FNAME()
|
H A D | pre_divrem_1.c | 109 r |= n1 >> (GMP_LIMB_BITS - shift); in mpn_preinv_divrem_1() 116 ((n1 << shift) | (n0 >> (GMP_LIMB_BITS - shift))), in mpn_preinv_divrem_1()
|
H A D | remove.c | 30 #if GMP_LIMB_BITS > 50 33 #define LOG GMP_LIMB_BITS
|
/dragonfly/contrib/gmp/ |
H A D | extract-dbl.c | 44 #define BITS_PER_PART GMP_LIMB_BITS in __gmp_extract_double() 164 rp[1] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double() 176 rp[1] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double() 189 rp[1] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double() 199 rp[2] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double() 217 rp[2] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double() 228 rp[2] = manh >> (GMP_LIMB_BITS - sc); in __gmp_extract_double() 229 rp[1] = (manh << sc) | (manl >> (GMP_LIMB_BITS - sc)); in __gmp_extract_double() 242 rp[2] = (manh >> (GMP_LIMB_BITS - sc)); in __gmp_extract_double() 261 rp[2] = (manh >> (GMP_LIMB_BITS - sc)); in __gmp_extract_double() [all …]
|
H A D | dumbmp.c | 65 #define HI(x) ((x) >> GMP_LIMB_BITS) 298 limb = bit / GMP_LIMB_BITS; in mpz_setbit() 299 bit %= GMP_LIMB_BITS; in mpz_setbit() 319 limb = bit / GMP_LIMB_BITS; in mpz_tstbit() 323 bit %= GMP_LIMB_BITS; in mpz_tstbit() 556 lcnt = bcnt / GMP_LIMB_BITS; in mpz_tdiv_q_2exp() 569 cnt = bcnt % GMP_LIMB_BITS; in mpz_tdiv_q_2exp() 573 tnc = GMP_LIMB_BITS - cnt; in mpz_tdiv_q_2exp() 604 bwhole = bcnt / GMP_LIMB_BITS; in mpz_tdiv_r_2exp() 605 bcnt %= GMP_LIMB_BITS; in mpz_tdiv_r_2exp() [all …]
|
H A D | gmp-impl.h | 180 #ifndef GMP_LIMB_BITS 181 #define GMP_LIMB_BITS (8 * SIZEOF_MP_LIMB_T) macro 2253 (cout) = __w >> (GMP_LIMB_BITS-1); \ 3089 && GMP_LIMB_BITS == 32 3097 && defined (__amd64__) && GMP_LIMB_BITS == 64 3127 #if GMP_LIMB_BITS == 8 3131 #if GMP_LIMB_BITS == 16 3137 #if GMP_LIMB_BITS == 32 3147 #if GMP_LIMB_BITS == 64 3382 #define W_TYPE_SIZE GMP_LIMB_BITS [all …]
|
/dragonfly/contrib/gmp/mpz/ |
H A D | pprime_p.c | 73 #if GMP_LIMB_BITS >= 4 in mpz_probab_prime_p() 76 #if GMP_LIMB_BITS >= 8 in mpz_probab_prime_p() 79 #if GMP_LIMB_BITS >= 16 in mpz_probab_prime_p() 82 #if GMP_LIMB_BITS >= 32 in mpz_probab_prime_p() 85 #if GMP_LIMB_BITS >= 64 in mpz_probab_prime_p()
|
H A D | fac_ui.c | 65 #if BITS_PER_ULONG == GMP_LIMB_BITS 136 z = GMP_LIMB_BITS - z - 1; /* z=floor(log_2(2n/e)) */ in mpz_fac_ui() 139 d[0] = (d[0] >> 2) | (d[1] << (GMP_LIMB_BITS - 2)); in mpz_fac_ui() 191 stt = GMP_LIMB_BITS - stt + 1 - APCONST; in mpz_fac_ui() 201 z = GMP_LIMB_BITS - z; in mpz_fac_ui() 371 n = GMP_LIMB_BITS - n - APCONST; in odd_product()
|
/dragonfly/contrib/diffutils/lib/ |
H A D | vasnprintf.c | 336 # define GMP_LIMB_BITS 32 macro 921 m.nlimbs = (LDBL_MANT_BIT + GMP_LIMB_BITS - 1) / GMP_LIMB_BITS; in decode_long_double() 937 # if (LDBL_MANT_BIT % GMP_LIMB_BITS) > GMP_LIMB_BITS / 2 in decode_long_double() 950 m.limbs[LDBL_MANT_BIT / GMP_LIMB_BITS] = (hi << (GMP_LIMB_BITS / 2)) | lo; in decode_long_double() 1009 m.nlimbs = (DBL_MANT_BIT + GMP_LIMB_BITS - 1) / GMP_LIMB_BITS; in decode_double() 1025 # if (DBL_MANT_BIT % GMP_LIMB_BITS) > GMP_LIMB_BITS / 2 in decode_double() 1038 m.limbs[DBL_MANT_BIT / GMP_LIMB_BITS] = (hi << (GMP_LIMB_BITS / 2)) | lo; in decode_double() 1356 while (y < (1.0L / (1 << (GMP_LIMB_BITS / 2)) / (1 << (GMP_LIMB_BITS / 2)))) in floorlog10l() 1358 y *= 1.0L * (1 << (GMP_LIMB_BITS / 2)) * (1 << (GMP_LIMB_BITS / 2)); in floorlog10l() 1447 while (y < (1.0 / (1 << (GMP_LIMB_BITS / 2)) / (1 << (GMP_LIMB_BITS / 2)))) in floorlog10() [all …]
|