Home
last modified time | relevance | path

Searched refs:GMP_LIMB_BITS (Results 1 – 25 of 52) sorted by relevance

123

/dragonfly/contrib/gmp/mpn/generic/
H A Dhgcd2.c101 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 Dperfpow.c85 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 Dget_d.c166 #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 Ddump.c39 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 Dperfsqr.c136 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 Dmod_1_2.c46 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 Dudiv_w_sdiv.c44 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 Dpowm_sec.c161 ((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 Dmod_1_3.c46 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 Dpow_1.c59 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 Dmod_1_4.c46 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 Ddivrem_2.c93 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 Dpowlo.c27 ((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 Dget_str.c149 #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 Dsqrtrem.c110 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 Dmod_1_1.c47 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 Dpopham.c74 #if GMP_LIMB_BITS > 32 in FNAME()
107 #if GMP_LIMB_BITS > 32 in FNAME()
H A Dpre_divrem_1.c109 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 Dremove.c30 #if GMP_LIMB_BITS > 50
33 #define LOG GMP_LIMB_BITS
/dragonfly/contrib/gmp/
H A Dextract-dbl.c44 #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 Ddumbmp.c65 #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 Dgmp-impl.h180 #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 Dpprime_p.c73 #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 Dfac_ui.c65 #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 Dvasnprintf.c336 # 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 …]

123