TimeAveraging.H

Go to the documentation of this file.

#ifndef TIMEAVERAGING_H
#define TIMEAVERAGING_H
 
#include "amr-wind/core/Factory.H"
#include "amr-wind/utilities/PostProcessing.H"
 
#include "AMReX_Vector.H"
 
#include <limits>
#include <memory>
#include <string>
#include <map>
 
namespace amr_wind {
 
class CFDSim;
class SimTime;
class Field;
 
namespace averaging {
 
class FieldTimeAverage
    : public Factory<FieldTimeAverage, CFDSim&, const std::string&>
{
public:
    static std::string base_identifier() { return "FieldTimeAverage"; }
 
    ~FieldTimeAverage() override = default;
 
    virtual void operator()(
        const SimTime& time,
        const amrex::Real filter_width,
        const amrex::Real elapsed_time) = 0;
 
    virtual const std::string& average_field_name() = 0;
};
 
class TimeAveraging : public PostProcessBase::Register<TimeAveraging>
{
public:
    static std::string identifier() { return "TimeAveraging"; }
 
    TimeAveraging(CFDSim& /*sim*/, std::string /*label*/);
 
    ~TimeAveraging() override;
 
    void pre_init_actions() override;
 
    void initialize() override;
 
    void post_advance_work() override;
 
    void post_regrid_actions() override {}
 
    const std::string& add_averaging(
        const std::string& field_name,
        const std::string& avg_type = "ReAveraging");
 
private:
    CFDSim& m_sim;
 
    const std::string m_label;
 
    amrex::Vector<std::unique_ptr<FieldTimeAverage>> m_averages;
 
    std::map<std::string, FieldTimeAverage*> m_registered;
 
    amrex::Real m_start_time{0.0};
 
    amrex::Real m_stop_time{std::numeric_limits<amrex::Real>::max()};
 
    amrex::Real m_filter{0.0};
};
 
} // namespace averaging
} // namespace amr_wind
 
#endif /* TIMEAVERAGING_H */