diff options
| author | Kacper <kacper@mail.openlinux.dev> | 2025-12-09 19:20:15 +0100 |
|---|---|---|
| committer | Kacper <kacper@mail.openlinux.dev> | 2025-12-09 19:20:15 +0100 |
| commit | 885f5974cdf65b59415837ae97f5a14ef1350670 (patch) | |
| tree | 66ac13de29c7f4932c5fcae11773df574e4e256a /lib/libm/rintf.c | |
| parent | 8f9e448b2ef6db7cd905540c21f3c5b190e7a1e7 (diff) | |
feat: add gzip and new headers
Diffstat (limited to 'lib/libm/rintf.c')
| -rw-r--r-- | lib/libm/rintf.c | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/lib/libm/rintf.c b/lib/libm/rintf.c new file mode 100644 index 00000000..5db9e0fe --- /dev/null +++ b/lib/libm/rintf.c @@ -0,0 +1,33 @@ +#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; +} |
