/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/ |
H A D | atan2.c | 126 mpfr_const_pi (tmp2, MPFR_RNDN); 127 mpfr_mul_ui (tmp2, tmp2, 3, MPFR_RNDN); /* Error <= 2 */ 128 mpfr_div_2ui (tmp2, tmp2, 2, MPFR_RNDN); 131 MPFR_PREC (dest) + (rnd_mode == MPFR_RNDN))) 171 MPFR_BLOCK (flags, div_inex = mpfr_div (tmp, y, x, MPFR_RNDN)); 187 /* In the case MPFR_RNDN with 2^(emin-2) < |y/x| < 2^(emin-1): 196 (this is what mpfr_underflow always does on MPFR_RNDN). 197 In the case MPFR_RNDN with |y/x| <= 2^(emin-2), we round 201 if (rnd_mode == MPFR_RNDN && MPFR_IS_ZERO (tmp)) 209 mpfr_atan (tmp, tmp, MPFR_RNDN); /* Erro [all...] |
H A D | cbrt.c | 90 n = MPFR_PREC (y) + (rnd_mode == MPFR_RNDN); 127 || (rnd_mode == MPFR_RNDN && mpz_tstbit (m, 0))) 135 rnd_mode=MPFR_RNDN and bit (n+1) from m is 1 */ 136 inexact += mpfr_set_z (y, m, MPFR_RNDN);
|
H A D | urandom.c | 61 || (emin == 1 && rnd_mode == MPFR_RNDN 96 || (rnd_mode == MPFR_RNDN && cnt == exp - emin - 1 128 || (rnd_mode == MPFR_RNDN && random_rounding_bit (rstate)))
|
H A D | asinh.c | 92 mpfr_sqrt (t, t, MPFR_RNDN); /* sqrt(x^2+1) */ 93 (neg ? mpfr_sub : mpfr_add) (t, t, x, MPFR_RNDN); /* sqrt(x^2+1)+x */ 94 mpfr_log (t, t, MPFR_RNDN); /* ln(sqrt(x^2+1)+x)*/
|
H A D | atanh.c | 99 mpfr_div (t, t, te, MPFR_RNDN); /* (1+xt)/(1-xt)*/ 100 mpfr_log (t, t, MPFR_RNDN); /* ln((1+xt)/(1-xt))*/ 101 mpfr_div_2ui (t, t, 1, MPFR_RNDN); /* (1/2)*ln((1+xt)/(1-xt))*/
|
H A D | const_catalan.c | 128 mpfr_mul (x, x, y, MPFR_RNDN); 133 mpfr_div (y, y, z, MPFR_RNDN); 134 mpfr_add (x, x, y, MPFR_RNDN); 135 mpfr_div_2ui (x, x, 3, MPFR_RNDN);
|
H A D | tanh.c | 105 mpfr_mul_2ui (te, x, 1, MPFR_RNDN); /* 2x */ 107 mpfr_exp (te, te, MPFR_RNDN); /* exp(2x) */ 111 mpfr_set4 (y, __gmpfr_one, MPFR_RNDN, sign); 123 mpfr_div (t, te, t, MPFR_RNDN); /* (exp(2x)-1)/(exp(2x)+1)*/
|
H A D | ui_pow.c | 34 inexact = mpfr_set_ui (t, n, MPFR_RNDN);
|
H A D | erf.c | 48 return mpfr_set_si (y, MPFR_INT_SIGN (x), MPFR_RNDN); 52 return mpfr_set (y, x, MPFR_RNDN); /* should keep the sign of x */ 128 if (rnd2 == MPFR_RNDN || rnd2 == MPFR_RNDU || rnd2 == MPFR_RNDA) 145 xf2 = mpfr_get_d (x, MPFR_RNDN); 206 mpfr_set_ui (s, 1, MPFR_RNDN); 207 mpfr_set_ui (t, 1, MPFR_RNDN); 217 mpfr_sub (s, s, u, MPFR_RNDN); 219 mpfr_add (s, s, u, MPFR_RNDN); 236 mpfr_div (s, s, t, MPFR_RNDN);
|
H A D | fma.c | 109 if (MPFR_UNLIKELY (mpfr_mul (u, x, y, MPFR_RNDN))) 114 It is an overflow iff u is an infinity (since MPFR_RNDN was used). 137 inexact = mpfr_div_2ui (u, x, 2, MPFR_RNDN); 139 inexact = mpfr_mul (u, u, y, MPFR_RNDN); 220 inexact = mpfr_mul_2ui (scaled_z, z, scale, MPFR_RNDN); 228 mpfr_mul_2ui (scaled_v, x, scale, MPFR_RNDN); 229 mpfr_mul (u, scaled_v, y, MPFR_RNDN); 234 mpfr_mul_2ui (scaled_v, y, scale, MPFR_RNDN); 235 mpfr_mul (u, x, scaled_v, MPFR_RNDN); 271 inex2 = mpfr_div_2ui (s, s, scale, MPFR_RNDN); [all...] |
H A D | fms.c | 111 if (MPFR_UNLIKELY (mpfr_mul (u, x, y, MPFR_RNDN))) 116 It is an overflow iff u is an infinity (since MPFR_RNDN was used). 139 inexact = mpfr_div_2ui (u, x, 2, MPFR_RNDN); 141 inexact = mpfr_mul (u, u, y, MPFR_RNDN); 222 inexact = mpfr_mul_2ui (scaled_z, z, scale, MPFR_RNDN); 230 mpfr_mul_2ui (scaled_v, x, scale, MPFR_RNDN); 231 mpfr_mul (u, scaled_v, y, MPFR_RNDN); 236 mpfr_mul_2ui (scaled_v, y, scale, MPFR_RNDN); 237 mpfr_mul (u, x, scaled_v, MPFR_RNDN); 273 inex2 = mpfr_div_2ui (s, s, scale, MPFR_RNDN); [all...] |
H A D | coth.c | 45 mpfr_sub_si (z, z, MPFR_SIGN(z) > 0 ? 1 : -1, MPFR_RNDN); \ 48 mpfr_add_si (z, z, MPFR_SIGN(z) > 0 ? 1 : -1, MPFR_RNDN); \
|
H A D | gen_inverse.h | 39 #define ACTION_ZERO(y) return mpfr_set_ui (y, 1, MPFR_RNDN) 85 return mpfr_underflow (y, (rnd_mode == MPFR_RNDN) ? 88 mpfr_ui_div (z, 1, z, MPFR_RNDN);
|
H A D | round_near_x.c | 173 MPFR_PREC (y) + (rnd == MPFR_RNDN))))) 217 if (rnd != MPFR_RNDN && !MPFR_IS_LIKE_RNDZ (rnd, MPFR_IS_NEG_SIGN(sign)))
|
/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/tests/ |
H A D | tmul_2exp.c | 69 inex1 = mpfr_set_ui_2exp (x, i, -4, MPFR_RNDN); 78 inex1 = mpfr_mul_2si (y, x, e, MPFR_RNDN); 80 inex1 = mpfr_set_ui (z1, 1 << k, MPFR_RNDN); 111 mpfr_out_str (stdout, 16, 0, z1, MPFR_RNDN); 114 mpfr_out_str (stdout, 16, 0, z2, MPFR_RNDN); 274 mpfr_set_str (w, val[k], 16, MPFR_RNDN); 276 if (mpfr_cmp_str (z, val[k+1], 16, MPFR_RNDN)) 281 mpfr_out_str (stdout, 16, 0, z, MPFR_RNDN); 286 if (mpfr_cmp_str (z, val[k+2], 16, MPFR_RNDN)) 291 mpfr_out_str (stdout, 16, 0, z, MPFR_RNDN); [all...] |
H A D | tsum.c | 110 mpfr_set_ui (somme, 0, MPFR_RNDN); 113 if (mpfr_add(somme, somme, tab[i], MPFR_RNDN)) 232 i = mpfr_sum (r, tabp, 0, MPFR_RNDN); 239 mpfr_set_ui (tab[0], 42, MPFR_RNDN); 240 i = mpfr_sum (r, tabp, 1, MPFR_RNDN); 247 mpfr_set_ui (tab[1], 17, MPFR_RNDN); 249 i = mpfr_sum (r, tabp, 3, MPFR_RNDN); 258 i = mpfr_sum (r, tabp, 3, MPFR_RNDN); 267 i = mpfr_sum (r, tabp, 3, MPFR_RNDN); 275 i = mpfr_sum (r, tabp, 2, MPFR_RNDN); [all...] |
H A D | tpow_all.c | 59 mpfr_out_str (stdout, 10, 0, z, MPFR_RNDN); 105 mpfr_out_str (stdout, 16, 0, (mpfr_ptr) px, MPFR_RNDN); 116 mpfr_out_str (stdout, 16, 0, z1, MPFR_RNDN); 120 mpfr_out_str (stdout, 16, 0, z2, MPFR_RNDN); 132 return mpfr_integer_p (x) && mpfr_fits_slong_p (x, MPFR_RNDN) && 133 (mpfr_get_si (x, MPFR_RNDN) & 1); 164 mpfr_fits_ulong_p (y, MPFR_RNDN)) 166 unsigned long yy = mpfr_get_ui (y, MPFR_RNDN); 180 mpfr_fits_ulong_p (x, MPFR_RNDN)) 182 unsigned long xx = mpfr_get_ui (x, MPFR_RNDN); [all...] |
H A D | tcan_round.c | 48 MPFR_RNDN, MPFR_RNDZ, p); 87 if (mpfr_can_round (x, 65, MPFR_RNDN, MPFR_RNDN, 54))
|
H A D | tset_q.c | 47 mpfr_out_str(stdout, 10, 0, x, MPFR_RNDN); 110 check(632549085, 1831935802, MPFR_RNDN, "3.4528998467600230393e-1"); 111 check (1, 1, MPFR_RNDN, "1.0");
|
H A D | tcheck.c | 53 mpfr_set_ui(a, 2, MPFR_RNDN); 59 mpfr_mul(a, a, a, MPFR_RNDN); 64 mpfr_set_ui(a, 2137, MPFR_RNDN); 116 mpfr_set_ui(a, 2137, MPFR_RNDN);
|
H A D | tcmpabs.c | 87 mpfr_set_ui(xx, 0, MPFR_RNDN); 88 mpfr_set_str (yy, "-0.1", 10, MPFR_RNDN); 93 mpfr_set_str (yy, "23489745.0329", 10, MPFR_RNDN); 103 mpfr_set_str (xx, "2346.09234", 10, MPFR_RNDN);
|
H A D | tdigamma.c | 39 mpfr_digamma (y, x, MPFR_RNDN);
|
H A D | tgeneric_ui.c | 64 mpfr_set_si (x, n == 0 ? 1 : -1, MPFR_RNDN); 78 mpfr_out_str (stdout, 2, prec, x, MPFR_RNDN); 86 mpfr_out_str (stdout, 2, prec, z, MPFR_RNDN); 89 mpfr_out_str (stdout, 2, prec, t, MPFR_RNDN);
|
H A D | tgmpop.c | 55 mpfr_add_q (y, x, z, MPFR_RNDN); 99 res = mpfr_add_q (y, x, z, MPFR_RNDN); 101 res = mpfr_sub_q (y, x, z, MPFR_RNDN); 104 res = mpfr_add_q (y, x, z, MPFR_RNDN); 106 res = mpfr_sub_q (y, x, z, MPFR_RNDN); 108 res = mpfr_div_q (y, x, z, MPFR_RNDN); 113 mpfr_set_ui (x, 42, MPFR_RNDN); 114 res = mpfr_add_q (y, x, z, MPFR_RNDN); 116 res = mpfr_sub_q (y, x, z, MPFR_RNDN); 118 res = mpfr_mul_q (y, x, z, MPFR_RNDN); [all...] |
H A D | tgeneric.c | 53 mpfr_out_str (stdout, 2, 0, (X), MPFR_RNDN); \ 58 mpfr_out_str (stdout, 2, 0, (U), MPFR_RNDN); \ 85 mpfr_out_str (stdout, 2, 0, (X), MPFR_RNDN); \ 90 mpfr_out_str (stdout, 2, 0, (U), MPFR_RNDN); \ 185 mpfr_set_si (x, n == 0 ? 1 : -1, MPFR_RNDN); 188 mpfr_set_si (u, randlimb () % 2 == 0 ? 1 : -1, MPFR_RNDN); 196 mpfr_set_si (x, n == 0 ? 1 : -1, MPFR_RNDN); 199 mpfr_set_si (u, randlimb () % 2 == 0 ? 1 : -1, MPFR_RNDN); 352 MPFR_RNDN); 355 MPFR_RNDN); [all...] |