X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flogbl.c;h=962973a7b99bfffdcb2c1677a3c0b2ba5f87d357;hb=246f1c811448f37a44b41cd8df8d0ef9736d95f4;hp=5d04abd764404dbb166ae01d18c25cba9b07d679;hpb=b69f695acedd4ce2798ef9ea28d834ceccc789bd;p=musl diff --git a/src/math/logbl.c b/src/math/logbl.c index 5d04abd7..962973a7 100644 --- a/src/math/logbl.c +++ b/src/math/logbl.c @@ -1,5 +1,4 @@ -#include -#include "libm.h" +#include #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*x); + return ilogbl(x); } #endif