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 "amr-wind/utilities/constants.H"
 
#include "AMReX_Vector.H"
 
#include <limits>
#include <memory>
#include <string>
#include <map>
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace amr_wind {
 
class CFDSim;
class SimTime;
class Field;
 
namespace averaging {

class FieldTimeAverage
    : public Factory<
          FieldTimeAverage,
          CFDSim&,
          const std::string&,
          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 avg_time_interval,
        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 output_actions() 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_rt};

    amrex::Real m_stop_time{std::numeric_limits<amrex::Real>::max()};

    amrex::Real m_filter{amr_wind::constants::LARGE_NUM};

    amrex::Real m_time_interval{-1.};

    amrex::Real m_accumulated_avg_time_interval{0.};
};
 
} // namespace averaging
} // namespace amr_wind
 
#endif /* TIMEAVERAGING_H */