From 885f5974cdf65b59415837ae97f5a14ef1350670 Mon Sep 17 00:00:00 2001 From: Kacper Date: Tue, 9 Dec 2025 19:20:15 +0100 Subject: feat: add gzip and new headers --- lib/libm/nexttowardf.c | 38 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 38 insertions(+) create mode 100644 lib/libm/nexttowardf.c (limited to 'lib/libm/nexttowardf.c') diff --git a/lib/libm/nexttowardf.c b/lib/libm/nexttowardf.c new file mode 100644 index 00000000..a342bc4c --- /dev/null +++ b/lib/libm/nexttowardf.c @@ -0,0 +1,38 @@ +#include "libm.h" + +float nexttowardf(float x, long double y) +{ + union { + float f; + uint32_t i; + } ux = { x }; + uint32_t e; + + if (isnan(x) || isnan(y)) + return x + y; + if (x == y) + return y; + if (x == 0) { + ux.i = 1; + if (signbit(y)) + ux.i |= 0x80000000; + } else if (x < y) { + if (signbit(x)) + ux.i--; + else + ux.i++; + } else { + if (signbit(x)) + 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; +} -- cgit v1.2.3