| #include <math.h> | |
| #include <stdint.h> | |
| float frexpf(float x, int* e) { | |
| union { | |
| float f; | |
| uint32_t i; | |
| } y = {x}; | |
| int ee = y.i >> 23 & 0xff; | |
| if (!ee) { | |
| if (x) { | |
| x = frexpf(x * 0x1p64, e); | |
| *e -= 64; | |
| } else | |
| *e = 0; | |
| return x; | |
| } else if (ee == 0xff) { | |
| return x; | |
| } | |
| *e = ee - 0x7e; | |
| y.i &= 0x807ffffful; | |
| y.i |= 0x3f000000ul; | |
| return y.f; | |
| } |