Bug fix for bounding boxes with OpenMP

awslabs · Feb 5, 2024 · 1e55136 · 1e55136
1 parent 25aeaff
commit 1e55136
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/palace/utils/geodata.cpp b/palace/utils/geodata.cpp
@@ -588,6 +588,7 @@ void GetAxisAlignedBoundingBox(const mfem::ParMesh &mesh, const mfem::Array<int>
   }
   else
   {
+    mesh.GetNodes()->HostRead();
     const int ref = mesh.GetNodes()->FESpace()->GetMaxElementOrder();
     auto BBUpdate = [&ref](mfem::GeometryRefiner refiner, mfem::Geometry::Type &geom,
                            mfem::ElementTransformation &T, mfem::DenseMatrix &pointmat,