Home
last modified time | relevance | path

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

12345678910>>...15

/netbsd/external/lgpl3/gmp/dist/
H A Dprimesieve.c58 #if GMP_LIMB_BITS > 61
60 #if GMP_LIMB_BITS == 64
73 #if GMP_LIMB_BITS > 30
75 #if GMP_LIMB_BITS == 32
85 #if GMP_LIMB_BITS > 15
89 #if GMP_LIMB_BITS > 7
102 if (off < GMP_LIMB_BITS) { \
253 limbs = bits / GMP_LIMB_BITS; in first_block_primesieve()
265 ASSERT (i < GMP_LIMB_BITS); in first_block_primesieve()
325 ASSERT (i < GMP_LIMB_BITS); in block_resieve()
[all …]
H A Dextract-dbl.c52 #define BITS_PER_PART GMP_LIMB_BITS in __gmp_extract_double()
172 rp[1] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double()
184 rp[1] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double()
197 rp[1] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double()
207 rp[2] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double()
225 rp[2] = manl >> (GMP_LIMB_BITS - sc); in __gmp_extract_double()
236 rp[2] = manh >> (GMP_LIMB_BITS - sc); in __gmp_extract_double()
237 rp[1] = (manh << sc) | (manl >> (GMP_LIMB_BITS - sc)); in __gmp_extract_double()
250 rp[2] = (manh >> (GMP_LIMB_BITS - sc)); in __gmp_extract_double()
269 rp[2] = (manh >> (GMP_LIMB_BITS - sc)); in __gmp_extract_double()
[all …]
H A Dasl.h43 #ifndef GMP_LIMB_BITS
44 #define GMP_LIMB_BITS 4 macro
47 #define GMP_NUMB_MASK (2 * (1ul << (GMP_LIMB_BITS - 1)) - 1)
77 GMP_ASSERT_ALWAYS (rhs < GMP_LIMB_BITS);
82 GMP_ASSERT_ALWAYS (rhs < GMP_LIMB_BITS);
/netbsd/external/lgpl3/gmp/dist/mpn/generic/
H A Dhgcd2_jacobi.c110 dh = (dh << 1) | (dl >> (GMP_LIMB_BITS - 1)); in div2()
122 dl = (dh << (GMP_LIMB_BITS - 1)) | (dl >> 1); in div2()
132 dh = (dh << 1) | (dl >> (GMP_LIMB_BITS - 1)); in div2()
138 dl = (dh << (GMP_LIMB_BITS - 1)) | (dl >> 1); in div2()
196 if (ah < (CNST_LIMB(1) << (GMP_LIMB_BITS / 2))) in mpn_hgcd2_jacobi()
198 ah = (ah << (GMP_LIMB_BITS / 2) ) + (al >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2_jacobi()
199 bh = (bh << (GMP_LIMB_BITS / 2) ) + (bl >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2_jacobi()
242 if (bh < (CNST_LIMB(1) << (GMP_LIMB_BITS / 2))) in mpn_hgcd2_jacobi()
244 ah = (ah << (GMP_LIMB_BITS / 2) ) + (al >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2_jacobi()
245 bh = (bh << (GMP_LIMB_BITS / 2) ) + (bl >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2_jacobi()
[all …]
H A Dhgcd2.c112 if (UNLIKELY ((d0 >> (GMP_LIMB_BITS - 3)) != 0) in div1()
282 nbi = n0 << ncnt >> (GMP_LIMB_BITS - NBITS); in div1()
283 dbi = d0 << ncnt >> (GMP_LIMB_BITS - NBITS); in div1()
379 dbi = d0 << dcnt >> (GMP_LIMB_BITS - NBITS); in div1()
425 n2 = n1 >> (GMP_LIMB_BITS - c); in div2()
426 n1 = (n1 << c) | (n0 >> (GMP_LIMB_BITS - c)); in div2()
428 d1 = (d1 << c) | (d0 >> (GMP_LIMB_BITS - c)); in div2()
564 ah = (ah << (GMP_LIMB_BITS / 2) ) + (al >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2()
565 bh = (bh << (GMP_LIMB_BITS / 2) ) + (bl >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2()
607 ah = (ah << (GMP_LIMB_BITS / 2) ) + (al >> (GMP_LIMB_BITS / 2)); in mpn_hgcd2()
[all …]
H A Dperfpow.c99 y = 2 + size / GMP_LIMB_BITS; in pow_equals()
133 rn = 1 + b / GMP_LIMB_BITS; in is_kth_power()
136 rp[rn - 1] &= (CNST_LIMB(1) << (b % GMP_LIMB_BITS)) - 1; in is_kth_power()
153 rn = 1 + (b - 1) / GMP_LIMB_BITS; in is_kth_power()
155 if ((b % GMP_LIMB_BITS) != 0) in is_kth_power()
156 rp[rn - 1] &= (CNST_LIMB(1) << (b % GMP_LIMB_BITS)) - 1; in is_kth_power()
195 mpn_binvert (ip, np, 1 + (b - 1) / GMP_LIMB_BITS, tp); in perfpow()
196 if (b % GMP_LIMB_BITS) in perfpow()
197 ip[(b - 1) / GMP_LIMB_BITS] &= (CNST_LIMB(1) << (b % GMP_LIMB_BITS)) - 1; in perfpow()
272 s = twos / GMP_LIMB_BITS; in mpn_perfect_power_p()
[all …]
H A Dget_d.c175 #if GMP_LIMB_BITS == 64 in mpn_get_d()
182 nbits = GMP_LIMB_BITS - lshift; in mpn_get_d()
204 #if GMP_LIMB_BITS == 32 in mpn_get_d()
217 nbits = GMP_LIMB_BITS - lshift - 21; in mpn_get_d()
223 nbits = GMP_LIMB_BITS - lshift; in mpn_get_d()
229 mlo = x << (GMP_LIMB_BITS - nbits - 11); in mpn_get_d()
288 #if GMP_LIMB_BITS > 53 in mpn_get_d()
298 lshift = GMP_LIMB_BITS - rshift; in mpn_get_d()
330 mhi += mlo >> (GMP_LIMB_BITS - lshift); in mpn_get_d()
336 mlo += x >> (GMP_LIMB_BITS - lshift); in mpn_get_d()
[all …]
H A Ddump.c49 if ((ptr[n] >> GMP_LIMB_BITS / 2) != 0) in mpn_dump()
51 printf ("%lX", (unsigned long) (ptr[n] >> GMP_LIMB_BITS / 2)); in mpn_dump()
52 printf ("%0*lX", (GMP_LIMB_BITS / 2 / 4), (unsigned long) ptr[n]); in mpn_dump()
62 printf ("%0*lX", (GMP_NUMB_BITS - GMP_LIMB_BITS / 2) / 4, in mpn_dump()
63 (unsigned long) (ptr[n] >> GMP_LIMB_BITS / 2)); in mpn_dump()
64 printf ("%0*lX", GMP_LIMB_BITS / 2 / 4, (unsigned long) ptr[n]); in mpn_dump()
H A Dgcd_22.c45 u0 = (u0 >> 1) | (u1 << (GMP_LIMB_BITS - 1)); in mpn_gcd_22()
48 v0 = (v0 >> 1) | (v1 << (GMP_LIMB_BITS - 1)); in mpn_gcd_22()
62 g.d1 = (u1 << 1) | (u0 >> (GMP_LIMB_BITS - 1)); in mpn_gcd_22()
72 ASSERT (c < GMP_LIMB_BITS - 1); in mpn_gcd_22()
94 if (UNLIKELY (c == GMP_LIMB_BITS)) in mpn_gcd_22()
101 u0 = (u0 >> c) | (u1 << (GMP_LIMB_BITS - c)); in mpn_gcd_22()
113 g.d1 = u0 >> (GMP_LIMB_BITS - 1); in mpn_gcd_22()
H A Dperfsqr.c146 ASSERT ((d) <= GMP_LIMB_BITS); \
163 ASSERT ((d) <= 2*GMP_LIMB_BITS); \
168 m = ((int) idx - GMP_LIMB_BITS < 0 ? (mlo) : (mhi)); \
169 idx %= GMP_LIMB_BITS; \
189 if (((sq_res_0x100[idx / GMP_LIMB_BITS] in mpn_perfect_square_p()
190 >> (idx % GMP_LIMB_BITS)) & 1) == 0) in mpn_perfect_square_p()
H A Ddiv_qr_2u_pi1.c56 r2 = np[nn-1] >> (GMP_LIMB_BITS - shift); in mpn_div_qr_2u_pi1()
57 r1 = (np[nn-1] << shift) | (np[nn-2] >> (GMP_LIMB_BITS - shift)); in mpn_div_qr_2u_pi1()
66 r1 |= r0 >> (GMP_LIMB_BITS - shift); in mpn_div_qr_2u_pi1()
72 rp[0] = (r1 >> shift) | (r2 << (GMP_LIMB_BITS - shift)); in mpn_div_qr_2u_pi1()
H A Dudiv_w_sdiv.c54 if (a1 < d - a1 - (a0 >> (GMP_LIMB_BITS - 1))) in mpn_udiv_w_sdiv()
62 sub_ddmmss (c1, c0, a1, a0, d >> 1, d << (GMP_LIMB_BITS - 1)); in mpn_udiv_w_sdiv()
66 q += (mp_limb_t) 1 << (GMP_LIMB_BITS - 1); in mpn_udiv_w_sdiv()
73 c0 = (a1 << (GMP_LIMB_BITS - 1)) + (a0 >> 1); in mpn_udiv_w_sdiv()
H A Dbsqrtinv.c58 mp_size_t bn, order[GMP_LIMB_BITS + 1]; in mpn_bsqrtinv()
63 bn = 1 + bnb / GMP_LIMB_BITS; in mpn_bsqrtinv()
85 bn = 1 + bnb / GMP_LIMB_BITS; in mpn_bsqrtinv()
H A Dpow_1.c69 cnt = GMP_LIMB_BITS; in mpn_pow_1()
91 for (i = GMP_LIMB_BITS - cnt - 1;;) in mpn_pow_1()
116 for (i = GMP_LIMB_BITS - cnt - 1;;) in mpn_pow_1()
/netbsd/external/lgpl3/gmp/dist/mini-gmp/tests/
H A Dt-invert.c27 #define GMP_LIMB_BITS (sizeof(mp_limb_t) * CHAR_BIT) macro
43 mpz_setbit (m, GMP_LIMB_BITS); in test_2by1()
48 mpz_setbit (t, 2* GMP_LIMB_BITS); in test_2by1()
79 mpz_setbit (m, GMP_LIMB_BITS); in test_3by2()
84 mpz_setbit (t, 3 * GMP_LIMB_BITS); in test_3by2()
114 if (GMP_LIMB_BITS == 64) in testmain()
124 mini_urandomb (u, GMP_LIMB_BITS); in testmain()
125 mpz_setbit (u, GMP_LIMB_BITS -1); in testmain()
132 mini_urandomb (u, 2*GMP_LIMB_BITS); in testmain()
133 mpz_setbit (u, 2*GMP_LIMB_BITS -1); in testmain()
H A Dt-aorsmul.c30 #define GMP_LIMB_BITS (sizeof(mp_limb_t) * CHAR_BIT) macro
31 #define MAXLIMBS ((MAXBITS + GMP_LIMB_BITS - 1) / GMP_LIMB_BITS)
/netbsd/external/lgpl3/gmp/dist/tune/
H A Dmodlinv.c43 #if GMP_LIMB_BITS <= 32
57 #if GMP_LIMB_BITS > 32 && GMP_LIMB_BITS <= 64
113 __count = GMP_LIMB_BITS-1; \
144 __count = GMP_LIMB_BITS-1; \
148 __inv = (__inv >> 1) | (__lowbit << (GMP_LIMB_BITS-1)); \
/netbsd/external/lgpl3/gmp/dist/tests/mpz/
H A Dt-set_f.c30 0, 1, GMP_LIMB_BITS, 2*GMP_LIMB_BITS, 5*GMP_LIMB_BITS in check_one()
103 mpz_rrandomb (z, RANDS, 2*GMP_LIMB_BITS); in check_various()
106 mpz_rrandomb (z, RANDS, 5*GMP_LIMB_BITS); in check_various()
H A Dt-cong_2exp.c123 mpz_errandomb (a, rands, 8*GMP_LIMB_BITS); in check_random()
124 mpz_errandomb (c, rands, 8*GMP_LIMB_BITS); in check_random()
125 d = urandom() % (8*GMP_LIMB_BITS); in check_random()
127 mpz_mul_2exp (a, a, urandom() % (2*GMP_LIMB_BITS)); in check_random()
128 mpz_mul_2exp (c, c, urandom() % (2*GMP_LIMB_BITS)); in check_random()
160 mp_bitcnt_t m = 10 * GMP_LIMB_BITS; in check_random_bits()
/netbsd/external/lgpl3/gmp/dist/mini-gmp/
H A Dmini-gmp.c658 tnc = GMP_LIMB_BITS - cnt; in mpn_lshift()
685 tnc = GMP_LIMB_BITS - cnt; in mpn_rshift()
725 i = bit / GMP_LIMB_BITS; in mpn_scan1()
735 i = bit / GMP_LIMB_BITS; in mpn_scan0()
1221 shift -= GMP_LIMB_BITS; in mpn_get_str_bits()
1335 if (shift >= GMP_LIMB_BITS) in mpn_set_str_bits()
1337 shift -= GMP_LIMB_BITS; in mpn_set_str_bits()
1748 m -= GMP_LIMB_BITS; in mpz_get_d()
2340 bit_index %= GMP_LIMB_BITS; in mpz_div_q_2exp()
2391 rn = (bit_index + GMP_LIMB_BITS - 1) / GMP_LIMB_BITS; in mpz_div_r_2exp()
[all …]
/netbsd/external/lgpl3/gmp/dist/mpz/
H A Dpprime_p.c85 #if GMP_LIMB_BITS >= 4 in mpz_probab_prime_p()
88 #if GMP_LIMB_BITS >= 8 in mpz_probab_prime_p()
91 #if GMP_LIMB_BITS >= 16 in mpz_probab_prime_p()
94 #if GMP_LIMB_BITS >= 32 in mpz_probab_prime_p()
97 #if GMP_LIMB_BITS >= 64 in mpz_probab_prime_p()
H A Dprimorial_ui.c66 __index = __i / GMP_LIMB_BITS; \
67 __mask = CNST_LIMB(1) << (__i % GMP_LIMB_BITS); \
74 __mask = __mask << 1 | __mask >> (GMP_LIMB_BITS-1); \
101 primesieve_size (mp_limb_t n) { return n_to_bit(n) / GMP_LIMB_BITS + 1; } in primesieve_size()
/netbsd/external/lgpl3/gmp/dist/tests/devel/
H A Dprimes.c65 primesieve_size (mp_limb_t n) { return n_to_bit(n) / GMP_LIMB_BITS + 1; } in primesieve_size()
86 __index = __i / GMP_LIMB_BITS; \
87 __mask = CNST_LIMB(1) << (__i % GMP_LIMB_BITS); \
94 __mask = __mask << 1 | __mask >> (GMP_LIMB_BITS-1); \
137 start = mpz_get_ui (g) + GMP_LIMB_BITS; in check_pprime()
150 …TRACE (printf ("LOOP =%li - %li\n", id_to_n (off+1), id_to_n (off + BLOCK_SIZE * GMP_LIMB_BITS)),3… in check_pprime()
151 LOOP_ON_SIEVE_BEGIN (prime, off, off + BLOCK_SIZE * GMP_LIMB_BITS - 1, in check_pprime()
174 off += BLOCK_SIZE * GMP_LIMB_BITS; in check_pprime()
/netbsd/external/lgpl3/mpfr/dist/examples/
H A Dversion.c212 printf ("GMP_LIMB_BITS = %d\n", (int) GMP_LIMB_BITS); in main()
217 if (mp_bits_per_limb != GMP_LIMB_BITS) in main()
219 if (GMP_LIMB_BITS != sizeof(mp_limb_t) * CHAR_BIT) in main()
224 if (c != GMP_LIMB_BITS) in main()
/netbsd/external/lgpl3/gmp/dist/mpn/cray/ieee/
H A Dinvert_limb.c86 #if GMP_LIMB_BITS == 32 in mpn_invert_limb()
93 #if GMP_LIMB_BITS == 64 in mpn_invert_limb()
110 th = (th << 2) | (tl >> GMP_LIMB_BITS - 2); in mpn_invert_limb()

12345678910>>...15