From 257d08624ffb144f4a38812dd079349ca2acb526 Mon Sep 17 00:00:00 2001
From: Patrick Shriwise <pshriwise@gmail.com>
Date: Thu, 12 Nov 2020 12:24:01 -0600
Subject: [PATCH] Initialize one libMesh point locator per thread.

---
 src/mesh.cpp | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/src/mesh.cpp b/src/mesh.cpp
index 194d3e72fb..24a574ba4b 100644
--- a/src/mesh.cpp
+++ b/src/mesh.cpp
@@ -2121,9 +2121,15 @@ void LibMesh::initialize() {
   libMesh::ExplicitSystem& eq_sys =
     equation_systems_->add_system<libMesh::ExplicitSystem>(eq_system_name_);
 
-  PL_ = m_->sub_point_locator();
-  PL_->set_contains_point_tol(FP_COINCIDENT);
-  PL_->enable_out_of_mesh_mode();
+  #pragma omp parallel
+  {
+    #pragma omp critical
+    {
+    PL_ = m_->sub_point_locator();
+    }
+    PL_->set_contains_point_tol(FP_COINCIDENT);
+    PL_->enable_out_of_mesh_mode();
+  }
 
   first_element_ = *m_->elements_begin();