summaryrefslogtreecommitdiff
path: root/src/math
diff options
context:
space:
mode:
authorSzabolcs Nagy <nsz@port70.net>2015-02-08 17:41:56 +0100
committerSzabolcs Nagy <nsz@port70.net>2015-02-08 17:41:56 +0100
commit3f92f92cb9c448ba82da5ff2cbbea692c7701c0c (patch)
tree80ec3415bf8ebaf5d70fc8a91aad354862ae299c /src/math
parent6e76e1540fc58a418494bf5eb832b556f9c5763e (diff)
downloadmusl-3f92f92cb9c448ba82da5ff2cbbea692c7701c0c.tar.gz
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')
-rw-r--r--src/math/__fpclassifyl.c5
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