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/rint.c | 28 ++++++++++++++++++++++++++++ 1 file changed, 28 insertions(+) create mode 100644 lib/libc/src/musl-math/rint.c (limited to 'lib/libc/src/musl-math/rint.c') diff --git a/lib/libc/src/musl-math/rint.c b/lib/libc/src/musl-math/rint.c new file mode 100644 index 0000000..fbba390 --- /dev/null +++ b/lib/libc/src/musl-math/rint.c @@ -0,0 +1,28 @@ +#include +#include +#include + +#if FLT_EVAL_METHOD==0 || FLT_EVAL_METHOD==1 +#define EPS DBL_EPSILON +#elif FLT_EVAL_METHOD==2 +#define EPS LDBL_EPSILON +#endif +static const double_t toint = 1/EPS; + +double rint(double x) +{ + union {double f; uint64_t i;} u = {x}; + int e = u.i>>52 & 0x7ff; + int s = u.i>>63; + double_t y; + + if (e >= 0x3ff+52) + return x; + if (s) + y = x - toint + toint; + else + y = x + toint - toint; + if (y == 0) + return s ? -0.0 : 0; + return y; +} -- cgit v1.2.3