PlaneSampler.H

Go to the documentation of this file.

#ifndef PLANESAMPLER_H
#define PLANESAMPLER_H
 
#include "src/utilities/sampling/SamplerBase.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::sampling {

class PlaneSampler : public SamplerBase::Register<PlaneSampler>
{
public:
    static std::string identifier() { return "PlaneSampler"; }
 
    explicit PlaneSampler(const CFDSim& /*unused*/);
 
    ~PlaneSampler() override;

    void initialize(const std::string& key) override;

    void check_bounds() override;

    void sampling_locations(SampleLocType& /*sample_locs*/) const override;

    void sampling_locations(
        SampleLocType& /*sample_locs*/,
        const amrex::Box& /*box*/) const override;
 
    void output_locations(SampleLocType& sample_locs) const override
    {
        sampling_locations(sample_locs);
    }
 
    void
    define_netcdf_metadata(const ncutils::NCGroup& /*unused*/) const override;
    void
    populate_netcdf_metadata(const ncutils::NCGroup& /*unused*/) const override;

    [[nodiscard]] std::string label() const override { return m_label; }
    std::string& label() override { return m_label; }

    [[nodiscard]] std::string sampletype() const override
    {
        return identifier();
    }

    [[nodiscard]] int id() const override { return m_id; }
    int& id() override { return m_id; }

    [[nodiscard]] long num_points() const override { return m_npts; }

    [[nodiscard]] long num_output_points() const override { return m_npts; }
 
private:
    const CFDSim& m_sim;
 
    amrex::Vector<amrex::Real> m_axis1;
    amrex::Vector<amrex::Real> m_axis2;
    amrex::Vector<amrex::Real> m_origin;
    amrex::Vector<amrex::Real> m_offset_vector{0.0_rt, 0.0_rt, 0.0_rt};
    amrex::Vector<amrex::Real> m_poffsets;
    amrex::Vector<int> m_npts_dir;
 
    std::string m_label;
 
    int m_id{-1};
    int m_npts{0};
};
 
} // namespace kynema_sgf::sampling
 
#endif /* PLANESAMPLER_H */