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 "solver.h"
#include "meshUtils.h"
#include "node.h"
#include "tetra.h"
 
#include "algebra/algebra.h"
 
const int DIM_PB_MAG = 2;
 
class LinAlgebra : public solver<DIM_PB_MAG>
    {
public:
    LinAlgebra(Settings &s , Mesh::mesh &my_msh )
        : solver<DIM_PB_MAG>(
            my_msh, s.paramTetra, s.paramFacette, "bicg_dir", s.TOL, s.verbose, s.MAXITER,
            [this](Mesh::Edge edge){ return msh->magNode[edge.first] && msh->magNode[edge.second]; }
        ),
        verbose(s.verbose),
        extSpaceField(my_msh.extSpaceField)
        {
        Xw.resize(2*NOD);
        base_projection();
        if (!s.recenter)
            { idx_dir = Nodes::IDX_UNDEF; }
        else
            { idx_dir = s.recentering_direction; }
 
        for (int i=0;i<NOD;i++)
            {
            if(!my_msh.magNode[i])
                {
                lvd.push_back(2*i);
                lvd.push_back(2*i+1);
                }
            }
        lvd.shrink_to_fit();
        }
 
    void checkBoundaryConditions(void) const {}
 
    void prepareElements(Eigen::Vector3d const &Hext , timing const &t_prm );
 
    void prepareElements(double const A_Hext , timing const &t_prm );
 
    void buildInitGuess(std::vector<double> &G) const;
 
    bool solve(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;
 
    std::vector<double> Xw;
 
    const int verbose;
 
    double DW_vz;
 
    double v_max;
 
    const std::vector< Eigen::Matrix<double,Nodes::DIM,Tetra::NPI> >& extSpaceField;
 
    std::vector<int> lvd;
    };// end class linAlgebra
#endif