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/ceilf.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/ceilf.c')
-rw-r--r-- | lib/libc/src/musl-math/ceilf.c | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/lib/libc/src/musl-math/ceilf.c b/lib/libc/src/musl-math/ceilf.c new file mode 100644 index 0000000..869835f --- /dev/null +++ b/lib/libc/src/musl-math/ceilf.c @@ -0,0 +1,27 @@ +#include "libm.h" + +float ceilf(float x) +{ + union {float f; uint32_t i;} u = {x}; + int e = (int)(u.i >> 23 & 0xff) - 0x7f; + uint32_t m; + + if (e >= 23) + return x; + if (e >= 0) { + m = 0x007fffff >> e; + if ((u.i & m) == 0) + return x; + FORCE_EVAL(x + 0x1p120f); + if (u.i >> 31 == 0) + u.i += m; + u.i &= ~m; + } else { + FORCE_EVAL(x + 0x1p120f); + if (u.i >> 31) + u.f = -0.0; + else if (u.i << 1) + u.f = 1.0; + } + return u.f; +} |