| #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; |
| |
| /* 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; |
| } |