diff options
Diffstat (limited to 'src/math/__fpclassifyl.c')
| -rw-r--r-- | src/math/__fpclassifyl.c | 7 |
1 files changed, 5 insertions, 2 deletions
diff --git a/src/math/__fpclassifyl.c b/src/math/__fpclassifyl.c index a5ad36f2..e4d231b5 100644 --- a/src/math/__fpclassifyl.c +++ b/src/math/__fpclassifyl.c @@ -7,8 +7,11 @@ int __fpclassifyl(long double x) { union ldshape u = { x }; 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) return u.bits.m & (uint64_t)-1>>1 ? FP_NAN : FP_INFINITE; return u.bits.m & (uint64_t)1<<63 ? FP_NORMAL : FP_NAN; |