|
@@ -7,8 +7,11 @@ int __fpclassifyl(long double x)
|
|
{
|
|
{
|
|
union ldshape u = { x };
|
|
union ldshape u = { x };
|
|
int e = u.bits.exp;
|
|
int e = u.bits.exp;
|
|
- if (!e)
|
|
|
|
- return u.bits.m ? FP_SUBNORMAL : FP_ZERO;
|
|
|
|
|
|
+ if (!e) {
|
|
|
|
+ if (u.bits.m >> 63) return FP_NAN;
|
|
|
|
+ else if (u.bits.m) return FP_SUBNORMAL;
|
|
|
|
+ else return FP_ZERO;
|
|
|
|
+ }
|
|
if (e == 0x7fff)
|
|
if (e == 0x7fff)
|
|
return u.bits.m & (uint64_t)-1>>1 ? FP_NAN : FP_INFINITE;
|
|
return u.bits.m & (uint64_t)-1>>1 ? FP_NAN : FP_INFINITE;
|
|
return u.bits.m & (uint64_t)1<<63 ? FP_NORMAL : FP_NAN;
|
|
return u.bits.m & (uint64_t)1<<63 ? FP_NORMAL : FP_NAN;
|