Refactor boundary assembly into better traits

benches/assembly_benchmark.rs

@@ -1,6 +1,6 @@
 use bempp::assembly::{boundary, boundary::BoundaryAssembler};
 use bempp::function::SerialFunctionSpace;
-use bempp::traits::FunctionSpace;
+use bempp::traits::{BoundaryAssembly, FunctionSpace};
 use criterion::{criterion_group, criterion_main, Criterion};
 use ndelement::ciarlet::LagrangeElementFamily;
 use ndelement::types::{Continuity, ReferenceCellType};

examples/assembly.rs

@@ -1,6 +1,6 @@
 use bempp::assembly::{boundary, boundary::BoundaryAssembler};
 use bempp::function::SerialFunctionSpace;
-use bempp::traits::FunctionSpace;
+use bempp::traits::{BoundaryAssembly, FunctionSpace};
 use ndelement::ciarlet::LagrangeElementFamily;
 use ndelement::types::{Continuity, ReferenceCellType};
 use ndgrid::shapes::regular_sphere;

examples/test_parallel_assembly.rs

@@ -5,9 +5,8 @@ use approx::assert_relative_eq;
 #[cfg(feature = "mpi")]
 use bempp::{
     assembly::boundary,
-    assembly::boundary::BoundaryAssembler,
     function::{ParallelFunctionSpace, SerialFunctionSpace},
-    traits::FunctionSpace,
+    traits::{BoundaryAssembly, FunctionSpace, ParallelBoundaryAssembly},
 };
 #[cfg(feature = "mpi")]
 use itertools::izip;

src/assembly.rs

@@ -2,14 +2,14 @@
 pub mod boundary;
 pub(crate) mod common;
 pub mod fmm_tools;
+pub mod kernels;
 pub mod potential;
 
 #[cfg(test)]
 mod test {
-    use super::boundary::BoundaryAssembler;
     use super::*;
     use crate::function::SerialFunctionSpace;
-    use crate::traits::FunctionSpace;
+    use crate::traits::{BoundaryAssembly, FunctionSpace};
     use cauchy::{c32, c64};
     use ndelement::ciarlet::CiarletElement;
     use ndelement::ciarlet::LagrangeElementFamily;
@@ -121,6 +121,16 @@ mod test {
            //};
     }
     macro_rules! create_assembler {
+        (Laplace, Hypersingular, $dtype:ident) => {
+            paste! {
+                boundary::HypersingularAssembler::<[<$dtype>], _, _>::new_laplace()
+            }
+        };
+        (Helmholtz, Hypersingular, $dtype:ident) => {
+            paste! {
+                boundary::HypersingularAssembler::<[<$dtype>], _, _>::new_helmholtz(3.0)
+            }
+        };
         (Laplace, $operator:ident, $dtype:ident) => {
             paste! {
                 boundary::[<$operator Assembler>]::<[<$dtype>], _>::new_laplace()