From d52b9345e935e26971639dcda8970694140f7298 Mon Sep 17 00:00:00 2001 From: enge Date: Thu, 28 Jun 2012 11:46:36 +0000 Subject: log: alternative implementation that avoids intermediate overflows It is probably slower (two calls to log) and should be combined with the previous approach. Problem of "underflow" (log of number close to 1) not yet solved. git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@1202 211d60ee-9f03-0410-a15a-8952a2c7a4e4 --- tests/log.dat | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) (limited to 'tests') diff --git a/tests/log.dat b/tests/log.dat index 0f590d0..6972cfc 100644 --- a/tests/log.dat +++ b/tests/log.dat @@ -166,10 +166,10 @@ + + 53 0x162E4301C0025p-23 2 0x1p0 2 0x1p67108864 2 0x1p67108864 U U + + 53 0xB1721802E8F1p-21 2 0x1p0 2 0x1p134217728 2 0x1p134217728 U U + + 53 0x2C5C85FF57581p-22 2 0x1p0 2 0x1p268435456 2 0x1p268435456 U U -# Due to intermediate overflow, the following result has real part infinity -# instead of the correct result. Since this may happen in other parts of the -# library as well, we do not consider it a bug for the moment. -# + + 53 0x58B90BFD4BCBFp-22 2 0x1p0 2 0x1p536870912 2 0x1p536870912 U U + +# Example leading to intermediate overflow in x^2+y^2 ++ + 53 0x58B90BFD4BCBFp-22 2 0x1p0 2 0x1p536870912 2 0x1p536870912 U U + # Due to intermediate "underflow" (see FIXME in log.c), the following example # loops. # log(-1 + i*eps) -- cgit v1.2.1