AMReX_MLAlgMG.H

Go to the documentation of this file.

#ifndef AMREX_ML_ALGMG_H_
#define AMREX_ML_ALGMG_H_
#include <AMReX_Config.H>
 
#include <AMReX_Geometry.H>
#include <AMReX_LayoutData.H>
#include <AMReX_MultiFab.H>
#include <AMReX_iMultiFab.H>
 
#include <memory>
 
namespace amrex {
 
template <typename T> class AlgMG;
class MLNodeLinOp;
 
class MLAlgMG
{
public:
 
    MLAlgMG (int mglev, BoxArray const& grids, DistributionMapping const& dmap,
             Geometry const& geom, iMultiFab const& owner_mask,
             iMultiFab const& dirichlet_mask, MLNodeLinOp const& linop);
 
    MLAlgMG (int mglev, BoxArray const& grids, DistributionMapping const& dmap,
             Geometry const& geom, FabFactory<FArrayBox> const& factory,
             iMultiFab const* overset_mask, Real ascalar, Real bscalar,
             MultiFab const* acoef, Array<MultiFab const*,AMREX_SPACEDIM> const& bcoef,
             MultiFab const* eb_bcoef,
             LayoutData<GpuArray<int,2*AMREX_SPACEDIM>> const& bctype,
             LayoutData<GpuArray<Real,2*AMREX_SPACEDIM>> const& bcl, int maxorder);
 
    ~MLAlgMG ();
    MLAlgMG (MLAlgMG const&) = delete;
    MLAlgMG (MLAlgMG&&) = delete;
    MLAlgMG& operator= (MLAlgMG const&) = delete;
    MLAlgMG& operator= (MLAlgMG&&) = delete;
 
    void solve (MultiFab& soln, MultiFab const& rhs, Real reltol, Real abstol, int maxiter);
 
    void applyVcycle (MultiFab& soln, MultiFab const& rhs);
 
    [[nodiscard]] AlgMG<Real>& solver () noexcept;
 
    struct Impl;
 
private:
    std::unique_ptr<Impl> m_impl;
};
 
}
 
#endif