Home
last modified time | relevance | path

Searched refs:MPFR_GET_EXP (Results 1 – 25 of 95) sorted by relevance

1234

/dragonfly/contrib/mpfr/src/
H A Dli2.c65 se = MPFR_GET_EXP (s); in li2_series()
87 - MPFR_GET_EXP (s); in li2_series()
89 se = MPFR_GET_EXP (s); in li2_series()
375 3 + MAX (1, -expo_l) + MPFR_GET_EXP (u)) - MPFR_GET_EXP (s); in mpfr_li2()
383 err = MAX (err, 3 + MPFR_GET_EXP (u)) - MPFR_GET_EXP (s); in mpfr_li2()
425 e1 = MPFR_GET_EXP (s); in mpfr_li2()
433 e2 = MPFR_GET_EXP (u); in mpfr_li2()
491 err = MAX (err, 2 + MPFR_GET_EXP (u)) - MPFR_GET_EXP (s); in mpfr_li2()
593 + MPFR_GET_EXP (s); in mpfr_li2()
598 err = MAX (err, 3 + MPFR_GET_EXP(w)) - MPFR_GET_EXP (s); in mpfr_li2()
[all …]
H A Deint.c51 if (MPFR_GET_EXP(x) <= - (mpfr_exp_t) w) in mpfr_eint_aux()
110 if (sizeinbase < MPFR_GET_EXP (errs)) in mpfr_eint_aux()
136 e = MPFR_GET_EXP (errs) - MPFR_GET_EXP (y); in mpfr_eint_aux()
166 for (k = 1; MPFR_GET_EXP(t) + (mpfr_exp_t) p > MPFR_GET_EXP(y); k++) in mpfr_eint_asympt()
174 mpfr_mul_2si (err, err, MPFR_GET_EXP(y) - MPFR_GET_EXP(t), MPFR_RNDU); in mpfr_eint_asympt()
177 mpfr_div_2si (err, err, MPFR_GET_EXP(y) - MPFR_GET_EXP(t), MPFR_RNDU); in mpfr_eint_asympt()
186 err_exp = MPFR_GET_EXP(err); in mpfr_eint_asympt()
265 if (MPFR_GET_EXP(x) == -1) /* 1/4 <= x < 1/2 */ in mpfr_eint()
288 te = MPFR_GET_EXP(tmp); in mpfr_eint()
296 te = MPFR_GET_EXP(tmp); in mpfr_eint()
[all …]
H A Derfc.c60 mpfr_mul_2si (err, err, MPFR_GET_EXP (y) - MPFR_GET_EXP (t), MPFR_RNDU); in mpfr_erfc_asympt()
62 mpfr_div_2si (err, err, MPFR_GET_EXP (y) - MPFR_GET_EXP (t), MPFR_RNDU); in mpfr_erfc_asympt()
63 if (MPFR_GET_EXP (t) + (mpfr_exp_t) prec <= MPFR_GET_EXP (y)) in mpfr_erfc_asympt()
114 exp_err = MPFR_GET_EXP (err); in mpfr_erfc_asympt()
222 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, __gmpfr_one, - MPFR_GET_EXP (x) - 1, in mpfr_erfc()
227 if (MPFR_GET_EXP (x) > 0) in mpfr_erfc()
228 prec += 2 * MPFR_GET_EXP(x); in mpfr_erfc()
238 2 * MPFR_GET_EXP (x) - 2 >= MPFR_INT_CEIL_LOG2 (prec)) in mpfr_erfc()
253 te = MPFR_GET_EXP (tmp); in mpfr_erfc()
262 err = MAX (te - MPFR_GET_EXP (tmp), 0) + 1; in mpfr_erfc()
H A Dai.c155 if ( (MPFR_IS_ZERO (x)) || (MPFR_GET_EXP (x) <= 0) ) in mpfr_ai1()
169 if (MPFR_GET_EXP (x) <= 0) in mpfr_ai1()
238 || (MPFR_GET_EXP (ti) + (mpfr_exp_t)prec + 3 <= MPFR_GET_EXP (s)); in mpfr_ai1()
240 || (MPFR_GET_EXP (tip1) + (mpfr_exp_t)prec + 3 <= MPFR_GET_EXP (s)); in mpfr_ai1()
290 - MPFR_GET_EXP (s); in mpfr_ai1()
378 if ( (MPFR_IS_ZERO (x)) || (MPFR_GET_EXP (x) <= 0) ) in mpfr_ai2()
392 if (MPFR_GET_EXP (x) <= 0) in mpfr_ai2()
526 MPFR_GET_EXP (u0) + (mpfr_exp_t)prec + 4 <= MPFR_GET_EXP (result); in mpfr_ai2()
528 MPFR_GET_EXP (u1) + (mpfr_exp_t)prec + 4 <= MPFR_GET_EXP (result); in mpfr_ai2()
537 + cond - MPFR_GET_EXP (result)); in mpfr_ai2()
[all …]
H A Dgamma.c42 return (MPFR_GET_EXP(x) <= 1) ? MPFR_PREC(x) + 2 - MPFR_GET_EXP(x) in mpfr_gamma_2_minus_x_exact()
44 : MPFR_GET_EXP(x) - 1); in mpfr_gamma_2_minus_x_exact()
52 return MPFR_PREC(x) - MPFR_GET_EXP(x); in mpfr_gamma_1_minus_x_exact()
53 else if (MPFR_GET_EXP(x) <= 0) in mpfr_gamma_1_minus_x_exact()
54 return MPFR_PREC(x) + 1 - MPFR_GET_EXP(x); in mpfr_gamma_1_minus_x_exact()
55 else if (MPFR_PREC(x) >= MPFR_GET_EXP(x)) in mpfr_gamma_1_minus_x_exact()
58 return MPFR_GET_EXP(x); in mpfr_gamma_1_minus_x_exact()
167 if (MPFR_GET_EXP (x) + 2 in mpfr_gamma()
179 MPFR_GET_EXP (x) == 1 - MPFR_EMAX_MAX && in mpfr_gamma()
271 expxp = MPFR_GET_EXP (xp); in mpfr_gamma()
[all …]
H A Dsubnormal.c53 || (MPFR_GET_EXP (y) >= in mpfr_subnormalize()
61 if (MPFR_GET_EXP (y) == __gmpfr_emin) in mpfr_subnormalize()
121 MPFR_ASSERTD (MPFR_GET_EXP (y) > __gmpfr_emin); in mpfr_subnormalize()
124 q = (mpfr_uexp_t) MPFR_GET_EXP (y) - __gmpfr_emin + 1; in mpfr_subnormalize()
130 MPFR_SET_EXP (dest, MPFR_GET_EXP (y)); in mpfr_subnormalize()
134 MPFR_SET_EXP (dest, MPFR_GET_EXP (dest) + 1)); in mpfr_subnormalize()
H A Dfms.c135 MPFR_GET_EXP (x) + MPFR_GET_EXP (y) >= __gmpfr_emax + 3) in mpfr_fms()
158 if (MPFR_GET_EXP (u) > MPFR_GET_EXP (z) && in mpfr_fms()
159 MPFR_GET_EXP (u) - MPFR_GET_EXP (z) > MPFR_PREC (u)) in mpfr_fms()
217 diffexp = MPFR_GET_EXP (z) - __gmpfr_emin; in mpfr_fms()
235 if (MPFR_GET_EXP (x) < MPFR_GET_EXP (y)) in mpfr_fms()
H A Dfma.c133 MPFR_GET_EXP (x) + MPFR_GET_EXP (y) >= __gmpfr_emax + 3) in mpfr_fma()
156 if (MPFR_GET_EXP (u) > MPFR_GET_EXP (z) && in mpfr_fma()
157 MPFR_GET_EXP (u) - MPFR_GET_EXP (z) > MPFR_PREC (u)) in mpfr_fma()
215 diffexp = MPFR_GET_EXP (z) - __gmpfr_emin; in mpfr_fma()
233 if (MPFR_GET_EXP (x) < MPFR_GET_EXP (y)) in mpfr_fma()
H A Dlngamma.c123 expo = MPFR_GET_EXP (x); in unit_bit()
297 err_u = MPFR_GET_EXP(s) + __gmpfr_ceil_log2 (d) + 1 - MPFR_GET_EXP(u); in GAMMA_FUNC()
302 err_s = MPFR_GET_EXP(s); /* 2-x <= 2^err_s */ in GAMMA_FUNC()
307 err_s += 3 - MPFR_GET_EXP(s); in GAMMA_FUNC()
344 err_s += 1 - MPFR_GET_EXP(v); in GAMMA_FUNC()
446 for (m = 2; MPFR_GET_EXP(v) + (mpfr_exp_t) w >= MPFR_GET_EXP(s); m++) in GAMMA_FUNC()
472 MPFR_ASSERTD(MPFR_GET_EXP(v) <= - (2 * m + 3)); in GAMMA_FUNC()
512 err_s = MPFR_GET_EXP(s); in GAMMA_FUNC()
550 err_t = MPFR_GET_EXP(t) + (mpfr_exp_t) in GAMMA_FUNC()
552 err_s = MPFR_GET_EXP(s) + (mpfr_exp_t) in GAMMA_FUNC()
[all …]
H A Dtanh.c62 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, xt, -2 * MPFR_GET_EXP (xt), 1, 0, in mpfr_tanh()
98 if (MPFR_GET_EXP (x) < 0) in mpfr_tanh()
99 Nt += -MPFR_GET_EXP (x); in mpfr_tanh()
121 d = MPFR_GET_EXP (te); /* For Error calculation */ in mpfr_tanh()
124 d = d - MPFR_GET_EXP (te); in mpfr_tanh()
138 if (MPFR_GET_EXP (t) == 1) in mpfr_tanh()
H A Dpow.c45 MPFR_GET_EXP (x) - 1) != 0); in mpfr_pow_is_exact()
122 expo = MPFR_GET_EXP (y); in is_odd()
227 err = MPFR_NOTZERO (t) && MPFR_GET_EXP (t) >= -1 ? in mpfr_pow_general()
228 MPFR_GET_EXP (t) + 3 : 1; in mpfr_pow_general()
231 if (MPFR_GET_EXP (k) > err) in mpfr_pow_general()
232 err = MPFR_GET_EXP (k); in mpfr_pow_general()
591 if (MPFR_IS_NEG (y) ? (MPFR_GET_EXP (x) > 1) : (MPFR_GET_EXP (x) < 0)) in mpfr_pow()
600 inex2 = mpfr_set_exp_t (tmp, MPFR_GET_EXP (x), MPFR_RNDN); in mpfr_pow()
634 if (y_is_integer && (MPFR_GET_EXP (y) <= 256)) in mpfr_pow()
649 mpfr_exp_t b = MPFR_GET_EXP (x) - 1; in mpfr_pow()
[all …]
H A Dsin.c70 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, x, -2 * MPFR_GET_EXP (x), 2, 0, in mpfr_sin()
82 expx = MPFR_GET_EXP (x); in mpfr_sin()
118 || MPFR_GET_EXP(xr) < (mpfr_exp_t) 3 - (mpfr_exp_t) m in mpfr_sin()
120 || MPFR_GET_EXP(c) < (mpfr_exp_t) 3 - (mpfr_exp_t) m) in mpfr_sin()
156 err = 2 * MPFR_GET_EXP (c) + (mpfr_exp_t) m - 3 - (reduce != 0); in mpfr_sin()
164 if (MPFR_GET_EXP (c) == 1) in mpfr_sin()
H A Dcmp2.c54 if (MPFR_GET_EXP (b) >= MPFR_GET_EXP (c)) in mpfr_cmp2()
57 diff_exp = (mpfr_uexp_t) MPFR_GET_EXP (b) - MPFR_GET_EXP (c); in mpfr_cmp2()
121 diff_exp = (mpfr_uexp_t) MPFR_GET_EXP (c) - MPFR_GET_EXP (b); in mpfr_cmp2()
H A Dacosh.c106 pln2 = Nt - MPFR_PREC_MIN < MPFR_GET_EXP (t) ? in mpfr_acosh()
107 MPFR_PREC_MIN : Nt - MPFR_GET_EXP (t); in mpfr_acosh()
116 exp_te = MPFR_GET_EXP (t); in mpfr_acosh()
130 d = exp_te - MPFR_GET_EXP (t); in mpfr_acosh()
136 err = 3 + MAX (1, d) - MPFR_GET_EXP (t); in mpfr_acosh()
H A Djn.c134 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (res, __gmpfr_one, -2 * MPFR_GET_EXP (z), in mpfr_jn()
163 (res, z, -2 * MPFR_GET_EXP (z), 3, 0, r, { in mpfr_jn()
214 if ((MPFR_GET_EXP (y) < 0 && absn > - expo.saved_emin) in mpfr_jn()
216 MPFR_GET_EXP (y) < expo.saved_emin / (mpfr_exp_t) absn)) in mpfr_jn()
261 exps = MPFR_IS_ZERO (s) ? MPFR_EMIN_MIN : MPFR_GET_EXP (s); in mpfr_jn()
280 exps = MPFR_IS_ZERO (s) ? MPFR_EMIN_MIN : MPFR_GET_EXP (t); in mpfr_jn()
284 exps = MPFR_IS_ZERO (s) ? MPFR_EMIN_MIN : MPFR_GET_EXP (s); in mpfr_jn()
288 if (MPFR_GET_EXP (t) + (mpfr_exp_t) prec <= exps && in mpfr_jn()
H A Dsum.c77 if (MPFR_GET_EXP (tab[i]) < min) in mpfr_sum_sort()
78 min = MPFR_GET_EXP(tab[i]); in mpfr_sum_sort()
79 if (MPFR_GET_EXP (tab[i]) > max) in mpfr_sum_sort()
80 max = MPFR_GET_EXP(tab[i]); in mpfr_sum_sort()
95 #define GET_EXP1(x) (MPFR_IS_ZERO (x) ? min : MPFR_GET_EXP (x))
128 #define GET_EXP2(x) (MPFR_IS_ZERO (x) ? MPFR_EMIN_MIN : MPFR_GET_EXP (x))
299 MPFR_GET_EXP (cur_sum) - prec + 2, in mpfr_sum()
H A Dsinh.c63 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, xt, -2 * MPFR_GET_EXP(xt), 2, 1, in mpfr_sinh()
84 if (MPFR_GET_EXP (x) < 0) in mpfr_sinh()
85 Nt -= 2*MPFR_GET_EXP (x); in mpfr_sinh()
149 d = MPFR_GET_EXP (t); in mpfr_sinh()
161 d = d - MPFR_GET_EXP (t) + 2; in mpfr_sinh()
H A Dadd.c101 if (MPFR_GET_EXP(b) < MPFR_GET_EXP(c)) in mpfr_add()
106 if (MPFR_GET_EXP(b) < MPFR_GET_EXP(c)) in mpfr_add()
H A Datanh.c64 if (MPFR_UNLIKELY (MPFR_GET_EXP (xt) > 0)) in mpfr_atanh()
66 if (MPFR_GET_EXP (xt) == 1 && mpfr_powerof2_raw (xt)) in mpfr_atanh()
78 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, xt, -2 * MPFR_GET_EXP (xt), 1, 1, in mpfr_atanh()
109 err = Nt - (MAX (4 - MPFR_GET_EXP (t), 0) + 1); in mpfr_atanh()
H A Dasin.c58 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (asin, x, -2 * MPFR_GET_EXP (x), 2, 1, in mpfr_asin()
95 MPFR_ASSERTD (MPFR_GET_EXP (xp) <= 0); in mpfr_asin()
96 MPFR_ASSERTD (MPFR_GET_EXP (x) <= 0); in mpfr_asin()
97 xp_exp = 2 - MPFR_GET_EXP (xp); in mpfr_asin()
H A Dcos.c163 expx = MPFR_GET_EXP (x); in mpfr_cos()
219 K = K0 + 1 + MAX(0, MPFR_GET_EXP(r)) / 2; in mpfr_cos()
224 MPFR_SET_EXP (r, MPFR_GET_EXP (r) - 2 * K); /* Can't overflow! */ in mpfr_cos()
233 MPFR_SET_EXP (s, MPFR_GET_EXP (s) + 1); /* Can't overflow */ in mpfr_cos()
237 MPFR_ASSERTD (MPFR_GET_EXP (s) <= 1); in mpfr_cos()
251 exps = MPFR_GET_EXP (s); in mpfr_cos()
H A Dagm.c142 e1 = MPFR_GET_EXP (op1); in mpfr_agm()
143 e2 = MPFR_GET_EXP (op2); in mpfr_agm()
242 err = MPFR_GET_EXP (vf) - MPFR_GET_EXP (v); /* 0 or 1 */ in mpfr_agm()
259 scale2 = - (((MPFR_GET_EXP (u) + MPFR_GET_EXP (v)) in mpfr_agm()
H A Dget_f.c88 sh = MPFR_GET_EXP(y) % GMP_NUMB_BITS; in mpfr_get_f()
110 EXP(x) = (MPFR_GET_EXP(y) + sh) / GMP_NUMB_BITS; in mpfr_get_f()
127 sh = MPFR_GET_EXP(z) % GMP_NUMB_BITS; in mpfr_get_f()
140 EXP(x) = (MPFR_GET_EXP(z) + sh) / GMP_NUMB_BITS; in mpfr_get_f()
H A Dlog.c97 else if (MPFR_UNLIKELY (MPFR_GET_EXP (a) == 1 && mpfr_cmp_ui (a, 1) == 0)) in mpfr_log()
123 m = (p + 1) / 2 - MPFR_GET_EXP (a) + 1; in mpfr_log()
137 cancel = MPFR_GET_EXP (tmp2) - MPFR_GET_EXP (tmp1); in mpfr_log()
H A Dzeta.c153 MPFR_ASSERTD (MPFR_IS_POS (s) && MPFR_GET_EXP (s) >= 0); in mpfr_zeta_pos()
167 if (MPFR_GET_EXP (s) > 3) in mpfr_zeta_pos()
170 err = MPFR_GET_EXP (s) - 1; in mpfr_zeta_pos()
183 dint = (mpfr_uexp_t) MPFR_GET_EXP (s); in mpfr_zeta_pos()
198 if (MPFR_GET_EXP (s1) <= -(mpfr_exp_t) ((mpfr_prec_t) (d-3)/2)) in mpfr_zeta_pos()
335 if (MPFR_GET_EXP (s) + 1 < - (mpfr_exp_t) MPFR_PREC(z)) in mpfr_zeta()
371 MPFR_EXP (tmp) = MPFR_GET_EXP (s) - 1; in mpfr_zeta()
392 if (MPFR_IS_POS (s) && MPFR_GET_EXP (s) >= 0) /* Case s >= 1/2 */ in mpfr_zeta()
404 precs1 = precs + 2 + MAX (0, - MPFR_GET_EXP (s)); in mpfr_zeta()

1234