diff options
author | Szabolcs Nagy <nsz@port70.net> | 2013-09-03 04:09:12 +0000 |
---|---|---|
committer | Szabolcs Nagy <nsz@port70.net> | 2013-09-05 11:30:07 +0000 |
commit | ee2ee92d62c43f6658d37ddea4c316d2089d0fe9 (patch) | |
tree | 1a3e5a63fd40fa763f87d5dd9efd021117ea1d11 /src/math/modff.c | |
parent | d1a2ead878c27ac4ec600740320f8b76e1f961e9 (diff) | |
download | musl-ee2ee92d62c43f6658d37ddea4c316d2089d0fe9.tar.gz musl-ee2ee92d62c43f6658d37ddea4c316d2089d0fe9.tar.xz musl-ee2ee92d62c43f6658d37ddea4c316d2089d0fe9.zip |
math: rewrite remainder functions (remainder, remquo, fmod, modf)
* results are exact * modfl follows truncl (raises inexact flag spuriously now) * modf and modff only had cosmetic cleanup * remainder is just a wrapper around remquo now * using iterative shift+subtract for remquo and fmod * ld80 and ld128 are supported as well
Diffstat (limited to 'src/math/modff.c')
-rw-r--r-- | src/math/modff.c | 29 |
1 files changed, 13 insertions, 16 deletions
diff --git a/src/math/modff.c b/src/math/modff.c index 90caf316..639514ef 100644 --- a/src/math/modff.c +++ b/src/math/modff.c @@ -2,36 +2,33 @@ float modff(float x, float *iptr) { - union {float x; uint32_t n;} u = {x}; + union {float f; uint32_t i;} u = {x}; uint32_t mask; - int e; - - e = (int)(u.n>>23 & 0xff) - 0x7f; + int e = (int)(u.i>>23 & 0xff) - 0x7f; /* no fractional part */ if (e >= 23) { *iptr = x; - if (e == 0x80 && u.n<<9 != 0) { /* nan */ + if (e == 0x80 && u.i<<9 != 0) { /* nan */ return x; } - u.n &= 0x80000000; - return u.x; + u.i &= 0x80000000; + return u.f; } /* no integral part */ if (e < 0) { - u.n &= 0x80000000; - *iptr = u.x; + u.i &= 0x80000000; + *iptr = u.f; return x; } mask = 0x007fffff>>e; - if ((u.n & mask) == 0) { + if ((u.i & mask) == 0) { *iptr = x; - u.n &= 0x80000000; - return u.x; + u.i &= 0x80000000; + return u.f; } - u.n &= ~mask; - *iptr = u.x; - STRICT_ASSIGN(float, x, x - *iptr); - return x; + u.i &= ~mask; + *iptr = u.f; + return x - u.f; } |