/dragonfly/contrib/mpfr/src/ |
H A D | get_str.c | 251 const mp_limb_t mpfr_l2b_2_0__tab[] = { 0x80000000 }; 281 const mp_limb_t mpfr_l2b_3_0__tab[] = { 0xcae00e00 }; 311 const mp_limb_t mpfr_l2b_4_0__tab[] = { 0x80000000 }; 2245 mp_limb_t *a; in mpfr_get_str() 2247 mp_limb_t *result; in mpfr_get_str() 2248 mp_limb_t *xp; in mpfr_get_str() 2249 mp_limb_t *reste; in mpfr_get_str() 2336 mp_limb_t *x1; in mpfr_get_str() 2340 count_leading_zeros (pow2, (mp_limb_t) b); in mpfr_get_str() 2438 mp_limb_t *x1; in mpfr_get_str() [all …]
|
H A D | mpfr-gmp.h | 77 #define MP_LIMB_T_MAX (~(mp_limb_t)0) 194 typedef mp_limb_t UWtype; 205 #define MP_BASE_AS_DOUBLE (4.0 * ((mp_limb_t) 1 << (GMP_NUMB_BITS - 2))) 291 __MPFR_DECLSPEC mp_size_t __gmpn_rootrem _MPFR_PROTO ((mp_limb_t*, 292 mp_limb_t*, mp_limb_t*, mp_size_t, mp_limb_t)); 298 __MPFR_DECLSPEC mp_limb_t __gmpn_sbpi1_divappr_q _MPFR_PROTO ((mp_limb_t*, 299 mp_limb_t*, mp_size_t, mp_limb_t*, mp_size_t, mp_limb_t)); 329 mp_limb_t dummy; \ 331 udiv_qrnnd (invxl, dummy, ~(xl), ~(mp_limb_t)0, xl); \ 343 mp_limb_t _v, _p, _t1, _t0, _mask; \ [all …]
|
H A D | set_d.c | 36 mp_limb_t manl; in __gmpfr_extract_double() 38 mp_limb_t manh; in __gmpfr_extract_double() 63 | ((mp_limb_t) x.s.manh << 43) | ((mp_limb_t) x.s.manl << 11)); in __gmpfr_extract_double() 72 manl = ((mp_limb_t) x.s.manh << 43) | ((mp_limb_t) x.s.manl << 11); in __gmpfr_extract_double() 123 manh = (mp_limb_t) d; in __gmpfr_extract_double() 124 manl = (mp_limb_t) ((d - manh) * MP_BASE_AS_DOUBLE); in __gmpfr_extract_double() 149 mp_limb_t tmpmant[MPFR_LIMBS_PER_DOUBLE]; in mpfr_set_d()
|
/dragonfly/contrib/gmp/mpn/generic/ |
H A D | hgcd2.c | 37 mp_limb_t n0, in div1() 38 mp_limb_t d0) in div1() 40 mp_limb_t q = 0; in div1() 91 mp_limb_t nh, mp_limb_t nl, in div2() 92 mp_limb_t dh, mp_limb_t dl) in div2() 94 mp_limb_t q = 0; in div2() 152 mp_limb_t nh, mp_limb_t nl, 153 mp_limb_t dh, mp_limb_t dl) 216 mpn_hgcd2 (mp_limb_t ah, mp_limb_t al, mp_limb_t bh, mp_limb_t bl, in mpn_hgcd2() 356 mp_limb_t r; in mpn_hgcd2() [all …]
|
H A D | gcdext_1.c | 48 mp_limb_t 50 mp_limb_t u, mp_limb_t v) in mpn_gcdext_1() 60 mp_limb_t s0 = 1; in mpn_gcdext_1() 61 mp_limb_t t0 = 0; in mpn_gcdext_1() 62 mp_limb_t s1 = 0; in mpn_gcdext_1() 64 mp_limb_t ug; in mpn_gcdext_1() 65 mp_limb_t vg; in mpn_gcdext_1() 66 mp_limb_t ugh; in mpn_gcdext_1() 67 mp_limb_t vgh; in mpn_gcdext_1() 266 mp_limb_t [all …]
|
H A D | mod_1.c | 72 static mp_limb_t 77 mp_limb_t dummy; in mpn_mod_1_unnorm() 136 mp_limb_t inv; in mpn_mod_1_unnorm() 154 static mp_limb_t 158 mp_limb_t n0, r; in mpn_mod_1_norm() 159 mp_limb_t dummy; in mpn_mod_1_norm() 201 mp_limb_t 220 mp_limb_t pre[4]; in mpn_mod_1() 233 mp_limb_t pre[4]; in mpn_mod_1() 239 mp_limb_t pre[5]; in mpn_mod_1() [all …]
|
H A D | mod_1_1.c | 33 mpn_mod_1_1p_cps (mp_limb_t cps[4], mp_limb_t b) in mpn_mod_1_1p_cps() 35 mp_limb_t bi; in mpn_mod_1_1p_cps() 36 mp_limb_t B1modb, B2modb; in mpn_mod_1_1p_cps() 57 mp_limb_t 58 mpn_mod_1_1p (mp_srcptr ap, mp_size_t n, mp_limb_t b, mp_limb_t bmodb[4]) in mpn_mod_1_1p() 60 mp_limb_t rh, rl, bi, q, ph, pl, r; in mpn_mod_1_1p() 61 mp_limb_t B1modb, B2modb; in mpn_mod_1_1p() 64 mp_limb_t mask; in mpn_mod_1_1p() 93 mask = -(mp_limb_t) (rh >= b); in mpn_mod_1_1p()
|
H A D | powm_sec.c | 74 mp_limb_t ul, lpl; \ 127 mp_limb_t ul, lpl; in mpn_local_sqr() 134 mp_limb_t cy; in mpn_local_sqr() 140 mp_limb_t cy; in mpn_local_sqr() 147 mp_limb_t cy; in mpn_local_sqr() 163 static inline mp_limb_t 167 mp_limb_t r; in getbits() 219 mp_limb_t minv; in mpn_powm_sec() 223 mp_limb_t expbits; in mpn_powm_sec() 309 mpn_tabselect (volatile mp_limb_t *rp, volatile mp_limb_t *tab, mp_size_t n, in mpn_tabselect() [all …]
|
H A D | mode1o.c | 103 mp_limb_t 104 mpn_modexact_1c_odd (mp_srcptr src, mp_size_t size, mp_limb_t d, in mpn_modexact_1c_odd() 105 mp_limb_t orig_c) in mpn_modexact_1c_odd() 107 mp_limb_t s, h, l, inverse, dummy, dmul, ret; in mpn_modexact_1c_odd() 108 mp_limb_t c = orig_c; in mpn_modexact_1c_odd() 194 mp_limb_t 195 mpn_modexact_1c_odd (mp_srcptr src, mp_size_t size, mp_limb_t d, mp_limb_t h) 197 mp_limb_t s, x, y, inverse, dummy, dmul, c1, c2; 198 mp_limb_t c = 0;
|
H A D | diveby3.c | 25 mp_limb_t 26 mpn_divexact_by3c (mp_ptr rp, mp_srcptr up, mp_size_t un, mp_limb_t c) in mpn_divexact_by3c() 28 mp_limb_t r; in mpn_divexact_by3c() 63 mp_limb_t 64 mpn_divexact_by3c (mp_ptr restrict rp, mp_srcptr restrict up, mp_size_t un, mp_limb_t c) in mpn_divexact_by3c() 66 mp_limb_t l, q, s; in mpn_divexact_by3c() 132 mp_limb_t 133 mpn_divexact_by3c (mp_ptr restrict rp, mp_srcptr restrict up, mp_size_t un, mp_limb_t cy) in mpn_divexact_by3c() 135 mp_limb_t s, sm, cl, q, qx, c2, c3; in mpn_divexact_by3c()
|
H A D | trialdiv.c | 44 mp_limb_t binv; 45 mp_limb_t lim; 49 mp_limb_t ppp; /* primes, multiplied together */ 50 mp_limb_t cps[7]; /* ppp values pre-computed for mpn_mod_1s_4p */ 67 mp_limb_t 70 mp_limb_t ppp; in mpn_trialdiv() 71 mp_limb_t *cps; in mpn_trialdiv() 74 mp_limb_t r, q; in mpn_trialdiv()
|
H A D | mod_1_3.c | 33 mpn_mod_1s_3p_cps (mp_limb_t cps[6], mp_limb_t b) in mpn_mod_1s_3p_cps() 35 mp_limb_t bi; in mpn_mod_1s_3p_cps() 36 mp_limb_t B1modb, B2modb, B3modb, B4modb; in mpn_mod_1s_3p_cps() 39 ASSERT (b <= (~(mp_limb_t) 0) / 3); in mpn_mod_1s_3p_cps() 72 mp_limb_t 73 mpn_mod_1s_3p (mp_srcptr ap, mp_size_t n, mp_limb_t b, mp_limb_t cps[6]) in mpn_mod_1s_3p() 75 mp_limb_t rh, rl, bi, q, ph, pl, ch, cl, r; in mpn_mod_1s_3p() 76 mp_limb_t B1modb, B2modb, B3modb, B4modb; in mpn_mod_1s_3p() 90 switch ((int) ((mp_limb_t) n * MODLIMB_INVERSE_3 >> (GMP_NUMB_BITS - 2))) in mpn_mod_1s_3p()
|
H A D | subcnd_n.c | 30 mp_limb_t 31 mpn_subcnd_n (mp_ptr rp, mp_srcptr up, mp_srcptr vp, mp_size_t n, mp_limb_t cnd) in mpn_subcnd_n() 33 mp_limb_t ul, vl, sl, rl, cy, cy1, cy2, mask; in mpn_subcnd_n() 39 mask = -(mp_limb_t) (cnd != 0); in mpn_subcnd_n() 61 mp_limb_t 62 mpn_subcnd_n (mp_ptr rp, mp_srcptr up, mp_srcptr vp, mp_size_t n, mp_limb_t cnd) in mpn_subcnd_n() 64 mp_limb_t ul, vl, rl, cy, mask; in mpn_subcnd_n() 70 mask = -(mp_limb_t) (cnd != 0); in mpn_subcnd_n()
|
H A D | perfpow.c | 39 mp_limb_t k, mp_bitcnt_t f, in pow_equals() 42 mp_limb_t *tp2; in pow_equals() 46 mp_limb_t h, l; in pow_equals() 81 mp_limb_t size; in pow_equals() 111 mp_limb_t k, mp_size_t bn, in binv_root() 203 mp_limb_t k, mp_srcptr yp, in is_kth_power() 207 mp_limb_t x, c; in is_kth_power() 255 mp_limb_t ub, mp_limb_t g, in perfpow() 258 mp_limb_t *yp, *tp, k = 0, *rp1; in perfpow() 326 mp_limb_t *nc, factor, g = 0; in mpn_perfect_power_p() [all …]
|
H A D | mod_1_2.c | 33 mpn_mod_1s_2p_cps (mp_limb_t cps[5], mp_limb_t b) in mpn_mod_1s_2p_cps() 35 mp_limb_t bi; in mpn_mod_1s_2p_cps() 36 mp_limb_t B1modb, B2modb, B3modb; in mpn_mod_1s_2p_cps() 39 ASSERT (b <= (~(mp_limb_t) 0) / 2); in mpn_mod_1s_2p_cps() 70 mp_limb_t 71 mpn_mod_1s_2p (mp_srcptr ap, mp_size_t n, mp_limb_t b, mp_limb_t cps[5]) in mpn_mod_1s_2p() 73 mp_limb_t rh, rl, bi, q, ph, pl, ch, cl, r; in mpn_mod_1s_2p() 74 mp_limb_t B1modb, B2modb, B3modb; in mpn_mod_1s_2p()
|
H A D | sqr_basecase.c | 40 mp_limb_t ul, lpl; \ 59 mp_limb_t cy; in mpn_sqr_basecase() 68 mp_limb_t ul, lpl; in mpn_sqr_basecase() 141 mp_limb_t cy; in mpn_sqr_basecase() 148 mp_limb_t x0, x1; in mpn_sqr_basecase() 152 mp_limb_t ul, lpl; in mpn_sqr_basecase() 193 mp_limb_t x0, x1; in mpn_sqr_basecase() 264 mp_limb_t ul, lpl; in mpn_sqr_basecase() 273 mp_limb_t cy; in mpn_sqr_basecase() 282 mp_limb_t cy; in mpn_sqr_basecase() [all …]
|
H A D | udiv_w_sdiv.c | 33 mp_limb_t 34 mpn_udiv_w_sdiv (mp_limb_t *rp, mp_limb_t a1, mp_limb_t a0, mp_limb_t d) in mpn_udiv_w_sdiv() 36 mp_limb_t q, r; in mpn_udiv_w_sdiv() 37 mp_limb_t c0, c1, b1; in mpn_udiv_w_sdiv() 56 q += (mp_limb_t) 1 << (GMP_LIMB_BITS - 1); in mpn_udiv_w_sdiv()
|
H A D | addmul_1.c | 31 mp_limb_t 32 mpn_addmul_1 (mp_ptr rp, mp_srcptr up, mp_size_t n, mp_limb_t vl) in mpn_addmul_1() 34 mp_limb_t ul, cl, hpl, lpl, rl; in mpn_addmul_1() 62 mp_limb_t 63 mpn_addmul_1 (mp_ptr rp, mp_srcptr up, mp_size_t n, mp_limb_t vl) in mpn_addmul_1() 65 mp_limb_t shifted_vl, ul, rl, lpl, hpl, prev_hpl, cl, xl, c1, c2, c3; in mpn_addmul_1() 98 mp_limb_t 99 mpn_addmul_1 (mp_ptr rp, mp_srcptr up, mp_size_t n, mp_limb_t vl) in mpn_addmul_1() 101 mp_limb_t shifted_vl, ul, rl, lpl, hpl, prev_hpl, xw, cl, xl; in mpn_addmul_1()
|
H A D | submul_1.c | 31 mp_limb_t 32 mpn_submul_1 (mp_ptr rp, mp_srcptr up, mp_size_t n, mp_limb_t vl) in mpn_submul_1() 34 mp_limb_t ul, cl, hpl, lpl, rl; in mpn_submul_1() 62 mp_limb_t 63 mpn_submul_1 (mp_ptr rp, mp_srcptr up, mp_size_t n, mp_limb_t vl) in mpn_submul_1() 65 mp_limb_t shifted_vl, ul, rl, lpl, hpl, prev_hpl, cl, xl, c1, c2, c3; in mpn_submul_1() 98 mp_limb_t 99 mpn_submul_1 (mp_ptr rp, mp_srcptr up, mp_size_t n, mp_limb_t vl) in mpn_submul_1() 101 mp_limb_t shifted_vl, ul, rl, lpl, hpl, prev_hpl, xw, cl, xl; in mpn_submul_1()
|
H A D | mod_1_4.c | 33 mpn_mod_1s_4p_cps (mp_limb_t cps[7], mp_limb_t b) in mpn_mod_1s_4p_cps() 35 mp_limb_t bi; in mpn_mod_1s_4p_cps() 36 mp_limb_t B1modb, B2modb, B3modb, B4modb, B5modb; in mpn_mod_1s_4p_cps() 39 ASSERT (b <= (~(mp_limb_t) 0) / 4); in mpn_mod_1s_4p_cps() 74 mp_limb_t 75 mpn_mod_1s_4p (mp_srcptr ap, mp_size_t n, mp_limb_t b, mp_limb_t cps[7]) in mpn_mod_1s_4p() 77 mp_limb_t rh, rl, bi, q, ph, pl, ch, cl, r; in mpn_mod_1s_4p() 78 mp_limb_t B1modb, B2modb, B3modb, B4modb, B5modb; in mpn_mod_1s_4p()
|
H A D | sbpi1_bdiv_qr.c | 43 mp_limb_t 46 mp_srcptr dp, mp_size_t dn, mp_limb_t dinv) in mpn_sbpi1_bdiv_qr() 50 mp_limb_t rh; in mpn_sbpi1_bdiv_qr() 51 mp_limb_t ql; in mpn_sbpi1_bdiv_qr() 67 mp_limb_t q; in mpn_sbpi1_bdiv_qr() 83 mp_limb_t q; in mpn_sbpi1_bdiv_qr() 102 mp_limb_t cy; in mpn_sbpi1_bdiv_qr()
|
/dragonfly/contrib/gmp/ |
H A D | gmp-impl.h | 82 __GMP_DECLSPEC mp_limb_t name __GMP_PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t)) 88 …__GMP_DECLSPEC mp_limb_t name __GMP_PROTO ((mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t, mp… 90 __GMP_DECLSPEC mp_limb_t name __GMP_PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t)) 310 mp_limb_t l; 774 __GMP_DECLSPEC mp_limb_t mpn_addmul_1c __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_li… 885 __GMP_DECLSPEC mp_limb_t mpn_mul_1c __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_limb_… 916 __GMP_DECLSPEC mp_limb_t mpn_submul_1c __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_li… 1299 __GMP_DECLSPEC mp_limb_t mpn_bdiv_dbm1c __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_l… 3766 __GMP_DECLSPEC int mpn_hgcd2 __GMP_PROTO ((mp_limb_t, mp_limb_t, mp_limb_t, mp_limb_t, struct hgcd_… 4178 mp_limb_t [all …]
|
H A D | dumbmp.c | 50 mp_limb_t *_mp_d; 90 mp_limb_t * 93 return (mp_limb_t *) xmalloc (n * sizeof (mp_limb_t)); in xmalloc_limbs() 159 PTR(r) = (mp_limb_t *) realloc (PTR(r), n * sizeof (mp_limb_t)); in mpz_realloc() 177 mpn_copyi (mp_limb_t *dst, mp_limb_t *src, int n) in mpn_copyi() 240 mp_limb_t hi; in mpz_sizeinbase() 292 mp_limb_t *rp; in mpz_setbit() 328 popc_limb (mp_limb_t a) in popc_limb() 360 mp_limb_t t, cy; in mpz_add() 413 mp_limb_t t, cy; in mpz_sub() [all …]
|
/dragonfly/contrib/gmp/mpz/ |
H A D | fac_ui.c | 29 static void ap_product_small __GMP_PROTO ((mpz_t, mp_limb_t, mp_limb_t, unsigned long, unsigned lon… 51 mp_limb_t __c; \ 112 mp_limb_t d[4]; in mpz_fac_ui() 133 umul_ppmm (d[1], d[0], (mp_limb_t) n, (mp_limb_t) FAC2OVERE); in mpz_fac_ui() 137 umul_ppmm (d[1], d[0], (mp_limb_t) n, (mp_limb_t) z); in mpz_fac_ui() 190 count_leading_zeros (stt, (mp_limb_t) n); in mpz_fac_ui() 235 popc_limb (i, (mp_limb_t) n); in mpz_fac_ui() 242 ap_product_small (mpz_t ret, mp_limb_t start, mp_limb_t step, in ap_product_small() 246 mp_limb_t b; in ap_product_small() 369 count_leading_zeros (n, (mp_limb_t) high); in odd_product() [all …]
|
H A D | aorsmul_i.c | 36 mp_limb_t __cy; \ 59 mpz_aorsmul_1 (mpz_ptr w, mpz_srcptr x, mp_limb_t y, mp_size_t sub) in mpz_aorsmul_1() 64 mp_limb_t cy; in mpz_aorsmul_1() 116 mp_limb_t cy2; in mpz_aorsmul_1() 158 mp_limb_t cy2; in mpz_aorsmul_1() 209 mpz_aorsmul_1 (w, t, (mp_limb_t) y >> GMP_NUMB_BITS, (mp_size_t) 0); in mpz_addmul_ui() 212 mpz_aorsmul_1 (w, t, (mp_limb_t) y & GMP_NUMB_MASK, (mp_size_t) 0); in mpz_addmul_ui() 217 mpz_aorsmul_1 (w, x, (mp_limb_t) y, (mp_size_t) 0); in mpz_addmul_ui() 237 mpz_aorsmul_1 (w, t, (mp_limb_t) y >> GMP_NUMB_BITS, (mp_size_t) -1); in mpz_submul_ui() 240 mpz_aorsmul_1 (w, t, (mp_limb_t) y & GMP_NUMB_MASK, (mp_size_t) -1); in mpz_submul_ui() [all …]
|