Fix OMP bug in AMReX-Codes#4238

WeiqunZhang · Dec 9, 2024 · e9f49de · e9f49de
1 parent 96db0a6
commit e9f49de
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/Src/EB/AMReX_EBFabFactory.cpp b/Src/EB/AMReX_EBFabFactory.cpp
@@ -22,6 +22,13 @@ EBFArrayBoxFactory::EBFArrayBoxFactory (const EB2::Level& a_level,
       m_parent(&a_level)
 {
     auto const& ebflags = getMultiEBCellFlagFab();
+    {
+        // If we do not do this here, there would a race condition when
+        // calling const_array() in getEBData, because const_arrays() is not
+        // thread safe.
+        auto const& ma = ebflags.const_arrays();
+        amrex::ignore_unused(ma);
+    }
 #ifdef AMREX_USE_GPU
     m_eb_data.resize(EBData::real_data_size*ebflags.local_size());
     Gpu::PinnedVector<Array4<Real const>> eb_data_hv;