diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2012-11-12 23:58:18 +0100 |
|---|---|---|
| committer | Szabolcs Nagy <nsz@port70.net> | 2012-11-12 23:58:18 +0100 |
| commit | be81f51030d90de8d5eb0b65dbdab3032d5cae23 (patch) | |
| tree | dcbee14143ff6e6d580cad060d54194d1db5480b /src/math/ilogb.c | |
| parent | 6f64505df37c41d7496612dffe86048765339a7a (diff) | |
| download | musl-be81f51030d90de8d5eb0b65dbdab3032d5cae23.tar.gz musl-be81f51030d90de8d5eb0b65dbdab3032d5cae23.tar.bz2 musl-be81f51030d90de8d5eb0b65dbdab3032d5cae23.tar.xz musl-be81f51030d90de8d5eb0b65dbdab3032d5cae23.zip | |
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.c | 8 |
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; } |