/dragonfly/contrib/mpc/src/ |
H A D | acosh.c | 34 int inex; in mpc_acosh() local 48 inex = mpc_acos (a, op, in mpc_acosh() 56 inex = MPC_INEX (MPC_INEX_IM (inex), -MPC_INEX_RE (inex)); in mpc_acosh() 60 inex = mpc_acos (a, op, in mpc_acosh() 68 inex = MPC_INEX (-MPC_INEX_IM (inex), MPC_INEX_RE (inex)); in mpc_acosh() 75 return inex; in mpc_acosh()
|
H A D | tan.c | 32 int inex; in mpc_tan() local 44 inex = mpc_set_si_si (rop, 0, in mpc_tan() 92 inex = MPC_INEX (0, inex_im); in mpc_tan() 122 inex = MPC_INEX (0, inex_im); in mpc_tan() 128 return inex; in mpc_tan() 228 inex = MPC_INEX(inex_re, inex_im); in mpc_tan() 237 inex = mpc_div (x, x, y, MPC_RNDZZ); in mpc_tan() 249 if (MPC_INEX_RE (inex)) in mpc_tan() 251 if (MPC_INEX_IM (inex)) in mpc_tan() 277 inex = mpc_set (rop, x, rnd); in mpc_tan() [all …]
|
H A D | mpc.h | 51 #define MPC_INEX_POS(inex) (((inex) < 0) ? 2 : ((inex) == 0) ? 0 : 1) argument 53 #define MPC_INEX_NEG(inex) (((inex) == 2) ? -1 : ((inex) == 0) ? 0 : 1) argument 62 #define MPC_INEX_RE(inex) MPC_INEX_NEG((inex) & 3) argument 63 #define MPC_INEX_IM(inex) MPC_INEX_NEG((inex) >> 2) argument 68 #define MPC_INEX1(inex) (inex & 15) argument 69 #define MPC_INEX2(inex) (inex >> 4) argument
|
H A D | pow_ui.c | 28 int inex; in mpc_pow_usi_naive() local 36 inex = mpc_pow (z, x, t, rnd); in mpc_pow_usi_naive() 39 return inex; in mpc_pow_usi_naive() 48 int inex; in mpc_pow_usi() local 117 inex = mpc_pow_usi_naive (z, x, y, sign, rnd); in mpc_pow_usi() 138 inex = mpc_set (z, t, rnd); in mpc_pow_usi() 151 inex = mpc_pow_usi_naive (z, x, y, sign, rnd); in mpc_pow_usi() 161 return inex; in mpc_pow_usi()
|
H A D | set_str.c | 28 int inex; in mpc_set_str() local 30 inex = mpc_strtoc (z, str, &p, base, rnd); in mpc_set_str() 32 if (inex != -1){ in mpc_set_str() 36 return inex; in mpc_set_str()
|
H A D | sqr.c | 33 int inex; in mpfr_fsss() local 44 inex = mpfr_sub (z, u, v, rnd); in mpfr_fsss() 49 inex = mpfr_mul_2ui (z, z, mpfr_get_emax (), rnd); in mpfr_fsss() 53 inex = (mpfr_signbit (u) ? 1 : -1); in mpfr_fsss() 57 inex = (mpfr_signbit (v) ? -1 : 1); in mpfr_fsss() 115 inex = mpfr_sub (z, u, v, rnd); in mpfr_fsss() 119 inex = overflow; in mpfr_fsss() 140 inex = mpfr_sub (z, u, v, rnd); in mpfr_fsss() 144 inex = underflow; in mpfr_fsss() 158 return inex; in mpfr_fsss()
|
H A D | sinh.c | 29 int inex; in mpc_sinh() local 39 inex = mpc_sin (sin_z, z, MPC_RND (MPC_RND_IM (rnd), MPC_RND_RE (rnd))); in mpc_sinh() 46 return MPC_INEX (MPC_INEX_IM (inex), MPC_INEX_RE (inex)); in mpc_sinh()
|
H A D | tanh.c | 29 int inex; in mpc_tanh() local 39 inex = mpc_tan (tan_z, z, MPC_RND (MPC_RND_IM (rnd), MPC_RND_RE (rnd))); in mpc_tanh() 46 return MPC_INEX (MPC_INEX_IM (inex), MPC_INEX_RE (inex)); in mpc_tanh()
|
H A D | asinh.c | 27 int inex; in mpc_asinh() local 39 inex = mpc_asin (a, z, in mpc_asinh() 54 return MPC_INEX (MPC_INEX_IM (inex), -MPC_INEX_RE (inex)); in mpc_asinh()
|
H A D | atanh.c | 27 int inex; in mpc_atanh() local 38 inex = mpc_atan (a, z, in mpc_atanh() 51 return MPC_INEX (MPC_INEX_IM (inex), -MPC_INEX_RE (inex)); in mpc_atanh()
|
H A D | mul.c | 116 int inex; in mul_real() local 134 return inex; in mul_real() 184 int inex; in mpfr_fmma() local 197 inex = mpfr_add (z, u, v, rnd); in mpfr_fmma() 275 inex = mpfr_add (z, u, v, rnd); in mpfr_fmma() 279 inex = overflow; in mpfr_fmma() 300 inex = mpfr_add (z, u, v, rnd); in mpfr_fmma() 304 inex = underflow; in mpfr_fmma() 320 return inex; in mpfr_fmma() 329 int overlap, inex; in mpc_mul_naive() local [all …]
|
/dragonfly/contrib/mpfr/src/ |
H A D | digamma.c | 126 int inex; in mpfr_digamma_reflection() local 144 inex = mpfr_integer_p (u); in mpfr_digamma_reflection() 146 if (inex) in mpfr_digamma_reflection() 190 inex = mpfr_set (y, v, rnd_mode); in mpfr_digamma_reflection() 197 return inex; in mpfr_digamma_reflection() 206 int inex; in mpfr_digamma_positive() local 275 return inex; in mpfr_digamma_positive() 281 int inex; in mpfr_digamma() local 353 inex = 1; in mpfr_digamma() 357 inex = -1; in mpfr_digamma() [all …]
|
H A D | div.c | 164 int inex; in mpfr_div() local 349 inex = -1; in mpfr_div() 360 inex = -1; in mpfr_div() 502 if (inex == 0) in mpfr_div() 694 inex = -1; in mpfr_div() 704 inex = 1; in mpfr_div() 728 inex = -1; in mpfr_div() 760 if (inex < 0) in mpfr_div() 765 inex = 1; /* always here */ in mpfr_div() 788 inex *= sign_quotient; in mpfr_div() [all …]
|
H A D | exceptions.c | 320 int inex; in mpfr_underflow() local 327 inex = -1; in mpfr_underflow() 332 inex = 1; in mpfr_underflow() 336 return sign > 0 ? inex : -inex; in mpfr_underflow() 344 int inex; in mpfr_overflow() local 350 inex = -1; in mpfr_overflow() 355 inex = 1; in mpfr_overflow() 359 return sign > 0 ? inex : -inex; in mpfr_overflow()
|
H A D | erf.c | 34 int inex, large; in mpfr_erf() local 98 inex = mpfr_prec_round (l, MPFR_PREC(y), rnd_mode); in mpfr_erf() 101 ok = SAME_SIGN (inex, inex2) && mpfr_cmp (l, h) == 0; in mpfr_erf() 131 inex = MPFR_INT_SIGN (x); in mpfr_erf() 132 mpfr_set_si (y, inex, rnd2); in mpfr_erf() 136 inex = -MPFR_INT_SIGN (x); in mpfr_erf() 148 inex = mpfr_erf_0 (y, x, xf2, rnd_mode); in mpfr_erf() 153 return mpfr_check_range (y, inex, rnd_mode); in mpfr_erf() 190 int inex; in mpfr_erf_0() local 254 inex = mpfr_set (res, s, rnd_mode); in mpfr_erf_0() [all …]
|
H A D | zeta.c | 149 int inex; in mpfr_zeta_pos() local 186 MPFR_ASSERTD (inex == 0); in mpfr_zeta_pos() 282 inex = mpfr_set (z, z_pre, rnd_mode); in mpfr_zeta_pos() 287 return inex; in mpfr_zeta_pos() 297 int inex; in mpfr_zeta() local 346 inex = 1; in mpfr_zeta() 351 inex = -1; in mpfr_zeta() 356 inex = 1; in mpfr_zeta() 360 inex = (signs > 0) ? 1 : -1; in mpfr_zeta() 393 inex = mpfr_zeta_pos (z, s, rnd_mode); in mpfr_zeta() [all …]
|
H A D | get_z.c | 29 int inex; in mpfr_get_z() local 49 inex = mpfr_rint (r, f, rnd); in mpfr_get_z() 50 MPFR_ASSERTN (inex != 1 && inex != -1); /* integral part of f is in mpfr_get_z() 60 return inex; in mpfr_get_z()
|
H A D | add1.c | 34 int sh, rb, fb, inex; in mpfr_add1() local 164 inex = mpfr_overflow (a, rnd_mode, MPFR_SIGN(a)); in mpfr_add1() 305 inex = mpfr_overflow (a, rnd_mode, MPFR_SIGN(a)); in mpfr_add1() 483 inex = 0; in mpfr_add1() 495 inex = MPFR_IS_NEG(a) ? 1 : -1; in mpfr_add1() 501 inex = MPFR_IS_POS(a) ? 1 : -1; in mpfr_add1() 507 inex = rb || fb ? (MPFR_IS_NEG(a) ? 1 : -1) : 0; in mpfr_add1() 513 inex = rb || fb ? (MPFR_IS_POS(a) ? 1 : -1) : 0; in mpfr_add1() 514 if (inex) in mpfr_add1() 525 inex = mpfr_overflow (a, rnd_mode, MPFR_SIGN(a)); in mpfr_add1() [all …]
|
H A D | gamma.c | 104 int inex = 0; /* 0 means: result gamma not set yet */ in mpfr_gamma() local 112 mpfr_get_prec (gamma), mpfr_log_prec, gamma, inex)); in mpfr_gamma() 184 MPFR_BLOCK (flags, inex = mpfr_ui_div (gamma, 1, x, rnd_mode)); in mpfr_gamma() 185 if (inex == 0) /* x is a power of two */ in mpfr_gamma() 189 inex = 1; in mpfr_gamma() 193 inex = -1; in mpfr_gamma() 209 return mpfr_check_range (gamma, inex, rnd_mode); in mpfr_gamma() 385 if (inex) in mpfr_gamma() 432 if (inex == 0) in mpfr_gamma() 433 inex = mpfr_set (gamma, GammaTrial, rnd_mode); in mpfr_gamma() [all …]
|
H A D | rint.c | 344 int inex; in mpfr_rint_round() local 352 inex = (MPFR_OVERFLOW (flags) in mpfr_rint_round() 357 return mpfr_check_range (r, inex, rnd_mode); in mpfr_rint_round() 371 int inex; in mpfr_rint_trunc() local 378 inex = mpfr_set (r, tmp, rnd_mode); in mpfr_rint_trunc() 381 return mpfr_check_range (r, inex, rnd_mode); in mpfr_rint_trunc() 395 int inex; in mpfr_rint_ceil() local 403 inex = (MPFR_OVERFLOW (flags) in mpfr_rint_ceil() 408 return mpfr_check_range (r, inex, rnd_mode); in mpfr_rint_ceil() 422 int inex; in mpfr_rint_floor() local [all …]
|
H A D | jn.c | 85 int inex; in mpfr_jn() local 100 mpfr_get_prec (res), mpfr_log_prec, res, inex)); in mpfr_jn() 135 2, 0, r, inex = _inexact; goto end); in mpfr_jn() 169 inex = - inex2; in mpfr_jn() 172 inex = inex2 != 0 ? inex2 : _inexact; in mpfr_jn() 185 inex = mpfr_jn_asympt (res, n, z, r); in mpfr_jn() 186 if (inex != 0) in mpfr_jn() 318 inex = ((n >= 0) || ((n & 1) == 0)) ? mpfr_set (res, s, r) in mpfr_jn() 325 return mpfr_check_range (res, inex, r); in mpfr_jn()
|
H A D | sin.c | 29 int inex; in mpfr_sin_fast() local 31 inex = mpfr_sincos_fast (y, NULL, x, rnd_mode); in mpfr_sin_fast() 32 inex = inex & 3; /* 0: exact, 1: rounded up, 2: rounded down */ in mpfr_sin_fast() 33 return (inex == 2) ? -1 : inex; in mpfr_sin_fast()
|
H A D | sub_ui.c | 35 int inex; in mpfr_sub_ui() local 43 mpfr_get_prec(y), mpfr_log_prec, y, inex)); in mpfr_sub_ui() 54 inex = mpfr_sub (y, x, uu, rnd_mode); in mpfr_sub_ui() 56 return mpfr_check_range (y, inex, rnd_mode); in mpfr_sub_ui()
|
H A D | set_z_exp.c | 32 int k, sign_z, inex; in mpfr_set_z_2exp() local 143 inex = MPFR_LIKELY ((sb | rb) != 0) ? -1 : 0; in mpfr_set_z_2exp() 147 inex = 1; in mpfr_set_z_2exp() 167 inex = 0; /* result is exact */ in mpfr_set_z_2exp() 172 if (rnd_mode == MPFR_RNDN && inex == 0 && mpfr_powerof2_raw (f)) in mpfr_set_z_2exp() 179 MPFR_RET (inex*sign_z); in mpfr_set_z_2exp()
|
/dragonfly/contrib/gdtoa/ |
H A D | gdtoa.c | 171 inex = 0; 470 inex = STRTOG_Inexlo; 481 inex = STRTOG_Inexhi; 491 inex = STRTOG_Inexlo; 595 inex = STRTOG_Inexlo; 599 inex = STRTOG_Inexhi; 638 inex = STRTOG_Inexhi; 693 inex = STRTOG_Inexhi; 733 inex = STRTOG_Inexhi; 745 inex = STRTOG_Inexlo; [all …]
|