From 30634165980deaeead9f4e71e0ff6a5c9d897776 Mon Sep 17 00:00:00 2001 From: Ian Moffett Date: Mon, 14 Jul 2025 22:39:34 -0400 Subject: usr: libc: Add math.h + musl math impl port Signed-off-by: Ian Moffett --- lib/libc/src/musl-math/frexpl.c | 30 ++++++++++++++++++++++++++++++ 1 file changed, 30 insertions(+) create mode 100644 lib/libc/src/musl-math/frexpl.c (limited to 'lib/libc/src/musl-math/frexpl.c') diff --git a/lib/libc/src/musl-math/frexpl.c b/lib/libc/src/musl-math/frexpl.c new file mode 100644 index 0000000..e05cf75 --- /dev/null +++ b/lib/libc/src/musl-math/frexpl.c @@ -0,0 +1,30 @@ +#include "libm.h" + +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 +long double frexpl(long double x, int *e) +{ + return frexp(x, e); +} +#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 +long double frexpl(long double x, int *e) +{ + union ldshape u = {x}; + int ee = u.i.se & 0x7fff; + + if (!ee) { + if (x) { + x = frexpl(x*0x1p120, e); + *e -= 120; + } else *e = 0; + return x; + } else if (ee == 0x7fff) { + *e = 0; + return x; + } + + *e = ee - 0x3ffe; + u.i.se &= 0x8000; + u.i.se |= 0x3ffe; + return u.f; +} +#endif -- cgit v1.2.3