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

dlls/msvcrt/math.c

@@ -983,8 +983,6 @@ float CDECL cosf( float x )
     }
 }
 
-extern float __expo2f(float x, float sign);
-
 /*********************************************************************
  *      expf (MSVCRT.@)
  */
@@ -1277,40 +1275,6 @@ float CDECL sinf( float x )
     }
 }
 
-/*********************************************************************
- *      sinhf (MSVCRT.@)
- */
-float CDECL sinhf( float x )
-{
-    UINT32 ui = *(UINT32*)&x;
-    float t, h, absx;
-
-    h = 0.5;
-    if (ui >> 31)
-        h = -h;
-    /* |x| */
-    ui &= 0x7fffffff;
-    absx = *(float*)&ui;
-
-    /* |x| < log(FLT_MAX) */
-    if (ui < 0x42b17217) {
-        t = expm1f(absx);
-        if (ui < 0x3f800000) {
-            if (ui < 0x3f800000 - (12 << 23))
-                return x;
-            return h * (2 * t - t * t / (t + 1));
-        }
-        return h * (t + t / (t + 1));
-    }
-
-    /* |x| > logf(FLT_MAX) or nan */
-    if (ui > 0x7f800000)
-        *(DWORD*)&t = *(DWORD*)&x | 0x400000;
-    else
-        t = __expo2f(absx, 2 * h);
-    return t;
-}
-
 static BOOL sqrtf_validate( float *x )
 {
     short c = _fdclass(*x);

libs/musl/src/math/sinhf.c

@@ -3,6 +3,7 @@
 float __cdecl sinhf(float x)
 {
 	union {float f; uint32_t i;} u = {.f = x};
+	uint32_t sign = u.i & 0x80000000;
 	uint32_t w;
 	float t, h, absx;
 
@@ -26,6 +27,10 @@ float __cdecl sinhf(float x)
 	}
 
 	/* |x| > logf(FLT_MAX) or nan */
+	if (w > 0x7f800000) {
+            u.i = w | sign | 0x400000;
+            return u.f;
+	}
 	t = __expo2f(absx, 2*h);
 	return t;
 }