diff options
| author | Ian Moffett <ian@osmora.org> | 2025-07-14 22:39:34 -0400 |
|---|---|---|
| committer | Ian Moffett <ian@osmora.org> | 2025-07-14 23:27:58 -0400 |
| commit | 30634165980deaeead9f4e71e0ff6a5c9d897776 (patch) | |
| tree | 18546a9d493a699eb3ea23ab5cc792d40cf82d7d /lib/libc/src/musl-math/scalbnf.c | |
| parent | 83acdbef7ae961c2ae296835df3f368573c89eeb (diff) | |
usr: libc: Add math.h + musl math impl portexpt
Signed-off-by: Ian Moffett <ian@osmora.org>
Diffstat (limited to 'lib/libc/src/musl-math/scalbnf.c')
| -rw-r--r-- | lib/libc/src/musl-math/scalbnf.c | 31 |
1 files changed, 31 insertions, 0 deletions
diff --git a/lib/libc/src/musl-math/scalbnf.c b/lib/libc/src/musl-math/scalbnf.c new file mode 100644 index 0000000..a5ad208 --- /dev/null +++ b/lib/libc/src/musl-math/scalbnf.c @@ -0,0 +1,31 @@ +#include <math.h> +#include <stdint.h> + +float scalbnf(float x, int n) +{ + union {float f; uint32_t i;} u; + float_t y = x; + + if (n > 127) { + y *= 0x1p127f; + n -= 127; + if (n > 127) { + y *= 0x1p127f; + n -= 127; + if (n > 127) + n = 127; + } + } else if (n < -126) { + y *= 0x1p-126f * 0x1p24f; + n += 126 - 24; + if (n < -126) { + y *= 0x1p-126f * 0x1p24f; + n += 126 - 24; + if (n < -126) + n = -126; + } + } + u.i = (uint32_t)(0x7f+n)<<23; + x = y * u.f; + return x; +} |