Fixed small bug

DiamonDinoia · DiamonDinoia · commit a44869704507 · 2024-06-28T20:26:56.000-04:00
diff --git a/src/spreadinterp.cpp b/src/spreadinterp.cpp
@@ -1276,7 +1276,7 @@ void interp_cube(FLT *FINUFFT_RESTRICT target, const FLT *du, const FLT *ker1,
         for (uint8_t dy{0}; dy < ns; ++dy) {
           const auto du_ptr = base_du_ptr + 2 * N1 * dy; // (see above)
           for (uint8_t l{0}; l < line_vectors; ++l) {
-            du_pts[l] = simd_type::load_unaligned(l * simd_size + base_du_ptr);
+            du_pts[l] = simd_type::load_unaligned(l * simd_size + du_ptr);
           }
           for (uint8_t l{0}; l < line_vectors; ++l) {
             line[l] = xsimd::fma(ker23_array[dy], du_pts[l], line[l]);

Original file line number	Diff line number	Diff line change
`@@ -1276,7 +1276,7 @@ void interp_cube(FLT FINUFFT_RESTRICT target, const FLT du, const FLT *ker1,`
`1276`	`1276`	`for (uint8_t dy{0}; dy < ns; ++dy) {`
`1277`	`1277`	`const auto du_ptr = base_du_ptr + 2 * N1 * dy; // (see above)`
`1278`	`1278`	`for (uint8_t l{0}; l < line_vectors; ++l) {`
`1279`		`- du_pts[l] = simd_type::load_unaligned(l * simd_size + base_du_ptr);`
	`1279`	`+ du_pts[l] = simd_type::load_unaligned(l * simd_size + du_ptr);`
`1280`	`1280`	`}`
`1281`	`1281`	`for (uint8_t l{0}; l < line_vectors; ++l) {`
`1282`	`1282`	`line[l] = xsimd::fma(ker23_array[dy], du_pts[l], line[l]);`