Add structure for resampling algorithms (#1265)

* First version of resampling structure

* Remove new warnings

* Added Doxygen comments

* Minor modifications

* Add do_resampling() in OneStep_sub1

* Apply suggestions from code review

Co-authored-by: Luca Fedeli <luca.fedeli.88@gmail.com>

* Use mutable in ResamplingTrigger and replace run with operator()

* Apply suggestions from code review

Co-authored-by: MaxThevenet <maxence.thevenet@desy.de>

* Make LevelingThinning operator() final

Co-authored-by: Luca Fedeli <luca.fedeli.88@gmail.com>
Co-authored-by: MaxThevenet <maxence.thevenet@desy.de>

1 files changed, 20 insertions, 0 deletions

diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
index e3aaa6002..1b4c243ac 100644
--- a/Source/Particles/PhysicalParticleContainer.cpp
+++ b/Source/Particles/PhysicalParticleContainer.cpp
@@ -122,6 +122,8 @@ PhysicalParticleContainer::PhysicalParticleContainer (AmrCore* amr_core, int isp
 
     pp.query("do_field_ionization", do_field_ionization);
 
+    pp.query("do_resampling", do_resampling);
+
     //check if Radiation Reaction is enabled and do consistency checks
     pp.query("do_classical_radiation_reaction", do_classical_radiation_reaction);
     //if the species is not a lepton, do_classical_radiation_reaction
@@ -1988,6 +1990,24 @@ PhysicalParticleContainer::getIonizationFunc (const WarpXParIter& pti,
                                 ion_atomic_number);
 }
 
+void PhysicalParticleContainer::resample (const Resampling& resampler, const int timestep)
+{
+    const amrex::Real global_numparts = TotalNumberOfParticles();
+
+    if (resampler.triggered(timestep, global_numparts))
+    {
+        for (int lev = 0; lev <= maxLevel(); lev++)
+        {
+            for (WarpXParIter pti(*this, lev); pti.isValid(); ++pti)
+            {
+                resampler(pti);
+            }
+        }
+    }
+
+}
+
+
 #ifdef WARPX_QED