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

[Alexandre Julliard]

Module: wine
Branch: master
Commit: e525730cc1d929ebd17c99ceb704476ea7f693ff
URL:    https://gitlab.winehq.org/wine/wine/-/commit/e525730cc1d929ebd17c99ceb704476...

Author: Alexandre Julliard 
Date:   Fri Apr  7 14:44:04 2023 +0200

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

---

dlls/msvcrt/math.c      | 81 +++----------------------------------------------
 dlls/msvcrt/msvcrt.spec |  2 +-
 2 files changed, 5 insertions(+), 78 deletions(-)

diff --git a/dlls/msvcrt/math.c b/dlls/msvcrt/math.c
index d8eab014e42..b44bc051453 100644
--- a/dlls/msvcrt/math.c
+++ b/dlls/msvcrt/math.c
@@ -292,87 +292,14 @@ int CDECL _isnanf( float num )

/*********************************************************************
  *      atanf (MSVCRT.@)
- *
- * Copied from musl: src/math/atanf.c
  */
-float CDECL atanf( float x )
+#if _MSVCR_VER == 0  /* other versions call atanf() directly */
+float CDECL MSVCRT_atanf( float x )
 {
-    static const float atanhi[] = {
-        4.6364760399e-01,
-        7.8539812565e-01,
-        9.8279368877e-01,
-        1.5707962513e+00,
-    };
-    static const float atanlo[] = {
-        5.0121582440e-09,
-        3.7748947079e-08,
-        3.4473217170e-08,
-        7.5497894159e-08,
-    };
-    static const float aT[] = {
-        3.3333328366e-01,
-        -1.9999158382e-01,
-        1.4253635705e-01,
-        -1.0648017377e-01,
-        6.1687607318e-02,
-    };
-
-    float w, s1, s2, z;
-    unsigned int ix, sign;
-    int id;
-
-#if _MSVCR_VER == 0
     if (isnan(x)) return math_error(_DOMAIN, "atanf", x, 0, x);
-#endif
-
-    ix = *(unsigned int*)&x;
-    sign = ix >> 31;
-    ix &= 0x7fffffff;
-    if (ix >= 0x4c800000) {  /* if |x| >= 2**26 */
-        if (isnan(x))
-            return x;
-        z = atanhi[3] + 7.5231638453e-37;
-        return sign ? -z : z;
-    }
-    if (ix < 0x3ee00000) {   /* |x| < 0.4375 */
-        if (ix < 0x39800000) {  /* |x| < 2**-12 */
-            if (ix < 0x00800000)
-                /* raise underflow for subnormal x */
-                fp_barrierf(x*x);
-            return x;
-        }
-        id = -1;
-    } else {
-        x = fabsf(x);
-        if (ix < 0x3f980000) {  /* |x| < 1.1875 */
-            if (ix < 0x3f300000) {  /*  7/16 <= |x| < 11/16 */
-                id = 0;
-                x = (2.0f * x - 1.0f) / (2.0f + x);
-            } else {                /* 11/16 <= |x| < 19/16 */
-                id = 1;
-                x = (x - 1.0f) / (x + 1.0f);
-            }
-        } else {
-            if (ix < 0x401c0000) {  /* |x| < 2.4375 */
-                id = 2;
-                x = (x - 1.5f) / (1.0f + 1.5f * x);
-            } else {                /* 2.4375 <= |x| < 2**26 */
-                id = 3;
-                x = -1.0f / x;
-            }
-        }
-    }
-    /* end of argument reduction */
-    z = x * x;
-    w = z * z;
-    /* break sum from i=0 to 10 aT[i]z**(i+1) into odd and even poly */
-    s1 = z * (aT[0] + w * (aT[2] + w * aT[4]));
-    s2 = w * (aT[1] + w * aT[3]);
-    if (id < 0)
-        return x - x * (s1 + s2);
-    z = atanhi[id] - ((x * (s1 + s2) - atanlo[id]) - x);
-    return sign ? -z : z;
+    return atanf( x );
 }
+#endif

static BOOL sqrtf_validate( float *x )
 {
diff --git a/dlls/msvcrt/msvcrt.spec b/dlls/msvcrt/msvcrt.spec
index 0b88ba97135..8bf1c87d73c 100644
--- a/dlls/msvcrt/msvcrt.spec
+++ b/dlls/msvcrt/msvcrt.spec
@@ -1257,7 +1257,7 @@
 @ cdecl atan(double)
 @ cdecl atan2(double double)
 @ cdecl -arch=!i386 asinf(float)
-@ cdecl -arch=!i386 atanf(float)
+@ cdecl -arch=!i386 atanf(float) MSVCRT_atanf
 @ cdecl -arch=!i386 atan2f(float float)
 @ cdecl -private atexit(ptr) MSVCRT_atexit  # not imported to avoid conflicts with Mingw
 @ cdecl atof(str)