Evaluator works

bempp · Jan 15, 2025 · a9f6f61 · a9f6f61
1 parent 69ac00c
commit a9f6f61
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/examples/laplace_evaluator.rs b/examples/laplace_evaluator.rs
@@ -29,7 +29,7 @@ fn main() {
 
     let mut rng = ChaCha8Rng::seed_from_u64(world.rank() as u64);
 
-    let grid = bempp::shapes::regular_sphere::<f64, _>(6, 1, &world);
+    let grid = bempp::shapes::regular_sphere::<f64, _>(5, 1, &world);
 
     let quad_degree = 6;
     // Get the number of cells in the grid.

diff --git a/src/evaluator_tools.rs b/src/evaluator_tools.rs
@@ -25,7 +25,7 @@ use rlst::{
             distributed_sparse_operator::DistributedCsrMatrixOperatorImpl,
             DistributedArrayVectorSpace,
         },
-        Operator,
+        zero_element, Operator,
     },
     rlst_array_from_slice2, rlst_dynamic_array1, rlst_dynamic_array2, rlst_dynamic_array3,
     rlst_dynamic_array4, Array, AsApply, DefaultIterator, DistributedCsrMatrix, DistributedVector,
@@ -577,6 +577,8 @@ where
         &self,
         x: &<Self::Domain as rlst::LinearSpace>::E,
     ) -> <Self::Range as rlst::LinearSpace>::E {
-        todo!()
+        let mut y = zero_element(self.range());
+        self.apply_extended(T::one(), x, T::zero(), &mut y);
+        y
     }
 }