diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2007-02-19 03:52:30 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2007-02-19 03:52:30 +0000 |
commit | 177eda24dfc9c95af80aeb8bfb4b2d51a891da45 (patch) | |
tree | 91fd1500cc012795930a2473c2b589756c090dc6 /log1p.c | |
parent | 4c9e9473ec7f62cbfda1a6f72db816ad0eb67627 (diff) | |
download | mpfr-177eda24dfc9c95af80aeb8bfb4b2d51a891da45.tar.gz |
fixed error in generic error for the logarithm, and propagated new bound
in algorithms.tex and source files
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@4372 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'log1p.c')
-rw-r--r-- | log1p.c | 4 |
1 files changed, 2 insertions, 2 deletions
@@ -123,8 +123,8 @@ mpfr_log1p (mpfr_ptr y, mpfr_srcptr x, mp_rnd_t rnd_mode) mpfr_log (t, t, GMP_RNDN); /* log(1+x)*/ /* estimation of the error */ - /*err=Nt-(__gmpfr_ceil_log2(1+pow(2,1-MPFR_GET_EXP(t))));*/ - err = Nt - (MAX (1 - MPFR_GET_EXP (t), 0) + 1); + /*err=Nt-(__gmpfr_ceil_log2(1+pow(2,2-MPFR_GET_EXP(t))));*/ + err = Nt - (MAX (2 - MPFR_GET_EXP (t), 0) + 1); if (MPFR_LIKELY (MPFR_CAN_ROUND (t, err, Ny, rnd_mode))) break; |