HomeSort by: relevance | last modified time | path
    Searched refs:err_s (Results 1 - 3 of 3) sorted by relevancy

  /src/external/lgpl3/mpfr/dist/src/
jyn_asympt.c 45 mpfr_t s, c, P, Q, t, iz, err_t, err_s, err_u; local
81 mpfr_init2 (err_s, 31);
110 mpfr_set_ui (err_s, 0, MPFR_RNDN); /* error on P and Q (sum of errors) */
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);
136 mpfr_mul_2ui (err_s, err_s, w, MPFR_RNDU)
    [all...]
lngamma.c 176 mpfr_exp_t err_s, err_t; local
361 err_s = MPFR_GET_EXP(s); /* 2-x <= 2^err_s */
366 err_s += 3 - MPFR_GET_EXP(s);
367 err_s = (err_s >= 0) ? err_s + 1 : 0;
368 /* the error on s is bounded by 2^err_s ulp(s), thus by
369 2^(err_s+1)*2^(-w)*|s| since ulp(s) <= 2^(1-w)*|s|.
371 |r|<=2^(-w), thus taking n=2^(err_s+1) we see tha
    [all...]
gamma_inc.c 357 mpfr_exp_t exp_s, new_exp_s, exp_t, err_s, logj; local
425 err_s = exp_s - new_exp_s;
428 is bounded by 2^err_s ulp(s) */
439 /* s = exp(-x) * (S +/- 2^err_s ulp(S)) * (1 + theta)^2.
440 = exp(-x) * (S +/- 2^err_s ulp(S)) * (1 +/- 3 ulp(1))
442 exp(-x) * [2^err_s*ulp(S) + S*3*ulp(1) + 2^err_s*ulp(S)*3*ulp(1)]
443 <= ulp(s) * [2^(err_s+1) + 6 + 1]
444 <= ulp(s) * 2^(err_s+2) as long as err_s >= 2. *
    [all...]

Completed in 23 milliseconds