Searched refs:precy (Results 1 - 13 of 13) sorted by relevance

/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/
H A Dgen_inverse.h46 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)))
H A Dtan.c30 mpfr_prec_t precy, m; local
63 precy = MPFR_PREC (y);
64 m = precy + MPFR_INT_CEIL_LOG2 (precy) + 13;
76 if (MPFR_LIKELY (MPFR_CAN_ROUND (c, m - 2, precy, rnd_mode)))
H A Dsin.c42 mpfr_prec_t precy, m; local
74 precy = MPFR_PREC (y);
76 if (precy >= MPFR_SINCOS_THRESHOLD)
79 m = precy + MPFR_INT_CEIL_LOG2 (precy) + 13;
154 if (MPFR_CAN_ROUND (c, err, precy, rnd_mode))
H A Dcos.c133 mpfr_prec_t K0, K, precy, m, k, l; local
166 precy = MPFR_PREC (y);
168 if (precy >= MPFR_SINCOS_THRESHOLD)
174 K0 = __gmpfr_isqrt (precy / 3);
175 m = precy + 2 * MPFR_INT_CEIL_LOG2 (precy) + 2 * K0;
250 if (MPFR_LIKELY (MPFR_CAN_ROUND (s, exps + m - k, precy, rnd_mode)))
256 than 2^(-precy) for directed rounding, or 2^(-precy-1) for rounding
259 if (m > k && (m - k >= precy
[all...]
H A Dget_f.c34 mpfr_prec_t precx, precy; local
80 precy = MPFR_PREC(y);
91 if (precy + sh <= precx) /* we can copy directly */
H A Dexp.c31 mpfr_prec_t precy; local
105 precy = MPFR_PREC (y);
107 /* if x < 2^(-precy), then exp(x) i.e. gives 1 +/- 1 ulp(1) */
108 if (MPFR_UNLIKELY (expx < 0 && (mpfr_uexp_t) (-expx) > precy))
148 if (MPFR_UNLIKELY (precy >= MPFR_EXP_THRESHOLD))
H A Dexp_2.c84 mpfr_prec_t q, precy; local
93 precy = MPFR_PREC(y);
122 K = (precy < MPFR_EXP_2_THRESHOLD) ? __gmpfr_isqrt ((precy + 1) / 2)
123 : __gmpfr_cuberoot (4*precy);
124 l = (precy - 1) / K + 1;
126 /* add K extra bits, i.e. failure probability <= 1/2^K = O(1/precy) */
127 q = precy + err + K + 5;
190 l = (precy < MPFR_EXP_2_THRESHOLD)
216 if (MPFR_LIKELY (MPFR_CAN_ROUND (s, q - K, precy, rnd_mod
[all...]
H A Dcoth.c46 if (MPFR_IS_ZERO(z) || MPFR_GET_EXP(z) <= - (mpfr_exp_t) precy) \
H A Dlngamma.c98 mpfr_prec_t precy, w; /* working precision */ local
191 precy = MPFR_PREC(y);
207 w = precy + MPFR_INT_CEIL_LOG2 (precy);
287 if (mpfr_can_round (s, w - err_s, MPFR_RNDN, MPFR_RNDZ, precy
300 w = precy + MPFR_INT_CEIL_LOG2 (precy);
484 while (MPFR_UNLIKELY (!MPFR_CAN_ROUND (s, w - err_s, precy, rnd)));
H A Dexp3.c49 mpfr_prec_t precy = MPFR_PREC(y), prec_i_have, prec_ptoj; local
78 for (i = 1; (prec_i_have < precy) && (i < n); i++)
130 diff = (mpfr_exp_t) prec_i_have - 2 * (mpfr_exp_t) precy;
138 diff = (mpfr_exp_t) prec_i_have - (mpfr_prec_t) precy;
H A Datan.c45 mpfr_prec_t precy = MPFR_PREC(y); local
94 p^i/2^(r*i) < 2^(-precy), i.e. r*i > precy + log2(p^i) */
121 if (accu[k-1] > precy)
127 we can stop when r*i > precy i.e. i > precy/r */
130 for (i = k = 0; (i < n) && (i <= precy / r); i += 2, k ++)
169 diff -= 2 * precy;
177 diff -= precy;
/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/tests/
H A Dtexp.c357 /* check case EXP(x) < -precy */
646 int precx, precy; local
678 * Moreover, since precy < MPFR_EXP_THRESHOLD (to avoid tests that take
688 * for precx = 96, precy = 16, mpfr_exp_3
698 for (precy = 16; precy <= 128; precy += 16)
700 mpfr_init2 (y, precy);
717 printf ("\nfor precx = %d, precy = %d, %s\n",
718 precx, precy, e
[all...]
H A Dtcomparisons.c39 mpfr_prec_t precx, precy; local
44 precy = (randlimb () % 17) * 11 + MPFR_PREC_MIN;
46 mpfr_set_prec (y, precy);

Completed in 74 milliseconds