1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
#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; }