spinAccumulationSolver.h

#ifndef spinAccumulationSolver_h
#define spinAccumulationSolver_h
 
#include <vector>
#include "config.h"
#include "node.h"
#include "tetra.h"
#include "electrostatSolver.h"
#include "solver.h"
#include "meshUtils.h"
 
const int DIM_PB_SPIN_ACC = 3;
 
class spinAcc : public solver<DIM_PB_SPIN_ACC>
    {
    public:
    spinAcc(const Settings &mySettings ,
            Mesh::mesh &_msh ,
            const double _tol ,
            const int max_iter );
 
    void boundaryConditions(void); // should be private
 
    bool compute(const Nodes::step k = Nodes::NEXT);
 
    double jn_multiplier = 1.0;
 
    std::vector<double> potential(const double factor) const;
 
    std::vector<Eigen::Vector3d> s;
 
    void checkBoundaryConditions(void) const override;
 
    private:
    std::vector<double> valDirichlet;
 
    std::vector<int> idxDirichlet;
 
    void fillDirichletData(const int k, Eigen::Vector3d &s_value);
 
    std::vector<double> V;
 
    double V_D = 0.0;
 
    const int precision = 8;
 
    double getMs(const Tetra::Tet &tet) const;
 
    double getSigma(const Tetra::Tet &tet) const;
 
    double getPolarizationRate(const Tetra::Tet &tet) const;
 
    double getDiffusionCst(const Tetra::Tet &tet) const;
 
    double getLsd(const Tetra::Tet &tet) const;
 
    double getLsf(const Tetra::Tet &tet) const;
 
    Eigen::Matrix<double,Nodes::DIM,Tetra::N> calc_u_nod(const Tetra::Tet &tet) const;
 
    Nodes::step uStep = Nodes::NEXT;
 
    void prepareExtraField(void) const;
 
    bool solve(void);
 
    void integrales(const Tetra::Tet &tet, Eigen::Matrix<double,DIM_PB*Tetra::N,DIM_PB*Tetra::N> &AE) const;
 
    void integrales(const Tetra::Tet &tet ,
                    std::vector<double> &BE ) const;
    };
 
#endif