| /src/external/lgpl3/mpfr/dist/src/ |
| print_rnd_mode.c | 38 case MPFR_RNDU:
39 return "MPFR_RNDU";
|
| erfc.c | 47 mpfr_ui_div (xx, 1, xx, MPFR_RNDU); /* upper bound for 1/(2x^2), err <= 2 */
48 mpfr_div_2ui (xx, xx, 1, MPFR_RNDU); /* exact */
54 mpfr_mul_ui (t, t, 2 * k - 1, MPFR_RNDU); /* err <= 4k-3 */
55 mpfr_mul (t, t, xx, MPFR_RNDU); /* err <= 4k */
60 mpfr_mul_2si (err, err, MPFR_GET_EXP (y) - MPFR_GET_EXP (t), MPFR_RNDU);
61 mpfr_add_ui (err, err, 14 * k, MPFR_RNDU); /* 2^(1-p) * t <= 2 ulp(t) */
62 mpfr_div_2si (err, err, MPFR_GET_EXP (y) - MPFR_GET_EXP (t), MPFR_RNDU);
66 mpfr_add_ui (err, err, 1, MPFR_RNDU);
75 mpfr_sqr (t, x, MPFR_RNDU); /* rel. err <= 2^(1-p) */
76 mpfr_div_2ui (err, err, 3, MPFR_RNDU); /* err/8 *
[all...] |
| eint.c | 107 mpfr_mul_2ui (eps, eps, w - 1, MPFR_RNDU);
109 mpfr_add_z (eps, eps, t, MPFR_RNDU);
111 mpfr_sub_z (eps, eps, t, MPFR_RNDU);
113 mpfr_mul_2si (eps, eps, sizeinbase - (w - 1) + e, MPFR_RNDU);
114 mpfr_div_ui (eps, eps, k, MPFR_RNDU);
115 mpfr_add_ui (eps, eps, 1, MPFR_RNDU);
125 mpfr_div_ui (erru, eps, k, MPFR_RNDU);
126 mpfr_add_ui (erru, erru, 1, MPFR_RNDU);
128 mpfr_add (errs, errs, erru, MPFR_RNDU);
140 mpfr_add_z (eps, eps, t, MPFR_RNDU);
[all...] |
| erf.c | 81 mpfr_sqr (l, x, MPFR_RNDU);
82 mpfr_div_ui (l, l, 3, MPFR_RNDU); /* upper bound on x^2/3 */
84 mpfr_const_pi (h, MPFR_RNDU); /* upper bound of Pi */
85 mpfr_sqrt (h, h, MPFR_RNDU); /* upper bound on sqrt(Pi) */
95 mpfr_div (h, x, h, MPFR_IS_POS(x) ? MPFR_RNDU : MPFR_RNDD);
126 if (rnd2 == MPFR_RNDN || rnd2 == MPFR_RNDU || rnd2 == MPFR_RNDA)
208 mpfr_sqr (y, x, MPFR_RNDU); /* err <= 1 ulp */
212 mpfr_set_ui (tauk, 0, MPFR_RNDU);
216 mpfr_mul (t, y, t, MPFR_RNDU);
217 mpfr_div_ui (t, t, k, MPFR_RNDU);
[all...] |
| jyn_asympt.c | 58 mpfr_set_ui (c, n, MPFR_RNDU);
59 mpfr_mul_ui (c, c, n, MPFR_RNDU);
60 mpfr_div_2ui (c, c, 1, MPFR_RNDU);
122 mpfr_mul_ui (err_t, t, 6 * k, MPFR_IS_POS(t) ? MPFR_RNDU : MPFR_RNDD);
125 mpfr_abs (err_u, t, MPFR_RNDU);
126 mpfr_mul_2ui (err_u, err_u, w, MPFR_RNDU); /* t * 2^w */
127 mpfr_add (err_u, err_u, err_t, MPFR_RNDU); /* max|t| * 2^w */
131 mpfr_div_2ui (err_s, err_s, w, MPFR_RNDU);
133 mpfr_add (err_s, err_s, t, MPFR_RNDU);
135 mpfr_sub (err_s, err_s, t, MPFR_RNDU);
[all...] |
| csc.c | 56 rnd_mode = (signx > 0) ? MPFR_RNDU : MPFR_RNDD; \
57 if (rnd_mode == MPFR_RNDU) \
|
| csch.c | 59 rnd_mode = (signx > 0) ? MPFR_RNDU : MPFR_RNDD; \
60 if (rnd_mode == MPFR_RNDU || (rnd_mode == MPFR_RNDZ && signx < 0)) \
|
| get_flt.c | 53 rnd_mode = negative ? MPFR_RNDD : MPFR_RNDU;
73 (rnd_mode == MPFR_RNDU ||
84 (rnd_mode == MPFR_RNDZ || rnd_mode == MPFR_RNDU ?
|
| bernoulli.c | 83 mpfr_set_ui_2exp (z, 251469612, -32, MPFR_RNDU); /* 1/e/2/pi <= z */
84 mpfr_mul_ui (z, z, n, MPFR_RNDU);
85 mpfr_log2 (z, z, MPFR_RNDU);
86 mpfr_mul_ui (z, z, n, MPFR_RNDU);
87 p = mpfr_get_ui (z, MPFR_RNDU); /* (n/e/2/pi)^n <= 2^p */
147 mpfr_const_pi (y, MPFR_RNDU);
149 mpfr_mul_2ui (y, y, 1, MPFR_RNDU);
151 mpfr_pow_ui (y, y, n, MPFR_RNDU);
180 mpfr_get_z (num, z, MPFR_RNDU);
|
| ai.c | 137 mpfr_abs (tmp_sp, x, MPFR_RNDU);
138 mpfr_pow_ui (tmp_sp, tmp_sp, 3, MPFR_RNDU);
139 mpfr_sqrt (tmp_sp, tmp_sp, MPFR_RNDU); /* tmp_sp ~ |x|^(3/2) */
142 mpfr_set_str (tmp2_sp, "0.96179669392597567", 10, MPFR_RNDU);
143 mpfr_mul (tmp2_sp, tmp_sp, tmp2_sp, MPFR_RNDU);
152 MPFR_BLOCK (flags, cond = mpfr_get_ui (tmp2_sp, MPFR_RNDU));
169 MPFR_BLOCK (flags, t = mpfr_get_ui (tmp2_sp, MPFR_RNDU));
207 mpfr_sqr (x3, x, MPFR_RNDU);
208 mpfr_mul (x3, x3, x, (MPFR_IS_POS (x)?MPFR_RNDU:MPFR_RNDD)); /* x3=x^3 */
211 x3u = mpfr_get_ui (x3, MPFR_RNDU); /* x3u >= ceil(x^3) *
[all...] |
| zeta.c | 321 mpfr_sub_ui (t, s, 1, MPFR_RNDU);
323 mpfr_ui_sub (t, 1, s, MPFR_RNDU);
325 mpfr_set_ui_2exp (u, 1, - precz - 14, MPFR_RNDU);
330 mpfr_mul_2ui (t, t, 1, MPFR_RNDU);
332 mpfr_set_ui_2exp (t, 1, precz + 14, MPFR_RNDU);
334 mpfr_add_ui (u, u, 1, MPFR_RNDU); /* u = 1+eps, rounded up */
335 mpfr_mul (t, t, u, MPFR_RNDU); /* t = max(1/eps,2*sd)*(1+eps) */
336 mpfr_add_ui (m1, t, 1, MPFR_RNDU);
338 mpfr_set_ui (t, 8, MPFR_RNDU);
340 mpfr_set (t, m1, MPFR_RNDU);
[all...] |
| compound.c | 34 || (s < 0 && (rnd_mode == MPFR_RNDA || rnd_mode == MPFR_RNDU)))
150 rnd1 = VSIGN (n) == MPFR_SIGN (x) ? MPFR_RNDD : MPFR_RNDU;
210 rnd2 = MPFR_IS_POS (u) ? MPFR_RNDD : MPFR_RNDU;
289 mpfr_ui_div (t, 1, x, MPFR_RNDU);
290 mpfr_add_ui (t, t, 1, MPFR_RNDU);
291 mpfr_pow_ui (t, t, n, MPFR_RNDU);
292 mpfr_sub_ui (t, t, 1, MPFR_RNDU);
298 || rnd_mode == MPFR_RNDU || rnd_mode == MPFR_RNDA)
|
| cot.c | 71 rnd_mode = (signx > 0) ? MPFR_RNDU : MPFR_RNDD; \
72 if (rnd_mode == MPFR_RNDU || (rnd_mode == MPFR_RNDZ && signx < 0)) \
|
| coth.c | 73 rnd_mode = (signx > 0) ? MPFR_RNDU : MPFR_RNDD; \
74 if (rnd_mode == MPFR_RNDU) \
|
| const_catalan.c | 125 mpfr_sqrt_ui (x, 3, MPFR_RNDU);
126 mpfr_add_ui (x, x, 2, MPFR_RNDU);
127 mpfr_log (x, x, MPFR_RNDU);
128 mpfr_const_pi (y, MPFR_RNDU);
132 mpfr_set_z (y, T, MPFR_RNDU);
|
| /src/external/lgpl3/mpfr/dist/tests/ |
| tlog.c | 94 check2("1.00089971802309629645", MPFR_RNDU, "8.99313519443722844508e-04");
98 check2("1.01979300812244555452", MPFR_RNDU, "1.95996734891603664741e-02");
102 check2("1.02900871924604464525", MPFR_RNDU, "2.85959303301472761438e-02");
106 check2("1.27832870030418943585", MPFR_RNDU, "2.45553521871417823608e-01");
110 check2("1.31706530746788241792", MPFR_RNDU, "2.75406009586277478185e-01");
114 check2("1.47116981099449883885", MPFR_RNDU, "3.86057874110010468272e-01");
118 check2("1.58405446812987782401", MPFR_RNDU, "4.59987679246663783150e-01");
122 check2("1.67192331263391547047", MPFR_RNDU, "5.13974647961076724911e-01");
126 check2("1.71101198068990645318", MPFR_RNDU, "5.37084997042120426691e-01");
130 check2("1.72634853551388700588", MPFR_RNDU, "5.46008504786553716670e-01")
[all...] |
| rnd_mode.c | 51 case MPFR_RNDU:
|
| tsqrt_ui.c | 50 check (2116118, MPFR_RNDU, "1.45468828276026215e3");
|
| tset_str.c | 86 mpfr_set_str (x, "1.00000000000000000006", 10, MPFR_RNDU);
240 mpfr_set_str (x, "9007199254740992", 10, MPFR_RNDU);
364 mpfr_set_str (y, "0.71128279983522479470@-160", 10, MPFR_RNDU);
382 mpfr_set_str (y, "0.jo08hg31hc5mmpj5mjjmgn55p2h35g@39", 27, MPFR_RNDU);
404 mpfr_set_str (y, "0.6e23ekb6acgh96abk10b6c9f2ka16i@45", 21, MPFR_RNDU);
435 mpfr_set_str (y, "0.jrchfhpp9en7hidqm9bmcofid9q3jg@39", 28, MPFR_RNDU);
630 CHECK53(y, "7e-303", MPFR_RNDU, x,
637 CHECK53(y, "93e-234", MPFR_RNDU, x,
644 CHECK53(y, "755e174", MPFR_RNDU, x,
651 CHECK53(y, "8699e-276", MPFR_RNDU, x
[all...] |
| tget_str.c | 58 s = mpfr_get_str (NULL, &e, 4, 2, x, MPFR_RNDU);
82 s = mpfr_get_str (NULL, &e, 3, 21, x, MPFR_RNDU);
90 s = mpfr_get_str (NULL, &e, 3, 20, x, MPFR_RNDU);
102 s = mpfr_get_str (NULL, &e, 3, 2, x, MPFR_RNDU);
162 s = mpfr_get_str (NULL, &e, 2, 2, x, MPFR_RNDU);
174 s = mpfr_get_str (NULL, &e, 10, 18, x, MPFR_RNDU);
180 s = mpfr_get_str (NULL, &e, 2, 4, x, MPFR_RNDU);
185 mpfr_out_str (stdout, 2, 4, x, MPFR_RNDU);
193 s = mpfr_get_str (NULL, &e, 8, 10, x, MPFR_RNDU);
203 s = mpfr_get_str (NULL, &e, 19, 12, x, MPFR_RNDU);
[all...] |
| texp.c | 86 test_exp(x, x, MPFR_RNDU);
152 mpfr_exp_2 (y, x, MPFR_RNDU);
153 mpfr_exp_3 (x, x, MPFR_RNDU);
392 test_exp (y, x, MPFR_RNDU);
550 if ((i >= 0 || emax < 0 || rnd == MPFR_RNDN || rnd == MPFR_RNDU) &&
624 mpfr_exp (y1, x, MPFR_RNDU);
628 mpfr_exp (y2, x, MPFR_RNDU);
629 mpfr_prec_round (y2, prec, MPFR_RNDU);
681 * but not necessarily in MPFR_RNDN and MPFR_RNDU (this is underflow
684 * - If eps > -2^(-p), no underflow in MPFR_RNDU
[all...] |
| tmul_ui.c | 120 mpfr_mul_ui (x, x, 3, MPFR_RNDU);
128 mpfr_mul_ui (x, x, 3, MPFR_RNDU);
136 mpfr_mul_ui (x, x, 3, MPFR_RNDU);
144 mpfr_mul_ui (x, x, 0, MPFR_RNDU);
147 mpfr_set_ui (x, 1, MPFR_RNDU);
148 mpfr_mul_ui (x, x, 0, MPFR_RNDU);
166 mpfr_mul_ui (x, x, 3, MPFR_RNDU);
|
| /src/external/lgpl3/mpc/dist/src/ |
| atan.c | 181 mpfr_sqr (z, y, MPFR_RNDU);
260 rnd1 = mpfr_sgn (mpc_realref (op)) > 0 ? MPFR_RNDD : MPFR_RNDU;
261 rnd2 = mpfr_sgn (mpc_realref (op)) < 0 ? MPFR_RNDU : MPFR_RNDD;
283 mpfr_atan2 (x, mpc_realref (op), a, MPFR_RNDU);
301 mpfr_sub (x, x, b, MPFR_RNDU);
307 mpfr_div_2ui (x, x, 1, MPFR_RNDU);
311 ok = mpfr_can_round (x, p - err, MPFR_RNDU, MPFR_RNDD,
348 mpfr_sqr (a, a, MPFR_RNDU);
349 mpfr_sqr (y, mpc_realref (op), MPFR_RNDU);
350 mpfr_add (a, a, y, MPFR_RNDU);
[all...] |
| sqrt.c | 42 const mpfr_rnd_t r = im_sgn ? MPFR_RNDD : MPFR_RNDU;
179 rnd_t = (im_cmp > 0 ? MPFR_RNDD : MPFR_RNDU);
181 rnd_t = (im_cmp > 0 ? MPFR_RNDU : MPFR_RNDD);
193 rnd_t = MPFR_RNDU;
199 rnd_t = MPFR_RNDU;
233 ok_w = mpfr_can_round (w, prec - 4, MPFR_RNDD, MPFR_RNDU,
300 if (rnd_w == MPFR_RNDU || rnd_w == MPFR_RNDA) {
308 if (rnd_w == MPFR_RNDU || rnd_w == MPFR_RNDA) {
316 if (rnd_w == MPFR_RNDU || rnd_w == MPFR_RNDA) {
341 and r = MPFR_RNDU
[all...] |
| /src/external/lgpl3/mpfr/dist/tune/ |
| bidimensional_sample.c | 130 mpfr_set_ui (temp, (unsigned long int)param.max_prec, MPFR_RNDU);
131 mpfr_div_ui (temp, temp, (unsigned long int)param.min_prec, MPFR_RNDU);
133 (unsigned long int)param.nb_points_prec, MPFR_RNDU);
134 incr_prec = mpfr_get_d (temp, MPFR_RNDU);
158 MPFR_RNDU);
164 MPFR_RNDU);
165 mpfr_exp2 (incr_x, incr_x, MPFR_RNDU);
171 MPFR_RNDU);
172 mpfr_exp2 (incr_x, incr_x, MPFR_RNDU);
183 mpfr_set_d (temp, param.min_x, MPFR_RNDU);
[all...] |
