log.dat: add example with overflow also with abs instead of norm

log.c: handle overflow again


git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@1212 211d60ee-9f03-0410-a15a-8952a2c7a4e4

1 files changed, 4 insertions, 1 deletions

diff --git a/src/log.c b/src/log.c
index 0af4b4e..70de44f 100644
--- a/src/log.c
+++ b/src/log.c
@@ -143,7 +143,10 @@ do {                                                            \
 
       mpc_abs (w, op, GMP_RNDN);
          /* error 0.5 ulp */
-      MPC_ASSERT (!mpfr_inf_p (w) && !mpfr_zero_p (w));
+      if (mpfr_inf_p (w))
+         /* intermediate overflow; the logarithm may be representable.
+            Intermediate underflow is impossible.                      */
+         break;
 
       mpfr_log (w, w, GMP_RNDN);
          /* generic error of log: (2^(- exp(w)) + 0.5) ulp */