-#include <limits.h>
-#include "libm.h"
+#include <math.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
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