CLN: remove cmath, closes pandas-dev#23209

Author: jbrockmendel

pandas/_libs/src/headers/cmath

@@ -9,15 +9,11 @@ from libc.stdlib cimport malloc, free
 import numpy as np
 cimport numpy as cnp
 from numpy cimport ndarray, int64_t, float64_t, float32_t
+from numpy.math cimport signbit, isnan, sqrtl as sqrt
+# sqrtl is square root specialized to "long double"
 cnp.import_array()
 
 
-cdef extern from "src/headers/cmath" namespace "std":
-    bint isnan(float64_t) nogil
-    bint notnan(float64_t) nogil
-    int signbit(float64_t) nogil
-    float64_t sqrt(float64_t x) nogil
-
 cimport pandas._libs.util as util
 from pandas._libs.util cimport numeric
 
@@ -384,21 +380,21 @@ def roll_count(ndarray[float64_t] values, int64_t win, int64_t minp,
                 count_x = 0.0
                 for j in range(s, e):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         count_x += 1.0
 
             else:
 
                 # calculate deletes
                 for j in range(start[i - 1], s):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         count_x -= 1.0
 
                 # calculate adds
                 for j in range(end[i - 1], e):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         count_x += 1.0
 
             if count_x >= minp:
@@ -429,15 +425,15 @@ cdef inline void add_sum(float64_t val, int64_t *nobs, float64_t *sum_x) nogil:
     """ add a value from the sum calc """
 
     # Not NaN
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] + 1
         sum_x[0] = sum_x[0] + val
 
 
 cdef inline void remove_sum(float64_t val, int64_t *nobs, float64_t *sum_x) nogil:
     """ remove a value from the sum calc """
 
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] - 1
         sum_x[0] = sum_x[0] - val
 
@@ -545,7 +541,7 @@ cdef inline void add_mean(float64_t val, Py_ssize_t *nobs, float64_t *sum_x,
     """ add a value from the mean calc """
 
     # Not NaN
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] + 1
         sum_x[0] = sum_x[0] + val
         if signbit(val):
@@ -556,7 +552,7 @@ cdef inline void remove_mean(float64_t val, Py_ssize_t *nobs, float64_t *sum_x,
                              Py_ssize_t *neg_ct) nogil:
     """ remove a value from the mean calc """
 
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] - 1
         sum_x[0] = sum_x[0] - val
         if signbit(val):
@@ -683,7 +679,7 @@ cdef inline void remove_var(float64_t val, float64_t *nobs, float64_t *mean_x,
     cdef:
         float64_t delta
 
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] - 1
         if nobs[0]:
             # a part of Welford's method for the online variance-calculation
@@ -772,7 +768,7 @@ def roll_var(ndarray[float64_t] values, int64_t win, int64_t minp,
                 val = values[i]
                 prev = values[i - win]
 
-                if notnan(val):
+                if not isnan(val):
                     if prev == prev:
 
                         # Adding one observation and removing another one
@@ -837,7 +833,7 @@ cdef inline void add_skew(float64_t val, int64_t *nobs,
     """ add a value from the skew calc """
 
     # Not NaN
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] + 1
 
         # seriously don't ask me why this is faster
@@ -852,7 +848,7 @@ cdef inline void remove_skew(float64_t val, int64_t *nobs,
     """ remove a value from the skew calc """
 
     # Not NaN
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] - 1
 
         # seriously don't ask me why this is faster
@@ -979,7 +975,7 @@ cdef inline void add_kurt(float64_t val, int64_t *nobs,
     """ add a value from the kurotic calc """
 
     # Not NaN
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] + 1
 
         # seriously don't ask me why this is faster
@@ -995,7 +991,7 @@ cdef inline void remove_kurt(float64_t val, int64_t *nobs,
     """ remove a value from the kurotic calc """
 
     # Not NaN
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] - 1
 
         # seriously don't ask me why this is faster
@@ -1116,7 +1112,7 @@ def roll_median_c(ndarray[float64_t] values, int64_t win, int64_t minp,
                 # setup
                 for j in range(s, e):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         nobs += 1
                         err = skiplist_insert(sl, val) != 1
                         if err:
@@ -1127,7 +1123,7 @@ def roll_median_c(ndarray[float64_t] values, int64_t win, int64_t minp,
                 # calculate adds
                 for j in range(end[i - 1], e):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         nobs += 1
                         err = skiplist_insert(sl, val) != 1
                         if err:
@@ -1136,7 +1132,7 @@ def roll_median_c(ndarray[float64_t] values, int64_t win, int64_t minp,
                 # calculate deletes
                 for j in range(start[i - 1], s):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         skiplist_remove(sl, val)
                         nobs -= 1
 
@@ -1510,7 +1506,7 @@ def roll_quantile(ndarray[float64_t, cast=True] values, int64_t win,
                 # setup
                 for j in range(s, e):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         nobs += 1
                         skiplist_insert(skiplist, val)
 
@@ -1519,14 +1515,14 @@ def roll_quantile(ndarray[float64_t, cast=True] values, int64_t win,
                 # calculate adds
                 for j in range(end[i - 1], e):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         nobs += 1
                         skiplist_insert(skiplist, val)
 
                 # calculate deletes
                 for j in range(start[i - 1], s):
                     val = values[j]
-                    if notnan(val):
+                    if not isnan(val):
                         skiplist_remove(skiplist, val)
                         nobs -= 1
 
@@ -1885,7 +1881,7 @@ cdef inline void remove_weighted_var(float64_t val,
     cdef:
         float64_t temp, q, r
 
-    if notnan(val):
+    if not isnan(val):
         nobs[0] = nobs[0] - 1
 
         if nobs[0]:
@@ -1955,7 +1951,7 @@ def roll_weighted_var(float64_t[:] values, float64_t[:] weights,
             w = weights[i % win_n]
             pre_w = weights[(i - win_n) % win_n]
 
-            if notnan(val):
+            if not isnan(val):
                 if pre_val == pre_val:
                     remove_weighted_var(pre_val, pre_w, &t,
                                         &sum_w, &mean, &nobs)