Move workaround to math_fmod_impl

Modules/_math.h

@@ -23,19 +23,3 @@ _Py_log1p(double x)
 }
 
 #define m_log1p _Py_log1p
-
-static double
-_Py_fmod(double x, double y)
-{
-    /* Some platforms (e.g. Windows 10 with MSC v.1916) loose sign
-       for zero result.  But C99+ says: "if y is nonzero, the result
-       has the same sign as x".
-     */
-    double r = fmod(x, y);
-
-    if (r == 0.0 && y != 0.0) {
-        r = copysign(r, x);
-    }
-
-    return r;
-}

Modules/mathmodule.c

@@ -2347,7 +2347,14 @@ math_fmod_impl(PyObject *module, double x, double y)
     if (isinf(y) && isfinite(x))
         return PyFloat_FromDouble(x);
     errno = 0;
-    r = _Py_fmod(x, y);
+    r = fmod(x, y);
+    /* Some platforms (e.g. Windows 10 with MSC v.1916) loose sign
+       for zero result.  But C99+ says: "if y is nonzero, the result
+       has the same sign as x".
+     */
+    if (r == 0.0 && y != 0.0) {
+        r = copysign(r, x);
+    }
     if (isnan(r)) {
         if (!isnan(x) && !isnan(y))
             errno = EDOM;