diff options
Diffstat (limited to 'lib/libc/src/musl-math/frexpf.c')
-rw-r--r-- | lib/libc/src/musl-math/frexpf.c | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/lib/libc/src/musl-math/frexpf.c b/lib/libc/src/musl-math/frexpf.c new file mode 100644 index 0000000..2dabe37 --- /dev/null +++ b/lib/libc/src/musl-math/frexpf.c @@ -0,0 +1,24 @@ +#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) { + *e = 0; + return x; + } + + *e = ee - 0x7e; + y.i &= 0x807ffffful; + y.i |= 0x3f000000ul; + return y.f; +} |