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 · Oct 28, 2022 · 11fae85 · 11fae85
1 parent e535795
commit 11fae85
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*/)
 {