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/rintf.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/rintf.c')
-rw-r--r-- | lib/libc/src/musl-math/rintf.c | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/lib/libc/src/musl-math/rintf.c b/lib/libc/src/musl-math/rintf.c new file mode 100644 index 0000000..9047688 --- /dev/null +++ b/lib/libc/src/musl-math/rintf.c @@ -0,0 +1,30 @@ +#include <float.h> +#include <math.h> +#include <stdint.h> + +#if FLT_EVAL_METHOD==0 +#define EPS FLT_EPSILON +#elif FLT_EVAL_METHOD==1 +#define EPS DBL_EPSILON +#elif FLT_EVAL_METHOD==2 +#define EPS LDBL_EPSILON +#endif +static const float_t toint = 1/EPS; + +float rintf(float x) +{ + union {float f; uint32_t i;} u = {x}; + int e = u.i>>23 & 0xff; + int s = u.i>>31; + float_t y; + + if (e >= 0x7f+23) + return x; + if (s) + y = x - toint + toint; + else + y = x + toint - toint; + if (y == 0) + return s ? -0.0f : 0.0f; + return y; +} |