GeometryRefinement.H

Go to the documentation of this file.

#ifndef GEOMETRYREFINEMENT_H
#define GEOMETRYREFINEMENT_H
 
#include "src/utilities/index_operations.H"
#include "src/utilities/tagging/RefinementCriteria.H"
 
namespace kynema_sgf {
namespace tagging {

class GeometryType
    : public Factory<GeometryType, const CFDSim&, const std::string&>
{
public:
    static std::string base_identifier() { return "GeometryType"; }
 
    GeometryType() = default;
 
    ~GeometryType() override = default;
 
    virtual void operator()(
        const amrex::Box&,
        const amrex::Geometry& geom,
        const amrex::Array4<amrex::TagBox::TagType>& tags,
        const tagging::TaggingOperator& op) const = 0;
 
    [[nodiscard]] virtual const amrex::RealBox& bound_box() const = 0;
};
 
} // namespace tagging

class GeometryRefinement
    : public RefinementCriteria::Register<GeometryRefinement>
{
public:
    static std::string identifier() { return "GeometryRefinement"; }
 
    explicit GeometryRefinement(const CFDSim& /*sim*/);
 
    ~GeometryRefinement() override = default;

    void initialize(const std::string& key) override;
 
    void
    operator()(int level, amrex::TagBoxArray& tags, amrex::Real time, int ngrow)
        override;
 
private:
    const CFDSim& m_sim;
 
    amrex::Vector<std::unique_ptr<tagging::GeometryType>> m_geom_refiners;

    int m_set_level{-1};
 
    // If the user hasn't specified a particular level, act on a range of levels
    int m_min_level{0};
    int m_max_level{32};
};
 
} // namespace kynema_sgf
 
#endif /* GEOMETRYREFINEMENT_H */