diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2016-12-27 02:33:59 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2016-12-27 02:33:59 +0000 |
commit | 31bbca629e7d379e51f9ddfc561c9c53483ad70d (patch) | |
tree | f86d9f6d89257c3762a7105f8c0becd9e9243709 /src/log.c | |
parent | dfcab06a6496da16b93c52131f0c88cd7b06ba9d (diff) | |
download | mpfr-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.c | 3 |
1 files changed, 3 insertions, 0 deletions
@@ -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. |