| Index: fusl/src/math/modff.c
|
| diff --git a/fusl/src/math/modff.c b/fusl/src/math/modff.c
|
| index 639514effafcc91997d53f4b8528e16b21440a9e..eae990eb0f856051584de27aa2fae452a8d58f99 100644
|
| --- a/fusl/src/math/modff.c
|
| +++ b/fusl/src/math/modff.c
|
| @@ -1,34 +1,36 @@
|
| #include "libm.h"
|
|
|
| -float modff(float x, float *iptr)
|
| -{
|
| - union {float f; uint32_t i;} u = {x};
|
| - uint32_t mask;
|
| - int e = (int)(u.i>>23 & 0xff) - 0x7f;
|
| +float modff(float x, float* iptr) {
|
| + union {
|
| + float f;
|
| + uint32_t i;
|
| + } u = {x};
|
| + uint32_t mask;
|
| + int e = (int)(u.i >> 23 & 0xff) - 0x7f;
|
|
|
| - /* no fractional part */
|
| - if (e >= 23) {
|
| - *iptr = x;
|
| - if (e == 0x80 && u.i<<9 != 0) { /* nan */
|
| - return x;
|
| - }
|
| - u.i &= 0x80000000;
|
| - return u.f;
|
| - }
|
| - /* no integral part */
|
| - if (e < 0) {
|
| - u.i &= 0x80000000;
|
| - *iptr = u.f;
|
| - return x;
|
| - }
|
| + /* no fractional part */
|
| + if (e >= 23) {
|
| + *iptr = x;
|
| + if (e == 0x80 && u.i << 9 != 0) { /* nan */
|
| + return x;
|
| + }
|
| + u.i &= 0x80000000;
|
| + return u.f;
|
| + }
|
| + /* no integral part */
|
| + if (e < 0) {
|
| + u.i &= 0x80000000;
|
| + *iptr = u.f;
|
| + return x;
|
| + }
|
|
|
| - mask = 0x007fffff>>e;
|
| - if ((u.i & mask) == 0) {
|
| - *iptr = x;
|
| - u.i &= 0x80000000;
|
| - return u.f;
|
| - }
|
| - u.i &= ~mask;
|
| - *iptr = u.f;
|
| - return x - u.f;
|
| + mask = 0x007fffff >> e;
|
| + if ((u.i & mask) == 0) {
|
| + *iptr = x;
|
| + u.i &= 0x80000000;
|
| + return u.f;
|
| + }
|
| + u.i &= ~mask;
|
| + *iptr = u.f;
|
| + return x - u.f;
|
| }
|
|
|