diff options
| author | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2012-06-28 11:46:36 +0000 |
|---|---|---|
| committer | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2012-06-28 11:46:36 +0000 |
| commit | d52b9345e935e26971639dcda8970694140f7298 (patch) | |
| tree | 8b5679fa4b51eec0e0ebd9523c8c94dc99cd4de4 /tests | |
| parent | dccc01961014c5dfdb18a4556a8f2b65e1932263 (diff) | |
| download | mpc-d52b9345e935e26971639dcda8970694140f7298.tar.gz | |
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
Diffstat (limited to 'tests')
| -rw-r--r-- | tests/log.dat | 8 |
1 files changed, 4 insertions, 4 deletions
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) |