diff options
| -rw-r--r-- | sys/src/libc/port/frexp.c | 13 |
1 files changed, 11 insertions, 2 deletions
diff --git a/sys/src/libc/port/frexp.c b/sys/src/libc/port/frexp.c index 669f12152..e2f4fb072 100644 --- a/sys/src/libc/port/frexp.c +++ b/sys/src/libc/port/frexp.c @@ -89,6 +89,16 @@ modf(double d, double *ip) FPdbleword x; int e; + x.x = d; + e = (x.hi >> SHIFT) & MASK; + if(e == MASK){ + *ip = d; + if(x.lo != 0 || (x.hi & 0xfffffL != 0)) /* NaN */ + return d; + /* ±Inf */ + x.hi &= 0x80000000L; + return x.x; + } if(d < 1) { if(d < 0) { x.x = modf(-d, ip); @@ -98,8 +108,7 @@ modf(double d, double *ip) *ip = 0; return d; } - x.x = d; - e = ((x.hi >> SHIFT) & MASK) - BIAS; + e -= BIAS; if(e <= SHIFT+1) { x.hi &= ~(0x1fffffL >> e); x.lo = 0; |