Home
last modified time | relevance | path

Searched refs:mpfr_exp_t (Results 1 – 25 of 276) sorted by relevance

12345678910>>...12

/netbsd/src/external/lgpl3/mpfr/dist/src/
Dexp_2.c28 mpfr_exp2_aux (mpz_t, mpfr_srcptr, mpfr_prec_t, mpfr_exp_t *);
30 mpfr_exp2_aux2 (mpz_t, mpfr_srcptr, mpfr_prec_t, mpfr_exp_t *);
31 static mpfr_exp_t
32 mpz_normalize (mpz_t, mpz_t, mpfr_exp_t);
33 static mpfr_exp_t
34 mpz_normalize2 (mpz_t, mpz_t, mpfr_exp_t, mpfr_exp_t);
39 static mpfr_exp_t
40 mpz_normalize (mpz_t rop, mpz_t z, mpfr_exp_t q) in mpz_normalize()
49 return (mpfr_exp_t) k - q; in mpz_normalize()
59 static mpfr_exp_t
[all …]
Dstrtofr.c39 mpfr_exp_t exp_base; /* number of digits before the point, + exponent
41 mpfr_exp_t exp_bin; /* binary exponent of the pxxx format for
397 mpfr_exp_t sum; in parse_string()
404 (mpfr_exp_t) read_exp; in parse_string()
406 mpfr_exp_t, mpfr_uexp_t, in parse_string()
425 (mpfr_exp_t) read_exp; in parse_string()
461 mpfr_exp_t exp; in parsed_string_to_mpfr()
594 exp = - ((mpfr_exp_t) diff_ysize * GMP_NUMB_BITS + count); in parsed_string_to_mpfr()
635 mpfr_exp_t tmp; in parsed_string_to_mpfr()
645 MPFR_SADD_OVERFLOW (tmp, pstr->exp_base, -(mpfr_exp_t) pstr_size, in parsed_string_to_mpfr()
[all …]
Dexceptions.c26 MPFR_THREAD_VAR (mpfr_exp_t, __gmpfr_emin, MPFR_EMIN_DEFAULT) in MPFR_THREAD_VAR() argument
27 MPFR_THREAD_VAR (mpfr_exp_t, __gmpfr_emax, MPFR_EMAX_DEFAULT) in MPFR_THREAD_VAR()
31 MPFR_COLD_FUNCTION_ATTR mpfr_exp_t in MPFR_THREAD_VAR()
40 mpfr_set_emin (mpfr_exp_t exponent) in mpfr_set_emin()
53 MPFR_COLD_FUNCTION_ATTR mpfr_exp_t
59 MPFR_COLD_FUNCTION_ATTR mpfr_exp_t
67 MPFR_COLD_FUNCTION_ATTR mpfr_exp_t
76 mpfr_set_emax (mpfr_exp_t exponent) in mpfr_set_emax()
89 MPFR_COLD_FUNCTION_ATTR mpfr_exp_t
95 MPFR_COLD_FUNCTION_ATTR mpfr_exp_t
[all …]
Dexp.c28 static MPFR_THREAD_ATTR mpfr_exp_t previous_emin;
31 static MPFR_THREAD_ATTR mpfr_exp_t previous_emax;
38 mpfr_exp_t expx; in mpfr_exp()
82 MPFR_TMP_INIT1(e_limb, e, sizeof (mpfr_exp_t) * CHAR_BIT); in mpfr_exp()
111 MPFR_TMP_INIT1(e_limb, e, sizeof (mpfr_exp_t) * CHAR_BIT - 1); in mpfr_exp()
137 mpfr_exp_t emin = __gmpfr_emin; in mpfr_exp()
138 mpfr_exp_t emax = __gmpfr_emax; in mpfr_exp()
Dsin.c41 mpfr_exp_t expx, err1, err; in mpfr_sin()
144 || MPFR_GET_EXP(xr) < (mpfr_exp_t) 3 - (mpfr_exp_t) m in mpfr_sin()
146 || MPFR_GET_EXP(c) < (mpfr_exp_t) 3 - (mpfr_exp_t) m) in mpfr_sin()
181 err = 2 * MPFR_GET_EXP (c) + (mpfr_exp_t) m - 3 - (reduce != 0); in mpfr_sin()
186 if (err < (mpfr_exp_t) MPFR_PREC (y)) in mpfr_sin()
Dmpfr.h186 typedef short mpfr_exp_t; typedef
189 typedef int mpfr_exp_t; typedef
192 typedef long mpfr_exp_t; typedef
209 typedef intmax_t mpfr_exp_t; typedef
216 #define MPFR_EMAX_DEFAULT ((mpfr_exp_t) (((mpfr_ulong) 1 << 30) - 1))
226 #define __MPFR_EXP_MAX ((mpfr_exp_t) (((mpfr_uexp_t) -1) >> 1))
235 mpfr_exp_t _mpfr_exp;
401 __MPFR_DECLSPEC mpfr_exp_t mpfr_get_emin (void);
402 __MPFR_DECLSPEC int mpfr_set_emin (mpfr_exp_t);
403 __MPFR_DECLSPEC mpfr_exp_t mpfr_get_emin_min (void);
[all …]
Dsum.c151 unsigned long n, mpfr_exp_t minexp, mpfr_exp_t maxexp, in sum_raw()
153 mpfr_exp_t *ep, mpfr_exp_t *minexpp, mpfr_exp_t *maxexpp) in sum_raw()
174 mpfr_exp_t maxexp2 = MPFR_EXP_MIN; in sum_raw()
189 mpfr_exp_t xe, vd; in sum_raw()
414 mpfr_exp_t e; /* exponent of the computed result */ in sum_raw()
415 mpfr_exp_t err; /* exponent of the error bound */ in sum_raw()
466 mpfr_exp_t diffexp; in sum_raw()
528 mpfr_exp_t maxexp, unsigned long rn) in sum_aux()
615 mpfr_exp_t minexp; /* exponent of the LSB of the block for sum_raw */ in sum_aux()
617 mpfr_exp_t e; /* temporary exponent of the result */ in sum_aux()
[all …]
Dli2.c39 mpfr_exp_t se, err; in li2_series()
83 if (MPFR_GET_EXP (w) <= se - (mpfr_exp_t) p) in li2_series()
91 if (MPFR_CAN_ROUND (s, (mpfr_exp_t) p - err, sump, rnd_mode)) in li2_series()
151 if ((MPFR_EXP (x) >= (mpfr_exp_t) w - MPFR_EXP (g)) && in mpfr_li2_asympt_pos()
198 if ((MPFR_EXP (x) >= (mpfr_exp_t) (w - 2) - MPFR_EXP (g)) && in mpfr_li2_asympt_neg()
218 mpfr_exp_t err; in mpfr_li2()
268 mpfr_exp_t expo_l; in mpfr_li2()
294 if (MPFR_CAN_ROUND (s, (mpfr_exp_t) m - err, yp, rnd_mode)) in mpfr_li2()
341 mpfr_exp_t expo_l; in mpfr_li2()
390 if (MPFR_CAN_ROUND (s, (mpfr_exp_t) m - err, yp, rnd_mode)) in mpfr_li2()
[all …]
Drem1.c59 mpfr_exp_t ex, ey; in mpfr_rem1()
107 mpfr_exp_t k = mpz_scan1 (my, 0); in mpfr_rem1()
121 if (ex + (mpfr_exp_t) mpz_sizeinbase (mx, 2) < in mpfr_rem1()
122 ey + (mpfr_exp_t) mpz_sizeinbase (my, 2)) in mpfr_rem1()
214 if (ex + (mpfr_exp_t) mpz_sizeinbase (r, 2) < in mpfr_rem1()
215 ey + (mpfr_exp_t) mpz_sizeinbase (my, 2)) in mpfr_rem1()
Dlog2p1.c31 static mpfr_exp_t
40 mpfr_exp_t e; in mpfr_log2p1_isexact()
56 mpfr_exp_t expx = MPFR_GET_EXP(x); in mpfr_log2p1_special()
57 mpfr_exp_t k = expx - 1, expk; in mpfr_log2p1_special()
139 mpfr_exp_t k; in mpfr_log2p1()
Droot.c51 mpfr_exp_t e, r, sh, f; in mpfr_rootn_ui()
134 r = e % (mpfr_exp_t) k; in mpfr_rootn_ui()
147 if ((mpfr_exp_t) size_m + r >= k * (mpfr_exp_t) n) in mpfr_rootn_ui()
150 f = (k * (mpfr_exp_t) n - (mpfr_exp_t) size_m - r) / k; in mpfr_rootn_ui()
185 MPFR_SET_EXP (y, MPFR_GET_EXP (y) + e / (mpfr_exp_t) k); in mpfr_rootn_ui()
235 mpfr_exp_t expt; in mpfr_root_aux()
Dfpif.c305 mpfr_exp_t exponent = MPFR_GET_EXP (x); in mpfr_fpif_store_exponent()
345 sizeof(mpfr_exp_t), exponent_size); in mpfr_fpif_store_exponent()
378 mpfr_exp_t exponent; in mpfr_fpif_read_exponent_from_file()
382 unsigned char buffer[sizeof(mpfr_exp_t)]; in mpfr_fpif_read_exponent_from_file()
407 exponent_size > sizeof(mpfr_exp_t))) in mpfr_fpif_read_exponent_from_file()
415 sizeof(mpfr_exp_t), exponent_size); in mpfr_fpif_read_exponent_from_file()
425 exponent = exp_sign_bit ? - (mpfr_exp_t) uexp : (mpfr_exp_t) uexp; in mpfr_fpif_read_exponent_from_file()
536 buf_size = MAX_VARIABLE_STORAGE(sizeof(mpfr_exp_t), mpfr_get_prec (x)); in mpfr_fpif_export()
Dmpn_exp.c30 nbits_mpfr_exp_t (mpfr_exp_t e) in nbits_mpfr_exp_t()
76 mpfr_mpn_exp (mp_limb_t *a, mpfr_exp_t *exp_r, int b, mpfr_exp_t e, size_t n) in mpfr_mpn_exp()
79 mpfr_exp_t f, h; in mpfr_mpn_exp()
145 mpfr_exp_t, mpfr_uexp_t, in mpfr_mpn_exp()
Dubf.c45 mpfr_exp_t ex = MPFR_GET_EXP (x); in mpfr_init_get_zexp()
56 MPFR_TMP_INIT1 (e_limb, e, sizeof (mpfr_exp_t) * CHAR_BIT); in mpfr_init_get_zexp()
113 mpfr_exp_t e; in mpfr_ubf_mul_exact()
216 mpfr_exp_t
245 mpfr_exp_t
249 mpfr_exp_t e; in mpfr_ubf_diff_exp()
Dyn.c45 static mpfr_exp_t
50 mpfr_exp_t e, emax; in mpfr_yn_s1()
84 static mpfr_exp_t
90 mpfr_exp_t exps, expU; in mpfr_yn_s3()
135 if (MPFR_EXP (u) + (mpfr_exp_t) MPFR_PREC (u) < MPFR_EXP (s) && in mpfr_yn_s3()
217 if (n == 0 && MPFR_EXP(z) < - (mpfr_exp_t) (MPFR_PREC(res) / 2)) in mpfr_yn()
271 if (n == 1 && MPFR_EXP(z) + 1 < - (mpfr_exp_t) MPFR_PREC(res)) in mpfr_yn()
275 mpfr_exp_t err1; in mpfr_yn()
302 err1 = (mpfr_exp_t) MPFR_PREC(y) - MPFR_EXP(y) + 1; in mpfr_yn()
323 mpfr_exp_t err1, err2, err3; in mpfr_yn()
Dpow.c27 # define MPFR_POW_EXP_THRESHOLD (MAX (sizeof(mpfr_exp_t) * CHAR_BIT, 256))
41 mpfr_exp_t d, b, i; in mpfr_pow_is_exact()
173 mpfr_exp_t err, exp_t; in mpfr_pow_general()
254 Ntmin = sizeof(mpfr_exp_t) * CHAR_BIT; in mpfr_pow_general()
545 mpfr_init2 (t, sizeof (mpfr_exp_t) * CHAR_BIT); in mpfr_pow()
578 MPFR_TMP_INIT1 (tmp_limb, tmp, sizeof (mpfr_exp_t) * CHAR_BIT); in mpfr_pow()
631 mpfr_exp_t b = MPFR_GET_EXP (x) - 1; in mpfr_pow()
635 sizeof(mpfr_exp_t) * CHAR_BIT); in mpfr_pow()
646 MPFR_ASSERTD (MPFR_EXP (y) > sizeof(mpfr_exp_t) * CHAR_BIT); in mpfr_pow()
679 mpfr_exp_t err, expx, logt; in mpfr_pow()
Derfc.c33 static mpfr_exp_t
39 mpfr_exp_t exp_err; in mpfr_erfc_asympt()
63 if (MPFR_GET_EXP (t) + (mpfr_exp_t) prec <= MPFR_GET_EXP (y)) in mpfr_erfc_asympt()
128 mpfr_exp_t te, err; in mpfr_erfc()
130 mpfr_exp_t emin = mpfr_get_emin (); in mpfr_erfc()
173 mpfr_exp_t e = MPFR_EXP(x); in mpfr_erfc()
Ddigamma.c35 static mpfr_exp_t
40 mpfr_exp_t e, exps, f, expu; in mpfr_digamma_approx()
83 if (expu < exps - (mpfr_exp_t) p) in mpfr_digamma_approx()
120 mpfr_exp_t e1, expv, expx, q; in mpfr_digamma_reflection()
162 e1 = MPFR_GET_EXP(t) - (mpfr_exp_t) p + 1; /* bound for t: err(t) <= 2^e1 */ in mpfr_digamma_reflection()
170 e1 += (mpfr_exp_t) p - MPFR_EXP(t); /* error is now 2^e1 ulps */ in mpfr_digamma_reflection()
212 mpfr_exp_t errt, erru, expt; in mpfr_digamma_positive()
397 if (MPFR_EXP(x) <= -2 * (mpfr_exp_t) MAX(MPFR_PREC(x), MPFR_PREC(y))) in mpfr_digamma()
Dmul_2ui.c40 mpfr_exp_t exp = MPFR_GET_EXP (y); in mpfr_mul_2ui()
41 mpfr_exp_t diffexp; in mpfr_mul_2ui()
48 MPFR_SET_EXP (y, exp + (mpfr_exp_t) n); in mpfr_mul_2ui()
/netbsd/src/external/lgpl3/mpc/dist/tests/
Drandom.c113 test_default_random (mpc_ptr z, mpfr_exp_t emin, mpfr_exp_t emax, in test_default_random()
153 mpfr_set_exp (mpc_realref (z), (mpfr_exp_t) gmp_urandomm_ui (rands, range) + emin); in test_default_random()
156 mpfr_set_exp (mpc_imagref (z), (mpfr_exp_t) gmp_urandomm_ui (rands, range) + emin); in test_default_random()
202 test_random_mpfr (mpfr_ptr x, mpfr_exp_t emin, mpfr_exp_t emax, in test_random_mpfr()
220 mpfr_set_exp (x, (mpfr_exp_t) gmp_urandomm_ui (rands, range) + emin); in test_random_mpfr()
245 test_random_mpc (mpc_ptr z, mpfr_exp_t emin, mpfr_exp_t emax, in test_random_mpc()
264 (mpfr_exp_t) gmp_urandomm_ui (rands, range) + emin); in test_random_mpc()
267 (mpfr_exp_t) gmp_urandomm_ui (rands, range) + emin); in test_random_mpc()
Dtgeneric.tpl28 mpfr_exp_t exp_min, mpfr_exp_t exp_max,
32 mpfr_exp_t exp_min,
33 mpfr_exp_t exp_max,
49 mpfr_exp_t exp_max)
54 mpfr_exp_t exp_min;
92 mpfr_exp_t exp_min, mpfr_exp_t exp_max,
213 special_random_mpc (mpc_ptr z, mpfr_exp_t exp_min, mpfr_exp_t exp_max,
239 mpfr_exp_t exp_min, mpfr_exp_t exp_max,
Dtacos.c35 mpfr_exp_t emin = mpfr_get_emin (); in test_reduced()
36 mpfr_exp_t emax = mpfr_get_emax (); in test_reduced()
67 mpfr_exp_t emin = mpfr_get_emin (); in test_reduced2()
68 mpfr_exp_t emax = mpfr_get_emax (); in test_reduced2()
/netbsd/src/external/lgpl3/mpfr/dist/tests/
Dmpfr-test.h209 void mpfr_random2 (mpfr_ptr, mp_size_t, mpfr_exp_t, gmp_randstate_t);
219 void set_emin (mpfr_exp_t);
220 void set_emax (mpfr_exp_t);
221 void tests_default_random (mpfr_ptr, int, mpfr_exp_t, mpfr_exp_t,
225 const char *, int, mpfr_exp_t, mpfr_exp_t,
288 mpfr_exp_t ex; in mpfr_print_raw()
/netbsd/src/external/lgpl3/mpc/dist/src/
Dsqr.c64 mpfr_exp_t ea, ec; in mpfr_fsss()
92 mpfr_exp_t emax = mpfr_get_emax (); in mpfr_fsss()
104 mpfr_set_exp (v, (mpfr_exp_t) mpz_get_ui (ev)); in mpfr_fsss()
111 mpfr_set_exp (u, (mpfr_exp_t) mpz_get_ui (eu)); in mpfr_fsss()
126 mpfr_exp_t emin = mpfr_get_emin (); in mpfr_fsss()
131 mpfr_set_exp (v, (mpfr_exp_t) mpz_get_si (ev)); in mpfr_fsss()
137 mpfr_set_exp (u, (mpfr_exp_t) mpz_get_si (eu)); in mpfr_fsss()
172 mpfr_exp_t emin; in mpc_sqr()
238 if (SAFE_ABS (mpfr_exp_t, in mpc_sqr()
240 > (mpfr_exp_t) MPC_MAX_PREC (op) / 2) { in mpc_sqr()
Dpow.c105 mpfr_exp_t ey; in fix_sign()
114 ey += (mpfr_exp_t) t; in fix_sign()
176 mpfr_exp_t ec, ed, ey; in mpc_pow_exact()
204 ey += (mpfr_exp_t) t; in mpc_pow_exact()
248 ec += (mpfr_exp_t) t; in mpc_pow_exact()
254 ec += (mpfr_exp_t) t; in mpc_pow_exact()
265 ec += (mpfr_exp_t) t; in mpc_pow_exact()
321 ec = -ec - (mpfr_exp_t) t; in mpc_pow_exact()
357 ed += (mpfr_exp_t) w; in mpc_pow_exact()
363 ed += (mpfr_exp_t) w; in mpc_pow_exact()
[all …]

12345678910>>...12