#include <math.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double logbl(long double x)
{
        return logb(x);
}
#else
long double logbl(long double x)
{
        if (!isfinite(x))
                return x * x;
        if (x == 0)
                return -1/(x*x);
        return ilogbl(x);
}
#endif