X-Git-Url: http://nsz.repo.hu/git/?p=musl;a=blobdiff_plain;f=src%2Fmath%2Fmodff.c;h=84d0b82a56250bc36e9ee7083f9f179e3e5ad90a;hp=bf6e4ced1a8d461946834c1ec26481db29106dcb;hb=f6ceccd92247575e4a35bc94f581a570b8052d43;hpb=9f58d06007818c43b9504e959ef81ab5f113b374 diff --git a/src/math/modff.c b/src/math/modff.c index bf6e4ced..84d0b82a 100644 --- a/src/math/modff.c +++ b/src/math/modff.c @@ -1,33 +1,37 @@ -#include "libm.h" +#include +#include float modff(float x, float *iptr) { - uint32_t u, mask; + union {float x; uint32_t n;} u = {x}; + uint32_t mask; int e; - GET_FLOAT_WORD(u, x); - e = (int)(u>>23 & 0xff) - 0x7f; + e = (int)(u.n>>23 & 0xff) - 0x7f; /* no fractional part */ if (e >= 23) { *iptr = x; - if (e == 0x80 && u<<9 != 0) /* nan */ + if (e == 0x80 && u.n<<9 != 0) { /* nan */ return x; - SET_FLOAT_WORD(x, u & 0x80000000); - return x; + } + u.n &= 0x80000000; + return u.x; } /* no integral part */ if (e < 0) { - SET_FLOAT_WORD(*iptr, u & 0x80000000); + u.n &= 0x80000000; + *iptr = u.x; return x; } mask = 0x007fffff>>e; - if ((u & mask) == 0) { + if ((u.n & mask) == 0) { *iptr = x; - SET_FLOAT_WORD(x, u & 0x80000000); - return x; + u.n &= 0x80000000; + return u.x; } - SET_FLOAT_WORD(*iptr, u & ~mask); + u.n &= ~mask; + *iptr = u.x; return x - *iptr; }