summaryrefslogtreecommitdiff
path: root/src/log.c
diff options
context:
space:
mode:
authorvlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2016-12-27 02:33:59 +0000
committervlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2016-12-27 02:33:59 +0000
commit31bbca629e7d379e51f9ddfc561c9c53483ad70d (patch)
treef86d9f6d89257c3762a7105f8c0becd9e9243709 /src/log.c
parentdfcab06a6496da16b93c52131f0c88cd7b06ba9d (diff)
downloadmpfr-31bbca629e7d379e51f9ddfc561c9c53483ad70d.tar.gz
[src/log.c] Added another FIXME concerning the error analysis.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@11092 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src/log.c')
-rw-r--r--src/log.c3
1 files changed, 3 insertions, 0 deletions
diff --git a/src/log.c b/src/log.c
index 0fbe58622..4b428070f 100644
--- a/src/log.c
+++ b/src/log.c
@@ -134,6 +134,9 @@ mpfr_log (mpfr_ptr r, mpfr_srcptr a, mpfr_rnd_t rnd_mode)
true. */
MPFR_ASSERTN (m >= LONG_MIN && m <= LONG_MAX);
+ /* FIXME: Why 1 ulp and not 1/2 ulp? Ditto with some other ones
+ below. The error concerning the AGM should be explained since
+ 4/s is inexact (one needs a bound on its derivative). */
mpfr_mul_2si (tmp2, a, m, MPFR_RNDN); /* s=a*2^m, err<=1 ulp */
MPFR_ASSERTD (MPFR_EXP (tmp2) >= (p + 3) / 2);
/* [FIXME] and one can have the equality, even if p is even.