summaryrefslogtreecommitdiff
path: root/src/log.c
diff options
context:
space:
mode:
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.