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/nextafterf.c | |
| parent | 8f9e448b2ef6db7cd905540c21f3c5b190e7a1e7 (diff) | |
feat: add gzip and new headers
Diffstat (limited to 'lib/libm/nextafterf.c')
| -rw-r--r-- | lib/libm/nextafterf.c | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/lib/libm/nextafterf.c b/lib/libm/nextafterf.c new file mode 100644 index 00000000..912e1578 --- /dev/null +++ b/lib/libm/nextafterf.c @@ -0,0 +1,33 @@ +#include "libm.h" + +float nextafterf(float x, float y) +{ + union { + float f; + uint32_t i; + } ux = { x }, uy = { y }; + uint32_t ax, ay, e; + + if (isnan(x) || isnan(y)) + return x + y; + if (ux.i == uy.i) + return y; + ax = ux.i & 0x7fffffff; + ay = uy.i & 0x7fffffff; + if (ax == 0) { + if (ay == 0) + return y; + ux.i = (uy.i & 0x80000000) | 1; + } else if (ax > ay || ((ux.i ^ uy.i) & 0x80000000)) + ux.i--; + else + ux.i++; + e = ux.i & 0x7f800000; + /* raise overflow if ux.f is infinite and x is finite */ + if (e == 0x7f800000) + FORCE_EVAL(x + x); + /* raise underflow if ux.f is subnormal or zero */ + if (e == 0) + FORCE_EVAL(x * x + ux.f * ux.f); + return ux.f; +} |
