linear_algebra.h

Go to the documentation of this file.

#ifndef linear_algebra_h
#define linear_algebra_h
 
#include <random>
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wmaybe-uninitialized"
#include <execution>
#pragma GCC diagnostic pop
 
#include "config.h"
 
#include "facette.h"
#include "feellgoodSettings.h"
#include "mesh.h"
#include "node.h"
#include "tetra.h"
 
class LinAlgebra
    {
public:
    inline LinAlgebra(Settings &s , Mesh::mesh &my_msh )
        : NOD(my_msh.getNbNodes()), MAXITER(s.MAXITER), TOL(s.TOL), ILU_tol(s.ILU_tol),
          ILU_fill_factor(s.ILU_fill_factor), verbose(s.verbose),
          prmTetra(s.paramTetra), prmFacette(s.paramFacette), refMsh(&my_msh)
        {
        Eigen::setNbThreads(s.solverNbTh);
        base_projection();
        if (!s.recenter)
            { idx_dir = Nodes::IDX_UNDEF; }
        else
            { idx_dir = s.recentering_direction; }
 
        if(s.getFieldType() == R4toR3)
            { setExtSpaceField(s); }
        }
 
    void prepareElements(Eigen::Vector3d const &Hext , timing const &t_prm );
 
    void prepareElements(double const A_Hext , timing const &t_prm );
 
    void buildInitGuess(Eigen::Ref<Eigen::VectorXd> G) const;
 
    int solver(timing const &t_prm );
 
    inline void set_DW_vz(double vz ) { DW_vz = vz; }
 
    inline double get_v_max(void) { return v_max; }
 
    void setExtSpaceField(Settings &s );
 
    void base_projection();
private:
    Nodes::index idx_dir;
 
    const int NOD;
 
    const int MAXITER;
 
    const double TOL;
 
    double ILU_tol;
 
    double ILU_fill_factor;
 
    const int verbose;
 
    const std::vector<Tetra::prm> &prmTetra;
 
    const std::vector<Facette::prm> &prmFacette;
 
    Mesh::mesh *refMsh;
 
    double DW_vz;
 
    double v_max;
 
    std::vector< Eigen::Matrix<double,Nodes::DIM,Tetra::NPI> > extSpaceField;
    };// end class linAlgebra
#endif