diff options
| author | Rich Felker <dalias@aerifal.cx> | 2012-11-15 08:37:58 -0500 |
|---|---|---|
| committer | Rich Felker <dalias@aerifal.cx> | 2012-11-15 08:37:58 -0500 |
| commit | 8bb181622222f2ee3462c8b021bcae4fcdbbd37a (patch) | |
| tree | 166a1be834c8cd5a4c2521d41d7ee506bd8c8e51 /src/math/ilogbl.c | |
| parent | update copyright file for recent contributions (diff) | |
| parent | fenv: return FE_TONEAREST in dummy fegetround (diff) | |
| download | grovel-8bb181622222f2ee3462c8b021bcae4fcdbbd37a.tar.gz grovel-8bb181622222f2ee3462c8b021bcae4fcdbbd37a.tar.xz | |
Merge remote-tracking branch 'nsz/math'
Diffstat (limited to 'src/math/ilogbl.c')
| -rw-r--r-- | src/math/ilogbl.c | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/src/math/ilogbl.c b/src/math/ilogbl.c index ed9ddcbc..1512934f 100644 --- a/src/math/ilogbl.c +++ b/src/math/ilogbl.c @@ -14,15 +14,19 @@ int ilogbl(long double x) int e = u.bits.exp; if (!e) { - if (m == 0) + if (m == 0) { + FORCE_EVAL(0/0.0f); return FP_ILOGB0; + } /* subnormal x */ for (e = -0x3fff+1; m < (uint64_t)1<<63; e--, m<<=1); return e; } - if (e == 0x7fff) + if (e == 0x7fff) { + FORCE_EVAL(0/0.0f); /* in ld80 msb is set in inf */ return m & (uint64_t)-1>>1 ? FP_ILOGBNAN : INT_MAX; + } return e - 0x3fff; } #endif |