add RPE

Author: tjhei

include/aspect/mesh_deformation/external_tool_interface.h

@@ -181,9 +181,7 @@ namespace aspect
 
         std::vector<Point<dim>> evaluation_points;
 
-        // Timo: Replace by whatever type you need here
-        class X {};
-        std::unique_ptr<X> remote_point_evaluator;
+        std::unique_ptr<Utilities::MPI::RemotePointEvaluation<dim, dim>> remote_point_evaluator;
     };
   }
 }

source/mesh_deformation/external_tool_interface.cc

@@ -83,12 +83,11 @@ namespace aspect
       // for later function calls describe the same number of points.
       this->evaluation_points = evaluation_points;
 
-      // Then also invalidate the previous evaluator:
-      remote_point_evaluator.reset();
-
-      // Timo:
-      // TODO: Implement set-up phase via MPIRemotePointEvaluation
-
+      // Set up RemotePointEvaluation with a Mapping of the undeformed mesh:
+      remote_point_evaluator = std::make_unique<Utilities::MPI::RemotePointEvaluation<dim, dim>>();
+      // TODO: does this need to be a higher order mapping for spherical problems?
+      static MappingQ<dim> mapping(1);
+      remote_point_evaluator->reinit(this->evaluation_points, this->get_triangulation(), mapping);
 
       // Finally, also ensure that upon mesh refinement, all of the
       // information set herein is invalidated:
@@ -97,8 +96,7 @@ namespace aspect
       {
         this->evaluation_points.clear();
         this->remote_point_evaluator.reset();
-      }
-              );
+      });
     }