electrostatSolver.h

Go to the documentation of this file.

#ifndef electrostat_h
#define electrostat_h
 
#include <iostream>
#include <map>
#include "solver.h"
#include "config.h"
#include "algebra/algebra.h"
 
const int DIM_PB_ELEC = 1;
 
class electrostatSolver: public solver<DIM_PB_ELEC>
    {
public:
    electrostatSolver(
            Mesh::mesh & _msh ,
            std::vector<Tetra::prm> & _pTetra ,
            std::vector<Facette::prm> & _pFac ,
            const double _tol ,
            const bool v ,
            const int max_iter ):
            solver<DIM_PB_ELEC>(_msh,_pTetra,_pFac,"cg_dir",_tol,v,max_iter) {}
 
    std::vector<double> V;
 
    void infos(void);
 
    void integrales(Tetra::Tet const &tet, Eigen::Ref<Eigen::Matrix<double,Tetra::N,Tetra::N> > AE);
 
    void integrales(Facette::Fac const &fac, std::vector<double> &BE);
 
    bool save(const std::string V_fileName , std::string const &metadata ) const;
 
    double getSigma(Tetra::Tet const &tet) const;
 
    double getCurrentDensity(Facette::Fac const &fac) const;
 
    void compute(const bool verbose , const std::string V_fileName );
 
    void checkBoundaryConditions(void) const;
 
private:
    const int precision = 8;
 
    bool solve(void);
    };  // end class electrostatSolver
#endif