From 3f92f92cb9c448ba82da5ff2cbbea692c7701c0c Mon Sep 17 00:00:00 2001 From: Szabolcs Nagy Date: Sun, 8 Feb 2015 17:41:56 +0100 Subject: [PATCH] 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. --- src/math/__fpclassifyl.c | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/math/__fpclassifyl.c b/src/math/__fpclassifyl.c index c2742297..481c0b94 100644 --- 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; + u.i.se = 0; 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 FP_NORMAL; } #endif