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/roundl.c | 37 +++++++++++++++++++++++++++++++++++++ 1 file changed, 37 insertions(+) create mode 100644 lib/libm/roundl.c (limited to 'lib/libm/roundl.c') diff --git a/lib/libm/roundl.c b/lib/libm/roundl.c new file mode 100644 index 00000000..cb148310 --- /dev/null +++ b/lib/libm/roundl.c @@ -0,0 +1,37 @@ +#include "libm.h" + +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 +long double roundl(long double x) +{ + return round(x); +} +#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 + +static const long double toint = 1 / LDBL_EPSILON; + +long double roundl(long double x) +{ + union ldshape u = { x }; + int e = u.i.se & 0x7fff; + long double y; + + if (e >= 0x3fff + LDBL_MANT_DIG - 1) + return x; + if (u.i.se >> 15) + x = -x; + if (e < 0x3fff - 1) { + FORCE_EVAL(x + toint); + return 0 * u.f; + } + y = x + toint - toint - x; + if (y > 0.5) + y = y + x - 1; + else if (y <= -0.5) + y = y + x + 1; + else + y = y + x; + if (u.i.se >> 15) + y = -y; + return y; +} +#endif -- cgit v1.2.3