/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/tests/ |
H A D | texpm1.c | 33 int ok = rnd_mode == MPFR_RNDN && mpfr_number_p (b) && mpfr_get_prec (a)>=53; 66 test_expm1 (y, x, MPFR_RNDN); 74 test_expm1 (y, x, MPFR_RNDN); 82 test_expm1 (y, x, MPFR_RNDN); 89 mpfr_set_ui (x, 0, MPFR_RNDN); 90 test_expm1 (y, x, MPFR_RNDN); 97 mpfr_neg (x, x, MPFR_RNDN); 98 test_expm1 (y, x, MPFR_RNDN); 108 i = test_expm1 (x, x, MPFR_RNDN); 131 i = test_expm1 (x, x, MPFR_RNDN); [all...] |
H A D | tmodf.c | 48 mpfr_out_str (stdout, 10, 0, xi, MPFR_RNDN); 57 mpfr_out_str (stdout, 10, 0, xf, MPFR_RNDN); 82 mpfr_modf (xi, xf, x, MPFR_RNDN); 88 mpfr_modf (xi, xf, x, MPFR_RNDN); 95 mpfr_modf (xi ,xf, x, MPFR_RNDN); 118 mpfr_set_str (x, "0.11111111", 2, MPFR_RNDN); 126 inexact = mpfr_modf (xi, xf, x, MPFR_RNDN); 134 mpfr_out_str (stdout, 2, 0, xi, MPFR_RNDN); 136 mpfr_out_str (stdout, 2, 0, xf, MPFR_RNDN); 160 inexact = mpfr_modf (xi, xf, x, MPFR_RNDN); [all...] |
H A D | tadd_ui.c | 41 mpfr_out_str(stdout, 10, 0, zz, MPFR_RNDN); 72 mpfr_add_ui (y, x, 2394875L, MPFR_RNDN); 77 mpfr_add_ui (y, x, 2394875L, MPFR_RNDN); 83 mpfr_add_ui (y, x, 2394875L, MPFR_RNDN); 106 check3 ("1.22191250737771397120e+20", 948002822, MPFR_RNDN, 110 check3 ("-2.0769715792901673e-5", 880524, MPFR_RNDN,
|
H A D | tconst_euler.c | 35 #define RAND_FUNCTION(x) mpfr_set_ui ((x), 0, MPFR_RNDN) 53 mpfr_const_euler (gamma, MPFR_RNDN); 66 (mpfr_const_euler) (y, MPFR_RNDN); 84 err = (rnd == MPFR_RNDN) ? yprec + 1 : yprec; 94 mpfr_out_str (stdout, 2, prec, z, MPFR_RNDN); 97 mpfr_out_str (stdout, 2, prec, t, MPFR_RNDN);
|
H A D | tsqr.c | 94 mpfr_sqr (y, x, MPFR_RNDN); 98 mpfr_sqr (y, x, MPFR_RNDN); 102 mpfr_sqr (y, x, MPFR_RNDN); 105 mpfr_set_ui (x, 0, MPFR_RNDN); 106 mpfr_sqr (y, x, MPFR_RNDN); 111 mpfr_set_ui (x, 1, MPFR_RNDN); 112 mpfr_div_2ui (x, x, 1, MPFR_RNDN); 114 mpfr_sqr (y, x, MPFR_RNDN);
|
H A D | tgamma.c | 45 mpfr_gamma (y, x, MPFR_RNDN); 53 mpfr_gamma (y, x, MPFR_RNDN); 61 mpfr_gamma (y, x, MPFR_RNDN); 68 mpfr_set_ui (x, 0, MPFR_RNDN); 69 mpfr_gamma (y, x, MPFR_RNDN); 76 mpfr_set_ui (x, 0, MPFR_RNDN); 77 mpfr_neg (x, x, MPFR_RNDN); 78 mpfr_gamma (y, x, MPFR_RNDN); 85 mpfr_set_ui (x, 1, MPFR_RNDN); 86 mpfr_gamma (y, x, MPFR_RNDN); [all...] |
H A D | tset_str.c | 59 printf("ERROR for mpfr_set_str (a, \"0.00000000001\", 10, MPFR_RNDN)\n" 70 mpfr_set_str (a, "2E1000000000", 10, MPFR_RNDN); 73 printf("ERROR for mpfr_set_str (a, \"2E1000000000\", 10, MPFR_RNDN);\n"); 116 mpfr_set_str (x, argv[1], base, MPFR_RNDN); 117 mpfr_out_str (stdout, obase, 0, x, MPFR_RNDN); 155 mpfr_set_str (y, "4.936a52bc17254@-133", 16, MPFR_RNDN); 169 mpfr_set_str (y, "0.fedcba98765434P-529", 16, MPFR_RNDN); 207 mpfr_set_str (x, "+243495834958.53452345E1", 10, MPFR_RNDN); 208 mpfr_set_str (x, "9007199254740993", 10, MPFR_RNDN); 256 if (mpfr_set_str (x, "@NaN@(garbage)", i, MPFR_RNDN) ! [all...] |
H A D | tget_str.c | 58 mpfr_set_ui (x, 2, MPFR_RNDN); 71 mpfr_set_ui (x, 17, MPFR_RNDN); 72 s = mpfr_get_str (NULL, &e, 3, 0, x, MPFR_RNDN); 74 s = mpfr_get_str (NULL, &e, 36, 0, x, MPFR_RNDN); 76 s = mpfr_get_str (NULL, &e, 62, 0, x, MPFR_RNDN); 80 mpfr_set_si (x, -1, MPFR_RNDN); 81 mpfr_div_2exp (x, x, 63, MPFR_RNDN); /* x = -2^(-63) */ 82 mpfr_add_ui (x, x, 1, MPFR_RNDN); /* x = 1 - 2^(-63) */ 83 mpfr_mul_2exp (x, x, 32, MPFR_RNDN); /* x = 2^32 - 2^(-31) */ 128 mpfr_set_str (x, "6.5", 10, MPFR_RNDN); [all...] |
H A D | teq.c | 41 mpfr_set (y, x, MPFR_RNDN); 105 mpfr_set_ui (y, 0, MPFR_RNDN); 121 mpfr_set_ui (y, 0, MPFR_RNDN); 124 mpfr_set_ui (x, 1, MPFR_RNDN); 125 mpfr_set_ui (y, 0, MPFR_RNDN); 129 mpfr_set_ui (x, 0, MPFR_RNDN); 130 mpfr_set_ui (y, 0, MPFR_RNDN); 133 mpfr_neg (y, y, MPFR_RNDN); 136 mpfr_neg (x, x, MPFR_RNDN); 139 mpfr_set_ui (x, 1, MPFR_RNDN); [all...] |
H A D | tpow_z.c | 43 mpfr_set_ui (x, 23, MPFR_RNDN); 45 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 49 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 53 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 60 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 64 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 69 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 73 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 79 res = mpfr_pow_z (y, x, z, MPFR_RNDN); 83 res = mpfr_pow_z (y, x, z, MPFR_RNDN); [all...] |
H A D | tcmp2.c | 59 mpfr_set_ui (y, 1, MPFR_RNDN); 62 mpfr_set_ui (x, 1, MPFR_RNDN); 63 mpfr_div_2exp (y, y, 1, MPFR_RNDN); /* y = 1/2^i */ 69 mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN); 71 mpfr_out_str (stdout, 2, 0, y, MPFR_RNDN); 76 mpfr_add (x, x, y, MPFR_RNDN); /* x = 1 + 1/2^i */ 81 mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN); 83 mpfr_out_str (stdout, 2, 0, y, MPFR_RNDN); 92 mpfr_set (y, x, MPFR_RNDN); 113 mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN); [all...] |
H A D | rnd_mode.c | 35 case MPFR_RNDN:
|
H A D | tmul_ui.c | 42 if (mpfr_mul_ui (z, x, u, MPFR_RNDN)) 66 mpfr_set_ui (x, 1, MPFR_RNDN); 151 mpfr_mul_ui (x, x, 2, MPFR_RNDN); 179 mpfr_mul_ui (x, y, 1335, MPFR_RNDN); 204 mpfr_mul_ui (y, x, 1, MPFR_RNDN); 209 mpfr_mul_ui (y, x, 6, MPFR_RNDN); 219 mpfr_set_str (x, "2143861251406875.0", 10, MPFR_RNDN); 220 mpfr_mul_ui (y, x, 23, MPFR_RNDN); 238 mpfr_mul_ui (y, x, 1, MPFR_RNDN); 251 mpfr_set_ui (x, 17, MPFR_RNDN); [all...] |
H A D | tcomparisons.c | 53 mpfr_neg (x, x, MPFR_RNDN); 55 mpfr_neg (y, y, MPFR_RNDN); 81 mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN); 83 mpfr_out_str (stdout, 2, 0, y, MPFR_RNDN); 107 mpfr_neg (x, x, MPFR_RNDN); 108 mpfr_set (y, x, MPFR_RNDN); /* exact -> x = y */ 115 mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN);
|
H A D | tconst_pi.c | 41 mpfr_const_pi (x, MPFR_RNDN); /* First one ! */ 42 mpfr_const_pi (y, MPFR_RNDN); /* Then the other - cache - */ 43 mpfr_prec_round (y, 20000, MPFR_RNDN); 70 #define RAND_FUNCTION(x) mpfr_set_ui ((x), 0, MPFR_RNDN) 105 mpfr_const_pi (x, MPFR_RNDN); 161 mpfr_out_str (stdout, 10, 0, x, MPFR_RNDN); putchar('\n'); 166 mpfr_const_pi (x, MPFR_RNDN);
|
H A D | tpow3.c | 67 mpfr_neg (s, s, MPFR_RNDN); 71 err = (rnd == MPFR_RNDN) ? yprec + 1 : yprec; 79 mpfr_out_str (stdout, 2, prec, x, MPFR_RNDN); 81 mpfr_out_str (stdout, 2, 0, s, MPFR_RNDN); 85 mpfr_out_str (stdout, 2, prec, z, MPFR_RNDN); 88 mpfr_out_str (stdout, 2, prec, t, MPFR_RNDN); 98 if ((rnd != MPFR_RNDN) && (compare * compare2 >= 0))
|
H A D | trint.c | 41 mpfr_rint (y, x, MPFR_RNDN); 45 mpfr_rint (y, x, MPFR_RNDN); 49 mpfr_rint (y, x, MPFR_RNDN); 52 mpfr_set_ui (x, 0, MPFR_RNDN); 53 mpfr_rint (y, x, MPFR_RNDN); 56 mpfr_set_ui (x, 0, MPFR_RNDN); 57 mpfr_neg (x, x, MPFR_RNDN); 58 mpfr_rint (y, x, MPFR_RNDN); 63 mpfr_set_ui (x, 1, MPFR_RNDN); 64 mpfr_mul_2exp (x, x, mp_bits_per_limb, MPFR_RNDN); [all...] |
/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/ |
H A D | yn.c | 55 mpfr_set_ui (s, 1, MPFR_RNDN); /* a[n] */ 60 mpfr_mul (s, s, y, MPFR_RNDN); 64 mpfr_add_z (s, s, f, MPFR_RNDN); 71 mpfr_div_z (s, s, f, MPFR_RNDN); 107 mpfr_fac_ui (t, n, MPFR_RNDN); 108 mpfr_div (t, c, t, MPFR_RNDN); /* c/n! */ 109 mpfr_mul_z (u, t, p, MPFR_RNDN); 110 mpfr_div_z (s, u, q, MPFR_RNDN); 116 mpfr_mul (t, t, y, MPFR_RNDN); 117 mpfr_div_ui (t, t, k, MPFR_RNDN); [all...] |
H A D | log10.c | 110 mpfr_set_ui (t, 10, MPFR_RNDN); /* 10 */ 112 mpfr_log (tt, a, MPFR_RNDN); /* log(a) */ 113 mpfr_div (t, tt, t, MPFR_RNDN); /* log(a)/log(10) */ 124 && mpfr_integer_p (t) && mpfr_fits_ulong_p (t, MPFR_RNDN) 125 && !mpfr_ui_pow_ui (tt, 10, mpfr_get_ui (t, MPFR_RNDN), MPFR_RNDN)
|
H A D | asin.c | 61 inexact = mpfr_abs (xp, x, MPFR_RNDN); 105 mpfr_sqr (xp, x, MPFR_RNDN); 106 mpfr_ui_sub (xp, 1, xp, MPFR_RNDN); 107 mpfr_sqrt (xp, xp, MPFR_RNDN); 108 mpfr_div (xp, x, xp, MPFR_RNDN); 109 mpfr_atan (xp, xp, MPFR_RNDN);
|
H A D | acosh.c | 103 mpfr_log (t, x, MPFR_RNDN); /* err(log) < 1/2 ulp(t) */ 107 mpfr_const_log2 (ln2, MPFR_RNDN); /* err(ln2) < 1/2 ulp(t) */ 108 mpfr_add (t, t, ln2, MPFR_RNDN); /* err <= 3/2 ulp(t) */ 122 mpfr_mul_2ui (t, t, 1, MPFR_RNDN); /* 2t */ 123 mpfr_sqrt (t, t, MPFR_RNDN); /* sqrt(2t) */ 129 mpfr_sqrt (t, t, MPFR_RNDN); /* sqrt(x^2-1) */ 130 mpfr_add (t, t, x, MPFR_RNDN); /* sqrt(x^2-1)+x */ 131 mpfr_log (t, t, MPFR_RNDN); /* ln(sqrt(x^2-1)+x) */
|
H A D | rem1.c | 33 a fmod (MPFR_RNDZ) or a remainder/remquo (MPFR_RNDN). 63 MPFR_ASSERTD (rnd_q == MPFR_RNDN || rnd_q == MPFR_RNDZ); 120 if (rnd_q == MPFR_RNDN) 137 else if (rnd_q == MPFR_RNDN) /* remainder case */ 158 else if (rnd_q == MPFR_RNDN) /* now 0 <= r < 2Y in the remainder case */ 172 inex = mpfr_set_ui (rem, 0, MPFR_RNDN); 175 mpfr_neg (rem, rem, MPFR_RNDN); 179 if (rnd_q == MPFR_RNDN) 187 ((rnd_q == MPFR_RNDN) && (compare == 0) && q_is_odd)); 192 if (quo && (rnd_q == MPFR_RNDN)) [all...] |
H A D | cmp_d.c | 32 res = mpfr_set_d (tmp, d, MPFR_RNDN);
|
H A D | cmp_ld.c | 32 res = mpfr_set_ld (tmp, d, MPFR_RNDN);
|
H A D | ai.c | 180 mpfr_set_ui (ti, 9, MPFR_RNDN); 181 mpfr_cbrt (ti, ti, MPFR_RNDN); 182 mpfr_mul (ti, ti, temp2, MPFR_RNDN); 183 mpfr_ui_div (ti, 1, ti , MPFR_RNDN); /* ti = 1/( Gamma (2/3)*9^(1/3) ) */ 185 mpfr_set_ui (tip1, 3, MPFR_RNDN); 186 mpfr_cbrt (tip1, tip1, MPFR_RNDN); 187 mpfr_mul (tip1, tip1, temp1, MPFR_RNDN); 188 mpfr_neg (tip1, tip1, MPFR_RNDN); 189 mpfr_div (tip1, x, tip1, MPFR_RNDN); /* tip1 = -x/(Gamma (1/3)*3^(1/3)) */ 191 mpfr_add (s, ti, tip1, MPFR_RNDN); [all...] |