| Index: fusl/src/math/modf.c
|
| diff --git a/fusl/src/math/modf.c b/fusl/src/math/modf.c
|
| index 1c8a1db90db7b4b61c5ec4271f184308690704db..01a1941e374b2d54526227e73b350b0a66e7522a 100644
|
| --- a/fusl/src/math/modf.c
|
| +++ b/fusl/src/math/modf.c
|
| @@ -1,34 +1,36 @@
|
| #include "libm.h"
|
|
|
| -double modf(double x, double *iptr)
|
| -{
|
| - union {double f; uint64_t i;} u = {x};
|
| - uint64_t mask;
|
| - int e = (int)(u.i>>52 & 0x7ff) - 0x3ff;
|
| +double modf(double x, double* iptr) {
|
| + union {
|
| + double f;
|
| + uint64_t i;
|
| + } u = {x};
|
| + uint64_t mask;
|
| + int e = (int)(u.i >> 52 & 0x7ff) - 0x3ff;
|
|
|
| - /* no fractional part */
|
| - if (e >= 52) {
|
| - *iptr = x;
|
| - if (e == 0x400 && u.i<<12 != 0) /* nan */
|
| - return x;
|
| - u.i &= 1ULL<<63;
|
| - return u.f;
|
| - }
|
| + /* no fractional part */
|
| + if (e >= 52) {
|
| + *iptr = x;
|
| + if (e == 0x400 && u.i << 12 != 0) /* nan */
|
| + return x;
|
| + u.i &= 1ULL << 63;
|
| + return u.f;
|
| + }
|
|
|
| - /* no integral part*/
|
| - if (e < 0) {
|
| - u.i &= 1ULL<<63;
|
| - *iptr = u.f;
|
| - return x;
|
| - }
|
| + /* no integral part*/
|
| + if (e < 0) {
|
| + u.i &= 1ULL << 63;
|
| + *iptr = u.f;
|
| + return x;
|
| + }
|
|
|
| - mask = -1ULL>>12>>e;
|
| - if ((u.i & mask) == 0) {
|
| - *iptr = x;
|
| - u.i &= 1ULL<<63;
|
| - return u.f;
|
| - }
|
| - u.i &= ~mask;
|
| - *iptr = u.f;
|
| - return x - u.f;
|
| + mask = -1ULL >> 12 >> e;
|
| + if ((u.i & mask) == 0) {
|
| + *iptr = x;
|
| + u.i &= 1ULL << 63;
|
| + return u.f;
|
| + }
|
| + u.i &= ~mask;
|
| + *iptr = u.f;
|
| + return x - u.f;
|
| }
|
|
|