MLPoissonSolver.h

Go to the documentation of this file.

#ifndef ML_POISSON_SOLVER_H_
#define ML_POISSON_SOLVER_H_

#include "Solvers/AmrPoissonSolver.h"
#include "Amr/AmrBoxLib.h"

class MLPoissonSolver : public AmrPoissonSolver< AmrBoxLib > {
    
private:
    typedef AmrBoxLib::AmrGeomContainer_t GeomContainer_t;
    typedef amrex::Vector<AmrBoxLib::AmrField_t*> AmrFieldContainer_pt;
    typedef amrex::Vector<const AmrBoxLib::AmrField_t*> const_AmrFieldContainer_pt;
    typedef AmrBoxLib::AmrGeometry_t AmrGeometry_t;
    typedef AmrBoxLib::AmrGrid_t AmrGrid_t;
    typedef AmrBoxLib::AmrProcMap_t           AmrProcMap_t;
    typedef AmrBoxLib::AmrScalarFieldContainer_t    AmrScalarFieldContainer_t;
    typedef AmrBoxLib::AmrVectorFieldContainer_t    AmrVectorFieldContainer_t;
    typedef AmrBoxLib::AmrGridContainer_t AmrGridContainer_t;
    typedef AmrBoxLib::AmrProcMapContainer_t AmrProcMapContainer_t;
    
public:
    
    MLPoissonSolver(AmrBoxLib* itsAmrObject_p);
    
    void solve(AmrScalarFieldContainer_t& rho,
               AmrScalarFieldContainer_t& phi,
               AmrVectorFieldContainer_t& efield,
               unsigned short baseLevel,
               unsigned short finestLevel,
               bool prevAsGuess = true);
    
    double getXRangeMin(unsigned short level = 0);
    double getXRangeMax(unsigned short level = 0);
    double getYRangeMin(unsigned short level = 0);
    double getYRangeMax(unsigned short level = 0);
    double getZRangeMin(unsigned short level = 0);
    double getZRangeMax(unsigned short level = 0);
    
    
    Inform &print(Inform &os) const;
    
    
private:
    
    void mlmg_m(AmrScalarFieldContainer_t& rho,
                AmrScalarFieldContainer_t& phi,
                AmrVectorFieldContainer_t &efield,
                int baseLevel,
                int finestLevel);
    
private:
    double reltol_m;                
    double abstol_m;                
};


inline Inform &operator<<(Inform &os, const MLPoissonSolver &fs) {
    return fs.print(os);
}

#endif