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/floorl.c | 34 ++++++++++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) create mode 100644 lib/libc/src/musl-math/floorl.c (limited to 'lib/libc/src/musl-math/floorl.c') diff --git a/lib/libc/src/musl-math/floorl.c b/lib/libc/src/musl-math/floorl.c new file mode 100644 index 0000000..16aaec4 --- /dev/null +++ b/lib/libc/src/musl-math/floorl.c @@ -0,0 +1,34 @@ +#include "libm.h" + +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 +long double floorl(long double x) +{ + return floor(x); +} +#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 + +static const long double toint = 1/LDBL_EPSILON; + +long double floorl(long double x) +{ + union ldshape u = {x}; + int e = u.i.se & 0x7fff; + long double y; + + if (e >= 0x3fff+LDBL_MANT_DIG-1 || x == 0) + return x; + /* y = int(x) - x, where int(x) is an integer neighbor of x */ + if (u.i.se >> 15) + y = x - toint + toint - x; + else + y = x + toint - toint - x; + /* special case because of non-nearest rounding modes */ + if (e <= 0x3fff-1) { + FORCE_EVAL(y); + return u.i.se >> 15 ? -1 : 0; + } + if (y > 0) + return x + y - 1; + return x + y; +} +#endif -- cgit v1.2.3