Line data    Source code

       1              : #ifndef AMREX_OPERATOR_DIAGONAL_H_
       2              : #define AMREX_OPERATOR_DIAGONAL_H_
       3              : 
       4              : #include <AMReX_MLCellLinOp.H>
       5              : #include <AMReX_Array.H>
       6              : #include <limits>
       7              : 
       8              : #include "Operator/Operator.H"
       9              : 
      10              : 
      11              : using namespace amrex;
      12              : 
      13              : namespace Operator
      14              : {
      15              : class Diagonal : public Operator<Grid::Node>
      16              : {
      17              : public:
      18              :     Diagonal () {}
      19              :     Diagonal (const Vector<Geometry>& a_geom,
      20              :         const Vector<BoxArray>& a_grids,
      21              :         const Vector<DistributionMapping>& a_dmap,
      22              :         const LPInfo& a_info);
      23              :     virtual ~Diagonal ();
      24              :     Diagonal (const Diagonal&) = delete;
      25              :     Diagonal (Diagonal&&) = delete;
      26              :     Diagonal& operator= (const Diagonal&) = delete;
      27              :     Diagonal& operator= (Diagonal&&) = delete;
      28              : protected:
      29              : 
      30              :     virtual void Fapply (int amrlev, int mglev, MultiFab& out, const MultiFab& in) const override final;
      31              :     virtual void Fsmooth (int amrlev, int mglev, MultiFab& sol, const MultiFab& rsh) const override final;
      32            0 :     virtual int getNComp() const override {return AMREX_SPACEDIM;};
      33            0 :     virtual void normalize (int /*amrlev*/, int /*mglev*/, MultiFab& /*mf*/) const override final {};
      34            0 :     virtual void reflux (int /*crse_amrlev*/,
      35              :                 MultiFab& /*res*/, const MultiFab& /*crse_sol*/, const MultiFab& /*crse_rhs*/,
      36            0 :                 MultiFab& /*fine_res*/, MultiFab& /*fine_sol*/, const MultiFab& /*fine_rhs*/) const override {};
      37              : 
      38              : };
      39              : }
      40              : #endif