diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2016-12-12 17:32:09 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2016-12-12 17:32:09 +0000 |
commit | ed8cf60e5755ada028226fd32f4296728a73b524 (patch) | |
tree | 851abce87d2c3aec86e8ca4880df80d9b8aa291f /src/cos.c | |
parent | 5b5002d48d93a08993accb6d762c3795692ceb04 (diff) | |
download | mpfr-ed8cf60e5755ada028226fd32f4296728a73b524.tar.gz |
[cos.c,exp_2.c,log.c] improved initial precision for Ziv's strategy
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@11031 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src/cos.c')
-rw-r--r-- | src/cos.c | 4 |
1 files changed, 2 insertions, 2 deletions
@@ -174,7 +174,7 @@ mpfr_cos (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) } K0 = __gmpfr_isqrt (precy / 3); - m = precy + 2 * MPFR_INT_CEIL_LOG2 (precy) + 2 * K0; + m = precy + 2 * MPFR_INT_CEIL_LOG2 (precy) + 2 * K0 + 4; if (expx >= 3) { @@ -245,7 +245,7 @@ mpfr_cos (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) l = 2 * l + 1; if (reduce) l += (K == 0) ? 4 : 1; - k = MPFR_INT_CEIL_LOG2 (l) + 2*K; + k = MPFR_INT_CEIL_LOG2 (l) + 2 * K; /* now the error is bounded by 2^(k-m) = 2^(EXP(s)-err) */ exps = MPFR_GET_EXP (s); |