summaryrefslogtreecommitdiff
path: root/src/math/ilogb.c
diff options
context:
space:
mode:
authorSzabolcs Nagy <nsz@port70.net>2012-11-12 23:58:18 +0100
committerSzabolcs Nagy <nsz@port70.net>2012-11-12 23:58:18 +0100
commitbe81f51030d90de8d5eb0b65dbdab3032d5cae23 (patch)
treedcbee14143ff6e6d580cad060d54194d1db5480b /src/math/ilogb.c
parent6f64505df37c41d7496612dffe86048765339a7a (diff)
downloadmusl-be81f51030d90de8d5eb0b65dbdab3032d5cae23.tar.gz
math: raise invalid flag in ilogb*.c on +-0, +-inf and nan
Diffstat (limited to 'src/math/ilogb.c')
-rw-r--r--src/math/ilogb.c8
1 files changed, 6 insertions, 2 deletions
diff --git a/src/math/ilogb.c b/src/math/ilogb.c
index 0a3a6a46..5a1819d8 100644
--- a/src/math/ilogb.c
+++ b/src/math/ilogb.c
@@ -8,13 +8,17 @@ int ilogb(double x)
if (!e) {
u.bits <<= 12;
- if (u.bits == 0)
+ if (u.bits == 0) {
+ FORCE_EVAL(0/0.0f);
return FP_ILOGB0;
+ }
/* subnormal x */
for (e = -0x3ff; u.bits < (uint64_t)1<<63; e--, u.bits<<=1);
return e;
}
- if (e == 0x7ff)
+ if (e == 0x7ff) {
+ FORCE_EVAL(0/0.0f);
return u.bits<<12 ? FP_ILOGBNAN : INT_MAX;
+ }
return e - 0x3ff;
}