Home
last modified time | relevance | path

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

123

/dragonfly/contrib/mpfr/src/
HDdigamma.c40 MPFR_ASSERTN(MPFR_IS_POS(x) && (MPFR_EXP(x) >= 2)); in mpfr_digamma_approx()
133 if (MPFR_EXP(x) < 0) in mpfr_digamma_reflection()
134 q = MPFR_PREC(x) + 1 - MPFR_EXP(x); in mpfr_digamma_reflection()
135 else if (MPFR_EXP(x) <= MPFR_PREC(x)) in mpfr_digamma_reflection()
138 q = MPFR_EXP(x); in mpfr_digamma_reflection()
160 e1 = MPFR_EXP(t) - (mpfr_exp_t) p + 1; /* bound for t: err(t) <= 2^e1 */ in mpfr_digamma_reflection()
163 if (MPFR_EXP(t) > 0) in mpfr_digamma_reflection()
164 e1 = e1 + 2 * MPFR_EXP(t) + 1; in mpfr_digamma_reflection()
168 e1 += (mpfr_exp_t) p - MPFR_EXP(t); /* error is now 2^e1 ulps */ in mpfr_digamma_reflection()
172 expv = MPFR_EXP(v); in mpfr_digamma_reflection()
[all …]
HDjyn_asympt.c128 if (MPFR_EXP(err_u) <= MPFR_EXP(Q)) in FUNCTION()
141 if (MPFR_EXP(err_u) <= MPFR_EXP(P)) in FUNCTION()
177 err = MPFR_EXP(c); in FUNCTION()
178 if (MPFR_EXP(s) > err) in FUNCTION()
179 err = MPFR_EXP(s); in FUNCTION()
196 err = MPFR_EXP(c); in FUNCTION()
197 if (MPFR_EXP(s) > err) in FUNCTION()
198 err = MPFR_EXP(s); in FUNCTION()
207 if (MPFR_EXP(s) > err) in FUNCTION()
208 err = MPFR_EXP(s); in FUNCTION()
[all …]
HDyn.c56 emax = MPFR_EXP(s); in mpfr_yn_s1()
65 e = MPFR_EXP(s); in mpfr_yn_s1()
73 return emax - MPFR_EXP(s); in mpfr_yn_s1()
111 exps = MPFR_EXP (s); in mpfr_yn_s3()
128 exps = MPFR_EXP (u); in mpfr_yn_s3()
132 exps = MPFR_EXP (s); in mpfr_yn_s3()
135 if (MPFR_EXP (u) + (mpfr_exp_t) MPFR_PREC (u) < MPFR_EXP (s) && in mpfr_yn_s3()
143 exps = expU - 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()
[all …]
HDswap.c45 e1 = MPFR_EXP(u); in mpfr_swap()
46 e2 = MPFR_EXP(v); in mpfr_swap()
47 MPFR_EXP(v) = e1; in mpfr_swap()
48 MPFR_EXP(u) = e2; in mpfr_swap()
HDpow_z.c61 rnd1 = (MPFR_EXP(x) >= 1) ? MPFR_RNDZ in mpfr_pow_pos_z()
63 rnd2 = (MPFR_EXP(x) >= 1) ? MPFR_RNDD : MPFR_RNDU; in mpfr_pow_pos_z()
129 MPFR_EXP (x) - 1) != 0); in mpfr_pow_pos_z()
230 MPFR_EXP (x) - 1) == 0)) in mpfr_pow_z()
232 mpfr_exp_t expx = MPFR_EXP (x); /* warning: x and y may be the same in mpfr_pow_z()
288 rnd1 = MPFR_EXP (x) < 1 ? MPFR_RNDZ : in mpfr_pow_z()
337 MPFR_EXP (x) - 1) != 0); in mpfr_pow_z()
HDfrac.c121 MPFR_EXP (t) = 0; /* should be re, but not necessarily in the range */ in mpfr_frac()
124 MPFR_EXP (r) += re; in mpfr_frac()
130 MPFR_EXP (r) = re; in mpfr_frac()
138 MPFR_EXP (r) ++; in mpfr_frac()
HDset.c35 MPFR_EXP (a) = MPFR_EXP (b); in mpfr_set4()
60 if (MPFR_UNLIKELY ( ++MPFR_EXP (a) > __gmpfr_emax)) in mpfr_set4()
HDconst_euler.c63 exp_S = MPFR_EXP(y); in mpfr_const_euler_internal()
70 err = 1 + MAX(exp_S + 2, MPFR_EXP(z)) - MPFR_EXP(y); in mpfr_const_euler_internal()
72 exp_S = MPFR_EXP(y); in mpfr_const_euler_internal()
79 err = err + exp_S - MPFR_EXP(y); in mpfr_const_euler_internal()
HDpow_si.c70 if (mpfr_cmp_si_2exp (x, MPFR_SIGN(x), MPFR_EXP(x) - 1) == 0) in mpfr_pow_si()
72 mpfr_exp_t expx = MPFR_EXP (x) - 1, expy; in mpfr_pow_si()
167 rnd1 = MPFR_EXP (x) < 1 ? MPFR_RNDZ : in mpfr_pow_si()
212 MPFR_EXP (x) - 1) != 0); in mpfr_pow_si()
HDmpfr-impl.h663 #define MPFR_EXP(x) ((x)->_mpfr_exp) macro
737 # define MPFR_SET_INVALID_EXP(x) ((void) (MPFR_EXP (x) = MPFR_EXP_INVALID))
739 # define MPFR_GET_EXP(x) MPFR_EXP (x)
740 # define MPFR_SET_EXP(x, exp) ((void) (MPFR_EXP (x) = (exp)))
755 #define MPFR_IS_NAN(x) (MPFR_EXP(x) == MPFR_EXP_NAN)
756 #define MPFR_SET_NAN(x) (MPFR_EXP(x) = MPFR_EXP_NAN)
757 #define MPFR_IS_INF(x) (MPFR_EXP(x) == MPFR_EXP_INF)
758 #define MPFR_SET_INF(x) (MPFR_EXP(x) = MPFR_EXP_INF)
759 #define MPFR_IS_ZERO(x) (MPFR_EXP(x) == MPFR_EXP_ZERO)
760 #define MPFR_SET_ZERO(x) (MPFR_EXP(x) = MPFR_EXP_ZERO)
[all …]
HDcheck.c70 if ((MPFR_EXP (x) < __gmpfr_emin) || (MPFR_EXP (x) > __gmpfr_emax)) in mpfr_check()
HDcmp_abs.c62 be = MPFR_EXP (b); in mpfr_cmpabs()
63 ce = MPFR_EXP (c); in mpfr_cmpabs()
HDli2.c155 MPFR_ASSERTN (MPFR_EXP (g) > MPFR_EXP (h)); in mpfr_li2_asympt_pos()
162 if ((MPFR_EXP (x) >= (mpfr_exp_t) w - MPFR_EXP (g)) && in mpfr_li2_asympt_pos()
199 MPFR_ASSERTN (MPFR_EXP (g) >= MPFR_EXP (h)); in mpfr_li2_asympt_neg()
205 if ((MPFR_EXP (x) >= (mpfr_exp_t) (w - 2) - MPFR_EXP (g)) && in mpfr_li2_asympt_neg()
HDrec_sqrt.c533 if (s == 0 && mpfr_cmp_ui_2exp (u, 1, MPFR_EXP(u) - 1) == 0) in mpfr_rec_sqrt()
548 MPFR_EXP(r) = - (MPFR_EXP(u) - 1 - s) / 2; in mpfr_rec_sqrt()
551 MPFR_EXP(r) ++; in mpfr_rec_sqrt()
HDstack_interface.c50 return MPFR_EXP (x); in mpfr_custom_get_exp()
86 MPFR_EXP (x) = e; in mpfr_custom_init_set()
HDlog1p.c117 if (MPFR_EXP(x) < 0) in mpfr_log1p()
118 Nt += -MPFR_EXP(x); in mpfr_log1p()
HDnext.c51 mpfr_exp_t exp = MPFR_EXP (x); in mpfr_nexttozero()
85 mpfr_exp_t exp = MPFR_EXP (x); in mpfr_nexttoinf()
HDagm.c261 MPFR_EXP (u) += scale2; in mpfr_agm()
262 MPFR_EXP (v) += scale2; in mpfr_agm()
303 MPFR_EXP (r) -= scaleop + scaleit; in mpfr_agm()
HDsin_cos.c155 || MPFR_EXP(xr) < (mpfr_exp_t) 3 - (mpfr_exp_t) m in mpfr_sin_cos()
156 || MPFR_EXP(c) < (mpfr_exp_t) 3 - (mpfr_exp_t) m) in mpfr_sin_cos()
611 mpfr_init2 (pi, (MPFR_EXP(x) > 0) ? w + MPFR_EXP(x) : w); in mpfr_sincos_fast()
645 errs = err - MPFR_EXP(ts); in mpfr_sincos_fast()
646 errc = err - MPFR_EXP(tc); in mpfr_sincos_fast()
HDlngamma.c188 if (MPFR_EXP(z0) <= - (mpfr_exp_t) MPFR_PREC(y)) in GAMMA_FUNC()
253 while (prec <= -MPFR_EXP(z0)); in GAMMA_FUNC()
687 if (MPFR_EXP(x) + 1 <= - (mpfr_exp_t) MPFR_PREC(y)) in mpfr_lgamma()
720 expl = MPFR_EXP (l); in mpfr_lgamma()
727 if (expl < MPFR_EXP(x) + (mpfr_exp_t) w) in mpfr_lgamma()
HDset_d.c214 MPFR_EXP(tmp) = __gmpfr_extract_double (tmpmant, d); in mpfr_set_d()
245 MPFR_EXP(tmp) -= (mpfr_exp_t) (cnt + k * GMP_NUMB_BITS); in mpfr_set_d()
HDget_exp.c30 return MPFR_EXP(x); /* do not use MPFR_GET_EXP of course... */ in mpfr_get_exp()
HDset_exp.c30 MPFR_EXP(x) = exponent; /* do not use MPFR_SET_EXP of course... */ in mpfr_set_exp()
HDinit2.c58 MPFR_EXP (x) = MPFR_EXP_INVALID; /* make sure that the exp field has a in mpfr_init2()
HDcsc.c49 if (MPFR_EXP(x) <= -2 * (mpfr_exp_t) MAX(MPFR_PREC(x), MPFR_PREC(y))) \

123