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, 12 insertions, 1 deletions

diff --git a/Source/Particles/MultiParticleContainer.cpp b/Source/Particles/MultiParticleContainer.cpp
index d619a5de4..e27a1d1d3 100644
--- a/Source/Particles/MultiParticleContainer.cpp
+++ b/Source/Particles/MultiParticleContainer.cpp
@@ -695,6 +695,18 @@ MultiParticleContainer::doCoulombCollisions ()
     }
 }
 
+void MultiParticleContainer::doResampling (const int timestep)
+{
+    WARPX_PROFILE("MPC::doResampling");
+
+    for (auto& pc : allcontainers)
+    {
+        if (!pc->do_resampling){ continue; }
+
+        pc->resample(m_resampler, timestep);
+    }
+}
+
 void MultiParticleContainer::CheckIonizationProductSpecies()
 {
     for (int i=0; i < static_cast<int>(species_names.size()); i++){
@@ -1033,7 +1045,6 @@ MultiParticleContainer::doQEDSchwinger ()
 // Get cell volume multiplied by temporal step. In 2D the transverse size is
 // chosen by the user in the input file.
     amrex::Geometry const & geom = warpx.Geom(level_0);
-    auto domain_box = geom.Domain();
 #if (AMREX_SPACEDIM == 2)
     const auto dVdt = geom.CellSize(0) * geom.CellSize(1)
         * m_qed_schwinger_y_size * warpx.getdt(level_0);