| /src/external/lgpl3/mpfr/dist/src/ |
| tan.c | 30 mpfr_prec_t precy, m; local
65 precy = MPFR_PREC (y);
66 m = precy + MPFR_INT_CEIL_LOG2 (precy) + 13;
78 if (MPFR_LIKELY (MPFR_CAN_ROUND (c, m - 2, precy, rnd_mode)))
|
| gen_inverse.h | 46 mpfr_prec_t precy; /* target precision */ local
66 precy = MPFR_PREC(y);
67 m = precy + MPFR_INT_CEIL_LOG2 (precy) + 3;
92 if (MPFR_LIKELY (MPFR_CAN_ROUND (z, m - 2, precy, rnd_mode)))
|
| exp.c | 39 mpfr_prec_t precy; local
132 precy = MPFR_PREC (y);
134 /* if x < 2^(-precy), then exp(x) gives 1 +/- 1 ulp(1) */
135 if (MPFR_UNLIKELY (expx < 0 && (mpfr_uexp_t) (-expx) > precy))
173 if (MPFR_UNLIKELY (precy >= MPFR_EXP_THRESHOLD))
|
| get_f.c | 36 mpfr_prec_t precx, precy; local
82 precy = MPFR_PREC(y);
93 if (precy + sh <= precx) /* we can copy directly */
|
| sin.c | 42 mpfr_prec_t precy, m; local
77 precy = MPFR_PREC (y);
79 if (precy >= MPFR_SINCOS_THRESHOLD)
87 m = precy + MPFR_INT_CEIL_LOG2 (MAX(precy,expx)) + 8;
95 At this point, precy < MPFR_SINCOS_THRESHOLD, so that both m and
182 if (MPFR_CAN_ROUND (c, err, precy, rnd_mode))
|
| fma.c | 107 mpfr_prec_t precx, precy; local
127 precy = MPFR_PREC (y);
135 if (precx == precy && e <= __gmpfr_emax && e > __gmpfr_emin)
217 MPFR_ASSERTN (precx + precy <= MPFR_PREC_MAX);
218 MPFR_GROUP_INIT_1 (group, precx + precy, u);
|
| sinu.c | 38 mpfr_prec_t precy, prec; local
107 precy = MPFR_GET_PREC (y);
111 prec = precy + MAX(expx, MPFR_INT_CEIL_LOG2 (precy)) + 8;
153 if (MPFR_CAN_ROUND (t, err, precy, rnd_mode))
|
| tanu.c | 32 mpfr_prec_t precy, prec; local
102 precy = MPFR_GET_PREC (y);
106 prec = precy + MAX(expx,MPFR_INT_CEIL_LOG2(precy)) + 8;
168 if (MPFR_CAN_ROUND (t, err, precy, rnd_mode))
|
| cos.c | 135 mpfr_prec_t K0, K, precy, m, k, l; local
170 precy = MPFR_PREC (y);
172 if (precy >= MPFR_SINCOS_THRESHOLD)
178 K0 = __gmpfr_isqrt (precy / 3);
179 m = precy + 2 * MPFR_INT_CEIL_LOG2 (precy) + 2 * K0 + 4;
254 if (MPFR_LIKELY (MPFR_CAN_ROUND (s, exps + m - k, precy, rnd_mode)))
260 than 2^(-precy) for directed rounding, or 2^(-precy-1) for rounding
263 if (m > k && (m - k >= precy + (rnd_mode == MPFR_RNDN))
[all...] |
| cosu.c | 32 mpfr_prec_t precy, prec; local
126 precy = MPFR_GET_PREC (y);
129 prec = precy + MAX(expx,MPFR_INT_CEIL_LOG2 (precy)) + 8;
172 if (MPFR_CAN_ROUND (t, err, precy, rnd_mode))
|
| exp3.c | 49 mpfr_prec_t precy = MPFR_PREC(y), prec_i_have, prec_ptoj; local
80 for (i = 1; (prec_i_have < precy) && (i < n); i++)
132 diff = (mpfr_exp_t) prec_i_have - 2 * (mpfr_exp_t) precy;
140 diff = (mpfr_exp_t) prec_i_have - (mpfr_prec_t) precy;
|
| exp_2.c | 82 mpfr_prec_t q, precy; local
95 precy = MPFR_PREC(y);
156 K = (precy < MPFR_EXP_2_THRESHOLD) ? __gmpfr_isqrt ((precy + 1) / 2) + 3
157 : __gmpfr_cuberoot (4*precy);
158 l = (precy - 1) / K + 1;
160 /* add K extra bits, i.e. failure probability <= 1/2^K = O(1/precy) */
161 q = precy + err + K + 10;
228 l = (precy < MPFR_EXP_2_THRESHOLD)
251 if (MPFR_LIKELY (MPFR_CAN_ROUND (s, q - err, precy, rnd_mode))
[all...] |
| atan.c | 100 mpfr_prec_t precy = MPFR_PREC(y); local
109 if (precy <= 192)
188 p^i/2^(r*i) < 2^(-precy), i.e. r*i > precy + log2(p^i) */
214 if (accu[k-1] > precy)
220 we can stop when r*i > precy i.e. i > precy/r */
223 if (precy / r <= n)
224 n = (precy / r) + 1;
263 diff -= 2 * precy;
[all...] |
| lngamma.c | 172 mpfr_prec_t precy, w; /* working precision */ local
312 precy = MPFR_PREC(y);
330 w = precy + MPFR_INT_CEIL_LOG2 (precy);
418 if (MPFR_CAN_ROUND (s, w - err_s, precy, rnd))
432 w = precy + MPFR_INT_CEIL_LOG2 (precy);
691 if (MPFR_LIKELY (MPFR_CAN_ROUND (s, w - err_s, precy, rnd)))
|
| /src/external/lgpl3/mpfr/dist/tests/ |
| tcomparisons.c | 36 mpfr_prec_t precx, precy; local
41 precy = (randlimb () % 17) * 11 + MPFR_PREC_MIN;
43 mpfr_set_prec (y, precy);
|
| tdigamma.c | 63 int i, precx, precy[NPREC] = { 200, 400, 520, 1416 }; local
71 mpfr_init2 (y[i], precy[i]);
72 mpfr_init2 (z, precy[0]);
|
| tcan_round.c | 232 int precx, precy, err; local
244 for (precy = precx - 4; precy <= precx + 4; precy++)
246 mpfr_set_prec (yinf, precy);
247 mpfr_set_prec (ysup, precy);
265 maxerr = precy + 3;
268 for (err = precy; err <= maxerr; err++)
292 (mpfr_rnd_t) rnd2, precy);
296 "precx=%d, precy=%d, i=%d, err=%d,
[all...] |
| tcompound.c | 403 mpfr_prec_t precy, prect; local
430 precy = MPFR_GET_PREC (y);
431 prect = precy + 20;
468 if (MPFR_CAN_ROUND (t, prect + 2 - err, precy, rnd_mode))
|
| texp.c | 366 /* check case EXP(x) < -precy */
659 int precx, precy; local
691 * Moreover, since precy < MPFR_EXP_THRESHOLD (to avoid tests that take
701 * for precx = 96, precy = 16, mpfr_exp_3
711 for (precy = 16; precy <= 128; precy += 16)
713 mpfr_init2 (y, precy);
731 printf ("\nfor precx = %d, precy = %d, %s\n",
732 precx, precy, e3 ? "mpfr_exp_3" : "mpfr_exp")
[all...] |
