From ed8cf60e5755ada028226fd32f4296728a73b524 Mon Sep 17 00:00:00 2001 From: zimmerma Date: Mon, 12 Dec 2016 17:32:09 +0000 Subject: [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 --- src/cos.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src/cos.c') diff --git a/src/cos.c b/src/cos.c index 36d2b0dfc..95adf4f83 100644 --- a/src/cos.c +++ b/src/cos.c @@ -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); -- cgit v1.2.1