Add rough laplace row major

src/helpers.rs

@@ -46,6 +46,49 @@ pub(crate) fn check_dimensions_evaluate<K: Kernel, T: RlstScalar>(
     );
 }
 
+
+pub(crate) fn check_dimensions_evaluate_row_major<K: Kernel, T: RlstScalar>(
+    kernel: &K,
+    eval_type: EvalType,
+    sources: &[T::Real],
+    targets: &[T::Real],
+    charges: &[T],
+    result: &[T],
+) {
+    assert!(
+        sources.len() % kernel.space_dimension() == 0,
+        "Length of sources {} is not a multiple of space dimension {}.",
+        sources.len(),
+        kernel.space_dimension()
+    );
+
+    assert!(
+        targets.len() % kernel.space_dimension() == 0,
+        "Length of targets {} is not a multiple of space dimension {}.",
+        sources.len(),
+        kernel.space_dimension()
+    );
+
+    let nsources = sources.len() / kernel.space_dimension();
+    let ntargets = targets.len() / kernel.space_dimension();
+
+    assert_eq!(
+        charges.len(),
+        nsources,
+        "Wrong dimension for `charges`. {} != {} ",
+        charges.len(),
+        nsources,
+    );
+
+    assert_eq!(
+        result.len(),
+        kernel.range_component_count(eval_type) * ntargets,
+        "Wrong dimension for `result`. {} != {} ",
+        result.len(),
+        ntargets * kernel.range_component_count(eval_type),
+    );
+}
+
 pub(crate) fn check_dimensions_assemble<K: Kernel, T: RlstScalar>(
     kernel: &K,
     eval_type: EvalType,

src/laplace_3d.rs

@@ -470,6 +470,8 @@ mod test {
             targets_y_eps[[index, 1]] += eps;
             targets_z_eps[[index, 2]] += eps;
         }
+        // println!("HERE FPP {:?} {:?} {:?} {:?}", targets_x_eps.data(), targets_x_eps[[0, 0]], targets_x_eps[[1, 0]], targets_x_eps[[2, 0]]);
+        // assert!(false);
 
         let mut expected = rlst_dynamic_array2!(f64, [4, ntargets]);