ResamplingTrigger.H

Go to the documentation of this file.

/* Copyright 2019-2020 Neil Zaim
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef WARPX_RESAMPLING_TRIGGER_H_
#define WARPX_RESAMPLING_TRIGGER_H_

#include "Utils/IntervalsParser.H"

#include <AMReX_ParmParse.H>

class ResamplingTrigger
{
public:

    ResamplingTrigger () = default;

    ResamplingTrigger (const std::string species_name);

    bool triggered (const int timestep, const amrex::Real global_numparts) const;

    void initialize_global_numcells () const;

private:
    // Intervals that define predetermined timesteps at which resampling is performed for all
    // species.
    IntervalsParser m_resampling_intervals;

    // Average number of particles per cell above which resampling is performed for a given species
    amrex::Real m_max_avg_ppc = std::numeric_limits<amrex::Real>::max();

    //Total number of simulated cells, summed over all mesh refinement levels.
    mutable amrex::Real m_global_numcells = amrex::Real(0.0);

    mutable bool m_initialized = false;
};

#endif //WARPX_RESAMPLING_TRIGGER_H_