diff options
author | Stefan Liebler <stli@linux.vnet.ibm.com> | 2014-04-29 15:43:36 +0200 |
---|---|---|
committer | Andreas Krebbel <krebbel@linux.vnet.ibm.com> | 2014-04-29 15:43:36 +0200 |
commit | 2ca180e97a16f3ee3ec9e7e5a5fb7eda3d863a02 (patch) | |
tree | 4b1d6f28281ad75b3abe35ee09684e1aae943012 /sysdeps/ieee754/dbl-64/s_log1p.c | |
parent | dc041bd4dbaf0f22b962647e0d0a1b8090d0e679 (diff) | |
download | glibc-2ca180e97a16f3ee3ec9e7e5a5fb7eda3d863a02.zip glibc-2ca180e97a16f3ee3ec9e7e5a5fb7eda3d863a02.tar.gz glibc-2ca180e97a16f3ee3ec9e7e5a5fb7eda3d863a02.tar.bz2 |
[BZ #16823] Fix log1pl returning wrong infinity sign
Diffstat (limited to 'sysdeps/ieee754/dbl-64/s_log1p.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/s_log1p.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_log1p.c b/sysdeps/ieee754/dbl-64/s_log1p.c index fd4dce5..c922148 100644 --- a/sysdeps/ieee754/dbl-64/s_log1p.c +++ b/sysdeps/ieee754/dbl-64/s_log1p.c @@ -110,7 +110,7 @@ __log1p (double x) if (__glibc_unlikely (ax >= 0x3ff00000)) /* x <= -1.0 */ { if (x == -1.0) - return -two54 / (x - x); /* log1p(-1)=+inf */ + return -two54 / zero; /* log1p(-1)=-inf */ else return (x - x) / (x - x); /* log1p(x<-1)=NaN */ } |