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 /src/math/__fpclassifyl.c | |
| parent | 6e76e1540fc58a418494bf5eb832b556f9c5763e (diff) | |
| download | musl-3f92f92cb9c448ba82da5ff2cbbea692c7701c0c.tar.gz 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.
Diffstat (limited to 'src/math/__fpclassifyl.c')
| -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 |