using fma where possible

DiamonDinoia · DiamonDinoia · commit 1c3e40e0324f · 2024-06-25T21:13:04.000-04:00
diff --git a/src/spreadinterp.cpp b/src/spreadinterp.cpp
@@ -832,7 +832,7 @@ void interp_line(FLT *FINUFFT_RESTRICT target, const FLT *du, const FLT *ker,
       out[1] = xsimd::fma(du[2 * j + 1], ker[dx], out[1]);
     }
   } else if (FINUFFT_UNLIKELY(i1 + ns >= N1)) { // wraps at right
-    for (int dx = 0; dx < N1 - i1; ++dx, ++j) {
+    for (UBIGINT dx = 0; dx < N1 - i1; ++dx, ++j) {
       out[0] = xsimd::fma(du[2 * j], ker[dx], out[0]);
       out[1] = xsimd::fma(du[2 * j + 1], ker[dx], out[1]);
     }

Original file line number	Diff line number	Diff line change
`@@ -832,7 +832,7 @@ void interp_line(FLT FINUFFT_RESTRICT target, const FLT du, const FLT *ker,`
`832`	`832`	`out[1] = xsimd::fma(du[2 * j + 1], ker[dx], out[1]);`
`833`	`833`	`}`
`834`	`834`	`} else if (FINUFFT_UNLIKELY(i1 + ns >= N1)) { // wraps at right`
`835`		`- for (int dx = 0; dx < N1 - i1; ++dx, ++j) {`
	`835`	`+ for (UBIGINT dx = 0; dx < N1 - i1; ++dx, ++j) {`
`836`	`836`	`out[0] = xsimd::fma(du[2 * j], ker[dx], out[0]);`
`837`	`837`	`out[1] = xsimd::fma(du[2 * j + 1], ker[dx], out[1]);`
`838`	`838`	`}`