Alexandre Julliard : msvcrt: Use the cosh() implementation from the bundled musl library.

[Alexandre Julliard]

Module: wine
Branch: master
Commit: 9fd9ebcb7a4a093990a8103cce6415b6e97299d7
URL:    https://gitlab.winehq.org/wine/wine/-/commit/9fd9ebcb7a4a093990a8103cce6415b...

Author: Alexandre Julliard 
Date:   Mon Apr  3 17:43:30 2023 +0200

msvcrt: Use the cosh() implementation from the bundled musl library.

With the changes from 60d178b410ba943c6c39b2ccf1deaec264d7d21a.

---

dlls/msvcrt/math.c        | 43 -------------------------------------------
 libs/musl/src/math/cosh.c |  5 +++++
 2 files changed, 5 insertions(+), 43 deletions(-)

diff --git a/dlls/msvcrt/math.c b/dlls/msvcrt/math.c
index 9741100b2b9..8bce6727ac4 100644
--- a/dlls/msvcrt/math.c
+++ b/dlls/msvcrt/math.c
@@ -2155,49 +2155,6 @@ double CDECL cos( double x )

extern double __expo2(double x, double sign);

-/*********************************************************************
- *		cosh (MSVCRT.@)
- *
- * Copied from musl: src/math/cosh.c
- */
-double CDECL cosh( double x )
-{
-    UINT64 ux = *(UINT64*)&x;
-    UINT64 sign = ux & 0x8000000000000000ULL;
-    UINT32 w;
-    double t;
-
-    /* |x| */
-    ux &= (uint64_t)-1 / 2;
-    x = *(double*)&ux;
-    w = ux >> 32;
-
-    /* |x| < log(2) */
-    if (w < 0x3fe62e42) {
-        if (w < 0x3ff00000 - (26 << 20)) {
-            fp_barrier(x + 0x1p120f);
-            return 1;
-        }
-        t = expm1(x);
-        return 1 + t * t / (2 * (1 + t));
-    }
-
-    /* |x| < log(DBL_MAX) */
-    if (w < 0x40862e42) {
-        t = exp(x);
-        /* note: if x>log(0x1p26) then the 1/t is not needed */
-        return 0.5 * (t + 1 / t);
-    }
-
-    /* |x| > log(DBL_MAX) or nan */
-    /* note: the result is stored to handle overflow */
-    if (ux > 0x7ff0000000000000ULL)
-        *(UINT64*)&t = ux | sign | 0x0008000000000000ULL;
-    else
-        t = __expo2(x, 1.0);
-    return t;
-}
-
 /* Copied from musl: src/math/exp_data.c */
 static const UINT64 exp_T[] = {
     0x0ULL, 0x3ff0000000000000ULL,
diff --git a/libs/musl/src/math/cosh.c b/libs/musl/src/math/cosh.c
index 384d9d8beac..8ab9417d2c2 100644
--- a/libs/musl/src/math/cosh.c
+++ b/libs/musl/src/math/cosh.c
@@ -7,6 +7,7 @@
 double __cdecl cosh(double x)
 {
    union {double f; uint64_t i;} u = {.f = x};
+        uint64_t sign = u.i & 0x8000000000000000ULL;
    uint32_t w;
    double t;

@@ -35,6 +36,10 @@ double __cdecl cosh(double x)

/* |x| > log(DBL_MAX) or nan */
    /* note: the result is stored to handle overflow */
+	if (w > 0x7ff00000) {
+		u.i |= sign | 0x0008000000000000ULL;
+		return u.f;
+	}
    t = __expo2(x, 1.0);
    return t;
 }