aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAndreas Jaeger <aj@suse.de>2003-12-28 18:57:35 +0000
committerAndreas Jaeger <aj@suse.de>2003-12-28 18:57:35 +0000
commit3737a0d70cac3ba09aa8f32709c4be8d70f41207 (patch)
tree59b4bff748ed586c59475201ba0662f6efd1f30e
parent6b6557e8b3b094132c619e3a2e00fe28422fd16f (diff)
downloadglibc-3737a0d70cac3ba09aa8f32709c4be8d70f41207.tar
glibc-3737a0d70cac3ba09aa8f32709c4be8d70f41207.tar.gz
glibc-3737a0d70cac3ba09aa8f32709c4be8d70f41207.tar.bz2
glibc-3737a0d70cac3ba09aa8f32709c4be8d70f41207.zip
Do not raise execptions for exp(NaN).
-rw-r--r--sysdeps/ieee754/dbl-64/e_exp.c4
1 files changed, 2 insertions, 2 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_exp.c b/sysdeps/ieee754/dbl-64/e_exp.c
index a6463e5874..717469e250 100644
--- a/sysdeps/ieee754/dbl-64/e_exp.c
+++ b/sysdeps/ieee754/dbl-64/e_exp.c
@@ -92,10 +92,10 @@ double __ieee754_exp(double x) {
if (n <= smallint) return 1.0;
if (n >= badint) {
- if (n > infint) return(zero/zero); /* x is NaN, return invalid */
+ if (n > infint) return(x+x); /* x is NaN */
if (n < infint) return ( (x>0) ? (hhuge*hhuge) : (tiny*tiny) );
/* x is finite, cause either overflow or underflow */
- if (junk1.i[LOW_HALF] != 0) return (zero/zero); /* x is NaN */
+ if (junk1.i[LOW_HALF] != 0) return (x+x); /* x is NaN */
return ((x>0)?inf.x:zero ); /* |x| = inf; return either inf or 0 */
}