WIP

Author: DiamonDinoia

src/spreadinterp.cpp

@@ -274,7 +274,7 @@ static void evaluate_kernel_vector(T *ker, T *args,
   }
   // Separate check from arithmetic (Is this really needed? doesn't slow down)
   for (int i = 0; i < N; i++)
-    if (abs(args[i]) >= (T)opts.ES_halfwidth) ker[i] = 0.0;
+    if (std::abs(args[i]) >= (T)opts.ES_halfwidth) ker[i] = 0.0;
 }
 
 template<typename T, uint8_t w, uint8_t upsampfact,
@@ -2195,7 +2195,7 @@ T evaluate_kernel(T x, const finufft_spread_opts &opts)
    Rescaled so max is 1, Barnett 7/21/24
 */
 {
-  if (abs(x) >= (T)opts.ES_halfwidth)
+  if (std::abs(x) >= (T)opts.ES_halfwidth)
     // if spreading/FT careful, shouldn't need this if, but causes no speed hit
     return 0.0;
   else
@@ -2212,7 +2212,7 @@ T evaluate_kernel_horner_kernel(T x, const finufft_spread_opts &opts)
    using generated piecewise polynomial approximation to the kernel.
 */
 {
-  if (abs(x) >= (T)opts.ES_halfwidth) {
+  if (std::abs(x) >= (T)opts.ES_halfwidth) {
     // if spreading/FT careful, shouldn't need this if, but causes no speed hit
     return 0.0;
   } else {