diff options
author | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2011-09-02 12:41:35 +0000 |
---|---|---|
committer | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2011-09-02 12:41:35 +0000 |
commit | e5ec494194c5585e2736354fd0bf284db0ba729b (patch) | |
tree | ee86d00fbf2b5e263d055c37a31eb326902e315d /BUGS | |
parent | 26d2d8b13d1b199493f45e63b27aa084f4cfa09c (diff) | |
download | mpc-e5ec494194c5585e2736354fd0bf284db0ba729b.tar.gz |
norm.c: rewritten to handle all cases of over- and underflow and to obtain a
correct ternary return value
norm.dat: new tests with over- and underflow
mpc-impl.h: added macro MPC_MIN
BUGS, TODO, mpc.texi updated to reflect these changes
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@1071 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'BUGS')
-rw-r--r-- | BUGS | 9 |
1 files changed, 6 insertions, 3 deletions
@@ -1,8 +1,7 @@ - Overflow and underflow are not considered in all functions, which might lead to some infinite loops. -- Except for mpc_mul, - no checks are made for intermediate over- or underflows, which may occur - in the middle of the algorithms although the final result may be +- Also, checks are not always made for intermediate over- or underflows, which + may occur in the middle of the algorithms although the final result may be representable. For instance, if in the computation of Im (cos(x+i*y)) = -sin(x)*sinh(y) an overflow occurs in sinh(y), the value sign(Im (cos(x+i*y))) * inf @@ -10,3 +9,7 @@ sin(x)*sinh(y) is representable. If furthermore an underflow occurred in sin(x) (which has not been observed in practice), then the return value would be NaN*(+-inf)=NaN. + +Currently, the following functions have been rewritten to solve these +problems: +mul, norm |