Diagonal.H

Go to the documentation of this file.

#ifndef AMREX_OPERATOR_DIAGONAL_H_
#define AMREX_OPERATOR_DIAGONAL_H_
 
#include <AMReX_MLCellLinOp.H>
#include <AMReX_Array.H>
#include <limits>
 
#include "Operator/Operator.H"
 
 
using namespace amrex;
 
namespace Operator
{
class Diagonal : public Operator<Grid::Node>
{
public:
    Diagonal () {}
    Diagonal (const Vector<Geometry>& a_geom,
        const Vector<BoxArray>& a_grids,
        const Vector<DistributionMapping>& a_dmap,
        const LPInfo& a_info);
    virtual ~Diagonal ();
    Diagonal (const Diagonal&) = delete;
    Diagonal (Diagonal&&) = delete;
    Diagonal& operator= (const Diagonal&) = delete;
    Diagonal& operator= (Diagonal&&) = delete;
protected:
 
    virtual void Fapply (int amrlev, int mglev, MultiFab& out, const MultiFab& in) const override final;
    virtual void Fsmooth (int amrlev, int mglev, MultiFab& sol, const MultiFab& rsh) const override final;
    virtual int getNComp() const override {return AMREX_SPACEDIM;};
    virtual void normalize (int /*amrlev*/, int /*mglev*/, MultiFab& /*mf*/) const override final {};
    virtual void reflux (int /*crse_amrlev*/,
                MultiFab& /*res*/, const MultiFab& /*crse_sol*/, const MultiFab& /*crse_rhs*/,
                MultiFab& /*fine_res*/, MultiFab& /*fine_sol*/, const MultiFab& /*fine_rhs*/) const override {};
 
};
}
#endif