jacsolve.h

Go to the documentation of this file.

#ifndef JACSOLVE_FILE
#define JACSOLVE_FILE
#include "geom.h"
#include "getmatrix.h"
#include <vector>
 
 
void JacobiSolve(CRS_Matrix const &SK, std::vector<double> const &f, std::vector<double> &u);
 
 
void JacobiSmoother(Matrix const &SK, std::vector<double> const &f, std::vector<double> &u,
                    std::vector<double> & r, int nsmooth=1, double const omega=1.0, bool zero=false);
 
void DiagPrecond(Matrix const &SK, std::vector<double> const &r, std::vector<double> &w,
                 double const omega=1.0);
 
 
 
class Multigrid
{
    public:
       Multigrid(Mesh const& cmesh, int nlevel);
 
       Multigrid(Multigrid const&)            = delete;
       Multigrid& operator=(Multigrid const&) = delete;
 
       ~Multigrid();
 
       size_t Nlevels() const
       {return _meshes.size(); }
 
       int Ndofs() const
       {return _meshes[Nlevels()-1].Nnodes(); }
 
       Mesh const& GetMesh(int lev) const
       { return _meshes[lev]; }
 
       std::vector<double> const&  GetSolution(int lev) const
       { return _u.at(lev); }
 
       void DefineOperators();
 
       void DefineOperator(size_t lev);
 
       void JacobiSolve(size_t lev);
 
       void MG_Step(size_t lev, int pre_smooth=1, bool const bzero=false, int nu=1);
 
       void MG_Solve(int pre_smooth=1, double eps=1e-6, int nu=1);
 
 
    private:
       gMesh_Hierarchy _meshes;                
       std::vector<FEM_Matrix>           _SK;  
       std::vector<std::vector<double>>  _u;   
       std::vector<std::vector<double>>  _f;   
       std::vector<std::vector<double>>  _d;   
       std::vector<std::vector<double>>  _w;   
       std::vector<BisectIntDirichlet> _Pc2f;  
 
};
 
 
 
 
#endif