AMReX_MFInterpolater.H

Go to the documentation of this file.

#ifndef AMREX_MF_INTERPOLATER_H_
#define AMREX_MF_INTERPOLATER_H_
#include <AMReX_Config.H>
 
#include <AMReX_Extension.H>
#include <AMReX_InterpBase.H>
 
namespace amrex {
 
class MultiFab;
class Geometry;
 
class MFInterpolater
    : public InterpBase
{
public:
    virtual void interp (MultiFab const& crsemf, int ccomp, MultiFab& finemf, int fcomp, int ncomp,
                         IntVect const& ng, Geometry const& cgeom, Geometry const& fgeom,
                         Box const& dest_domain, IntVect const& ratio,
                         Vector<BCRec> const& bcs, int bcscomp) = 0;
};
 
class MFPCInterp final
    : public MFInterpolater
{
public:
    Box CoarseBox (Box const& fine, int ratio) override;
    Box CoarseBox (Box const& fine, IntVect const& ratio) override;
 
    void interp (MultiFab const& crsemf, int ccomp, MultiFab& finemf, int fcomp, int ncomp,
                         IntVect const& ng, Geometry const& cgeom, Geometry const& fgeom,
                         Box const& dest_domain, IntVect const& ratio,
                         Vector<BCRec> const& bcs, int bcscomp) override;
};
 
class MFCellConsLinInterp
    : public MFInterpolater
{
public:
    explicit MFCellConsLinInterp (bool do_linear_limiting_)
        : do_linear_limiting(do_linear_limiting_) {}
 
    Box CoarseBox (Box const& fine, int ratio) override;
    Box CoarseBox (Box const& fine, IntVect const& ratio) override;
 
    void interp (MultiFab const& crsemf, int ccomp, MultiFab& finemf, int fcomp, int ncomp,
                         IntVect const& ng, Geometry const& cgeom, Geometry const& fgeom,
                         Box const& dest_domain, IntVect const& ratio,
                         Vector<BCRec> const& bcs, int bcomp) override;
protected:
    bool do_linear_limiting = true;
};
 
class MFCellConsLinMinmaxLimitInterp
    : public MFInterpolater
{
public:
    Box CoarseBox (Box const& fine, int ratio) override;
    Box CoarseBox (Box const& fine, IntVect const& ratio) override;
 
    void interp (MultiFab const& crsemf, int ccomp, MultiFab& finemf, int fcomp, int ncomp,
                         IntVect const& ng, Geometry const& cgeom, Geometry const& fgeom,
                         Box const& dest_domain, IntVect const& ratio,
                         Vector<BCRec> const& bcs, int bcomp) override;
};
 
class MFCellBilinear final
    : public MFInterpolater
{
public:
    Box CoarseBox (Box const& fine, int ratio) override;
    Box CoarseBox (Box const& fine, IntVect const& ratio) override;
 
    void interp (MultiFab const& crsemf, int ccomp, MultiFab& finemf, int fcomp, int ncomp,
                         IntVect const& ng, Geometry const& cgeom, Geometry const& fgeom,
                         Box const& dest_domain, IntVect const& ratio,
                         Vector<BCRec> const& bcs, int bcscomp) override;
};
 
/*
 * \brief [Bi|Tri] linear interpolation on nodal data
 */
class MFNodeBilinear final
    : public MFInterpolater
{
public:
    Box CoarseBox (Box const& fine, int ratio) override;
    Box CoarseBox (Box const& fine, IntVect const& ratio) override;
 
    void interp (MultiFab const& crsemf, int ccomp, MultiFab& finemf, int fcomp, int ncomp,
                         IntVect const& ng, Geometry const& cgeom, Geometry const& fgeom,
                         Box const& dest_domain, IntVect const& ratio,
                         Vector<BCRec> const& bcs, int bcscomp) override;
};
 
extern AMREX_EXPORT MFPCInterp          mf_pc_interp;
extern AMREX_EXPORT MFCellConsLinInterp mf_cell_cons_interp;
extern AMREX_EXPORT MFCellConsLinInterp mf_lincc_interp;
extern AMREX_EXPORT MFCellConsLinMinmaxLimitInterp mf_linear_slope_minmax_interp;
extern AMREX_EXPORT MFCellBilinear      mf_cell_bilinear_interp;
extern AMREX_EXPORT MFNodeBilinear      mf_node_bilinear_interp;
 
}
 
#endif