Searched refs:precy (Results 1 - 13 of 13) sorted by relevance
/netbsd-6-1-5-RELEASE/external/lgpl3/mpfr/dist/ |
H A D | 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)))
|
H A D | tan.c | 30 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 D | sin.c | 42 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 D | cos.c | 133 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 D | get_f.c | 34 mpfr_prec_t precx, precy; local 80 precy = MPFR_PREC(y); 91 if (precy + sh <= precx) /* we can copy directly */
|
H A D | exp.c | 31 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 D | exp_2.c | 84 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 D | coth.c | 46 if (MPFR_IS_ZERO(z) || MPFR_GET_EXP(z) <= - (mpfr_exp_t) precy) \
|
H A D | lngamma.c | 98 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 D | exp3.c | 49 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 D | atan.c | 45 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 D | texp.c | 357 /* 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 D | tcomparisons.c | 39 mpfr_prec_t precx, precy; local 44 precy = (randlimb () % 17) * 11 + MPFR_PREC_MIN; 46 mpfr_set_prec (y, precy);
|
Completed in 74 milliseconds