summaryrefslogtreecommitdiff
path: root/src/cos.c
diff options
context:
space:
mode:
authorzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2016-12-12 17:32:09 +0000
committerzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2016-12-12 17:32:09 +0000
commited8cf60e5755ada028226fd32f4296728a73b524 (patch)
tree851abce87d2c3aec86e8ca4880df80d9b8aa291f /src/cos.c
parent5b5002d48d93a08993accb6d762c3795692ceb04 (diff)
downloadmpfr-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.c4
1 files changed, 2 insertions, 2 deletions
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);