Clear the point locator after displacing the mesh

Displacing the mesh invalidates all the bounding boxes on the point locator, so we need to make sure we clear the locator so we don't use bad bounding box data Closes idaholab#22534
lindsayad · Nov 1, 2022 · 0c5b2a1 · 0c5b2a1
1 parent c2c6059
commit 0c5b2a1
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 0 deletions.
diff --git a/framework/include/loops/UpdateDisplacedMeshThread.h b/framework/include/loops/UpdateDisplacedMeshThread.h
@@ -34,6 +34,7 @@ class UpdateDisplacedMeshThread : public ThreadedNodeLoop<NodeRange, NodeRange::
   virtual void onNode(NodeRange::const_iterator & nd) override;
 
   void join(const UpdateDisplacedMeshThread & /*y*/);
+  void post() override;
 
 protected:
   void init();

diff --git a/framework/src/loops/UpdateDisplacedMeshThread.C b/framework/src/loops/UpdateDisplacedMeshThread.C
@@ -13,6 +13,7 @@
 #include "DisplacedProblem.h"
 #include "MooseMesh.h"
 
+#include "libmesh/mesh_base.h"
 #include "libmesh/numeric_vector.h"
 #include "libmesh/transient_system.h"
 #include "libmesh/explicit_system.h"
@@ -142,6 +143,12 @@ UpdateDisplacedMeshThread::onNode(NodeRange::const_iterator & nd)
   }
 }
 
+void
+UpdateDisplacedMeshThread::post()
+{
+  _displaced_problem.mesh().getMesh().clear_point_locator();
+}
+
 void
 UpdateDisplacedMeshThread::join(const UpdateDisplacedMeshThread & /*y*/)
 {