Add assert that we are not inside threaded region when InitRandom is …

…called. (#3204)
AMReX-Codes · Mar 15, 2023 · fdcbb7a · fdcbb7a
1 parent d54067a
commit fdcbb7a
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/Src/Base/AMReX_Random.cpp b/Src/Base/AMReX_Random.cpp
@@ -67,6 +67,12 @@ InitRandom (ULong cpu_seed, int nprocs, ULong gpu_seed)
     nthreads = OpenMP::get_max_threads();
     generators.resize(nthreads);
 
+#ifdef AMREX_USE_OMP
+    if (omp_in_parallel()) {
+        amrex::Abort("It is not safe to call amrex::InitRandom inside a threaded region.");
+    }
+#endif
+
 #ifdef AMREX_USE_OMP
 #pragma omp parallel
 #endif