From 22ec468b9cec80aa3ffdf58dc853b8cf70f5236c Mon Sep 17 00:00:00 2001 From: enge Date: Sat, 30 Jun 2012 13:51:19 +0000 Subject: 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 --- src/log.c | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) (limited to 'src') 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 */ -- cgit v1.2.1