improved expm1/log1p for small argument

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@3487 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 1 insertions, 1 deletions

diff --git a/expm1.c b/expm1.c
index 914afa3e0..dec6b90b5 100644
--- a/expm1.c
+++ b/expm1.c
@@ -73,7 +73,7 @@ mpfr_expm1 (mpfr_ptr y, mpfr_srcptr x , mp_rnd_t rnd_mode)
 
     /* compute the precision of intermediary variable */
     /* the optimal number of bits : see algorithms.tex */
-    Nt = Ny + MPFR_INT_CEIL_LOG2 (Nt) + 5;
+    Nt = Ny + MPFR_INT_CEIL_LOG2 (Ny) + 6;
 
     /* if |x| is smaller than 2^(-e), we will loose about e bits in the
        subtraction exp(x) - 1 */