/dragonfly/contrib/mpfr/src/ |
H A D | li2.c | 65 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 D | eint.c | 51 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 D | erfc.c | 60 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 D | ai.c | 155 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 D | gamma.c | 42 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 D | subnormal.c | 53 || (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 D | fms.c | 135 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 D | fma.c | 133 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 D | lngamma.c | 123 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 D | tanh.c | 62 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 D | pow.c | 45 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 D | sin.c | 70 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 D | cmp2.c | 54 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 D | acosh.c | 106 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 D | jn.c | 134 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 D | sum.c | 77 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 D | sinh.c | 63 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 D | add.c | 101 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 D | atanh.c | 64 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 D | asin.c | 58 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 D | cos.c | 163 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 D | agm.c | 142 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 D | get_f.c | 88 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 D | log.c | 97 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 D | zeta.c | 153 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()
|