diff options
author | Szabolcs Nagy <nsz@port70.net> | 2015-02-08 17:41:56 +0100 |
---|---|---|
committer | Szabolcs Nagy <nsz@port70.net> | 2015-02-08 17:41:56 +0100 |
commit | 3f92f92cb9c448ba82da5ff2cbbea692c7701c0c (patch) | |
tree | 80ec3415bf8ebaf5d70fc8a91aad354862ae299c | |
parent | 6e76e1540fc58a418494bf5eb832b556f9c5763e (diff) | |
download | musl-3f92f92cb9c448ba82da5ff2cbbea692c7701c0c.tar.gz musl-3f92f92cb9c448ba82da5ff2cbbea692c7701c0c.tar.bz2 musl-3f92f92cb9c448ba82da5ff2cbbea692c7701c0c.tar.xz musl-3f92f92cb9c448ba82da5ff2cbbea692c7701c0c.zip |
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.
-rw-r--r-- | src/math/__fpclassifyl.c | 5 |
1 files 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 |