jyapayne/musl
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

math: fix logb*.c exceptions now that ilogb raises invalid

Browse source
This commit is contained in:
Szabolcs Nagy 2012-11-13 00:45:51 +01:00
commit 52bcf3f5fb
3 changed files with 17 additions and 25 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/math/logbl.c
View file

@ -1,4 +1,3 @@
#include <limits.h>
#include "libm.h"
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double logbl(long double x)
@ -8,12 +7,10 @@ long double logbl(long double x)
#else
long double logbl(long double x)
{
int i = ilogbl(x);
if (i == FP_ILOGB0)
return -1.0/fabsl(x);
if (i == FP_ILOGBNAN || i == INT_MAX)
if (!isfinite(x))
return x * x;
return i;
if (x == 0)
return -1/(x+0);
return ilogbl(x);
}
#endif