OpenE2K
/
glibc
2
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

[BZ #16823] Fix log1pl returning wrong infinity sign

This commit is contained in:
Stefan Liebler 2014-04-29 15:43:36 +02:00 committed by Andreas Krebbel
parent dc041bd4db
commit 2ca180e97a
4 changed files with 11 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
ChangeLog
View File

@ -1,3 +1,11 @@
2014-04-29 Stefan Liebler <stli@linux.vnet.ibm.com>
[BZ #16823]
* sysdeps/ieee754/ldbl-128/s_log1pl.c (__log1pl):
Always divide by positive zero when computing -Inf result.
* sysdeps/ieee754/dbl-64/s_log1p.c (__log1p): Likewise.
* sysdeps/ieee754/flt-32/s_log1pf.c (__log1pf): Likewise.
2014-04-28 Adhemerval Zanella <azanella@linux.vnet.ibm.com>
* sysdeps/powerpc/fpu/fclrexcpt.c (__feclearexcept): Do not update

2
sysdeps/ieee754/dbl-64/s_log1p.c
View File

@ -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 */
}

2
sysdeps/ieee754/flt-32/s_log1pf.c
View File

@ -42,7 +42,7 @@ __log1pf(float x)
k = 1;
if (hx < 0x3ed413d7) { /* x < 0.41422 */
if(ax>=0x3f800000) { /* x <= -1.0 */
if(x==(float)-1.0) return -two25/(x-x); /* log1p(-1)=+inf */
if(x==(float)-1.0) return -two25/zero; /* log1p(-1)=-inf */
else return (x-x)/(x-x); /* log1p(x<-1)=NaN */
}
if(ax<0x31000000) { /* |x| < 2**-29 */

2
sysdeps/ieee754/ldbl-128/s_log1pl.c
View File

@ -150,7 +150,7 @@ __log1pl (long double xm1)
if (x <= 0.0L)
{
if (x == 0.0L)
return (-1.0L / (x - x));
return (-1.0L / zero); /* log1p(-1) = -inf */
else
return (zero / (x - x));
}