/netbsd/external/lgpl3/gmp/dist/ |
H A D | primesieve.c | 58 #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 D | extract-dbl.c | 52 #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 D | asl.h | 43 #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 D | hgcd2_jacobi.c | 110 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 D | hgcd2.c | 112 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 D | perfpow.c | 99 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 D | get_d.c | 175 #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 D | dump.c | 49 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 D | gcd_22.c | 45 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 D | perfsqr.c | 146 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 D | div_qr_2u_pi1.c | 56 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 D | udiv_w_sdiv.c | 54 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 D | bsqrtinv.c | 58 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 D | pow_1.c | 69 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 D | t-invert.c | 27 #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 D | t-aorsmul.c | 30 #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 D | modlinv.c | 43 #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 D | t-set_f.c | 30 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 D | t-cong_2exp.c | 123 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 D | mini-gmp.c | 658 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 D | pprime_p.c | 85 #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 D | primorial_ui.c | 66 __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 D | primes.c | 65 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 D | version.c | 212 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 D | invert_limb.c | 86 #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()
|