New Diagnostics: FieldMaximum (#988)

* Add MaxField reduced diagnostic

* Add test of the MaxField reduced diag

* Minor modifications

* Fix test input file

* Use read_raw_data.py module in test

* Merge FieldEnergy and MaxField reduced diags

* Add specific error message if old syntax is used

* Fix column header numbers

* Update copyright

* Change name of diag in input file from Field to FieldReduced

* Add comment when allocating MultiFAB

* Add BackwardCompatibility () function

* Update Source/Diagnostics/ReducedDiags/MultiReducedDiags.cpp

Co-authored-by: MaxThevenet <mthevenet@lbl.gov>

* Add FieldReduced.cpp in CMakeLists.txt

* Apply suggestions from code review

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Remove magic numbers and update MFIter construct

* Apply suggestions from code review

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Fix typo

* Do a reduction on the fly instead of allocating a MultiFAB

* Interpolate fields to cell center before computing maximum value of |E| and |B|

* Reseparate FieldReduced into FieldEnergy and MaxField

* Fix few typos

* Replace MaxField with FieldMaximum

Co-authored-by: MaxThevenet <mthevenet@lbl.gov>
Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

1 files changed, 3 insertions, 9 deletions

diff --git a/Source/Particles/MultiParticleContainer.cpp b/Source/Particles/MultiParticleContainer.cpp
index 9bcb19b2e..ec4abdea9 100644
--- a/Source/Particles/MultiParticleContainer.cpp
+++ b/Source/Particles/MultiParticleContainer.cpp
@@ -1088,17 +1088,11 @@ MultiParticleContainer::doQEDSchwinger ()
     const MultiFab & By = warpx.getBfield(level_0,1);
     const MultiFab & Bz = warpx.getBfield(level_0,2);
 
-    MFItInfo info;
-    if (TilingIfNotGPU()) {
-        info.EnableTiling();
-    }
 #ifdef _OPENMP
-    info.SetDynamic(WarpX::do_dynamic_scheduling);
-#pragma omp parallel
+#pragma omp parallel if (amrex::Gpu::notInLaunchRegion())
 #endif
-
-    for (MFIter mfi(Ex, info); mfi.isValid(); ++mfi )
-    {
+     for (MFIter mfi(Ex, TilingIfNotGPU()); mfi.isValid(); ++mfi )
+     {
         // Make the box cell centered to avoid creating particles twice on the tile edges
         const Box& box = enclosedCells(mfi.nodaltilebox());