electrostatSolver.h

Go to the documentation of this file.

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