msvcrt: Match native cosf results in -Pi/4 - Pi/4 range.

Signed-off-by: Piotr Caban <piotr@codeweavers.com> Signed-off-by: Alexandre Julliard <julliard@winehq.org>

msvcrt: Match native cosf results in -Pi/4 - Pi/4 range.
9008cd2f · Piotr Caban · Alexandre Julliard · ee7b5ebc · 9008cd2f
Commit 9008cd2f authored Jun 21, 2021 by Piotr Caban Committed by Alexandre Julliard Jun 21, 2021
Hide whitespace changes
Inline Side-by-side

Showing with 9 additions and 8 deletions

math.c dlls/msvcrt/math.c +9 -8

No files found.
--- a/dlls/msvcrt/math.c
+++ b/dlls/msvcrt/math.c
@@ -680,16 +680,17 @@ static float __sindf(double x)
 /* Copied from musl: src/math/__cosdf.c */
 static float __cosdf(double x)
 {
-    static const double C0 = -0x1ffffffd0c5e81.0p-54,
+    static const double C0 = -0x1.0000000000000p-1,
-        C1 = 0x155553e1053a42.0p-57,
+        C1 = 0x1.5555555555555p-5,
-        C2 = -0x16c087e80f1e27.0p-62,
+        C2 = -0x1.6c16c16c16c17p-10,
-        C3 = 0x199342e0ee5069.0p-68;
+        C3 = 0x1.a01a01a01a01ap-16,
-    double r, w, z;
+        C4 = -0x1.27e4fb7789f5cp-22;
+    double z;
    z = x * x;
-    w = z * z;
+    if (x > -7.8163146972656250e-03 && x < 7.8163146972656250e-03)
-    r = C2 + z * C3;
+        return 1 + C0 * z;
-    return ((1.0 + z * C0) + w * C1) + (w * z) * r;
+    return 1.0 + z * (C0 + z * (C1 + z * (C2 + z * (C3 + z * C4))));
 }
 static const UINT64 exp2f_T[] = {