From 30634165980deaeead9f4e71e0ff6a5c9d897776 Mon Sep 17 00:00:00 2001 From: Ian Moffett Date: Mon, 14 Jul 2025 22:39:34 -0400 Subject: usr: libc: Add math.h + musl math impl port Signed-off-by: Ian Moffett --- lib/libc/src/musl-math/cosl.c | 39 +++++++++++++++++++++++++++++++++++++++ 1 file changed, 39 insertions(+) create mode 100644 lib/libc/src/musl-math/cosl.c (limited to 'lib/libc/src/musl-math/cosl.c') diff --git a/lib/libc/src/musl-math/cosl.c b/lib/libc/src/musl-math/cosl.c new file mode 100644 index 0000000..79c41c7 --- /dev/null +++ b/lib/libc/src/musl-math/cosl.c @@ -0,0 +1,39 @@ +#include "libm.h" + +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 +long double cosl(long double x) { + return cos(x); +} +#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 +long double cosl(long double x) +{ + union ldshape u = {x}; + unsigned n; + long double y[2], hi, lo; + + u.i.se &= 0x7fff; + if (u.i.se == 0x7fff) + return x - x; + x = u.f; + if (x < M_PI_4) { + if (u.i.se < 0x3fff - LDBL_MANT_DIG) + /* raise inexact if x!=0 */ + return 1.0 + x; + return __cosl(x, 0); + } + n = __rem_pio2l(x, y); + hi = y[0]; + lo = y[1]; + switch (n & 3) { + case 0: + return __cosl(hi, lo); + case 1: + return -__sinl(hi, lo, 1); + case 2: + return -__cosl(hi, lo); + case 3: + default: + return __sinl(hi, lo, 1); + } +} +#endif -- cgit v1.2.3