diff options
author | Ian Moffett <ian@osmora.org> | 2025-07-14 22:39:34 -0400 |
---|---|---|
committer | Ian Moffett <ian@osmora.org> | 2025-07-14 22:39:34 -0400 |
commit | 7450f4f2bd19b1c6de0218a44aaba53639985768 (patch) | |
tree | 8b63e43a8a64810154b05c1176bd547a20e99d60 /lib/libc/src/musl-math/floorl.c | |
parent | 83acdbef7ae961c2ae296835df3f368573c89eeb (diff) |
usr: libc: Add math.h + musl math impl port
Signed-off-by: Ian Moffett <ian@osmora.org>
Diffstat (limited to 'lib/libc/src/musl-math/floorl.c')
-rw-r--r-- | lib/libc/src/musl-math/floorl.c | 34 |
1 files changed, 34 insertions, 0 deletions
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 |