math: fix __fpclassifyl(-0.0) for IEEE binary128

The sign bit was not cleared before checking for 0 so -0.0
was misclassified as FP_SUBNORMAL instead of FP_ZERO.

diff --git a/src/math/__fpclassifyl.c b/src/math/__fpclassifyl.c
index c274229..481c0b9 100644 (file)
--- a/src/math/__fpclassifyl.c
+++ b/src/math/__fpclassifyl.c
@@ -24,12 +24,11 @@ int __fpclassifyl(long double x)
 {
        union ldshape u = {x};
        int e = u.i.se & 0x7fff;
 {
        union ldshape u = {x};
        int e = u.i.se & 0x7fff;

+       u.i.se = 0;

        if (!e)
                return u.i2.lo | u.i2.hi ? FP_SUBNORMAL : FP_ZERO;
        if (!e)
                return u.i2.lo | u.i2.hi ? FP_SUBNORMAL : FP_ZERO;

-       if (e == 0x7fff) {
-               u.i.se = 0;
+       if (e == 0x7fff)

                return u.i2.lo | u.i2.hi ? FP_NAN : FP_INFINITE;
                return u.i2.lo | u.i2.hi ? FP_NAN : FP_INFINITE;

-       }

        return FP_NORMAL;
 }
 #endif
        return FP_NORMAL;
 }
 #endif