spinAcc Class Reference

#include <spinAccumulationSolver.h>

Inheritance diagram for spinAcc:

Public Member Functions
	spinAcc (const Settings &mySettings, Mesh::mesh &_msh, const double _tol, const int max_iter)
void	boundaryConditions (void)
bool	compute (void)
void	checkBoundaryConditions (void) const override
Public Member Functions inherited from solver< DIM_PB_SPIN_ACC >
	solver (Mesh::mesh &_msh, const std::vector< Tetra::prm > &_pTetra, const std::vector< Triangle::prm > &_pTri, const std::string &name, const double _tol, const bool v, const int max_iter, const std::function< bool(Mesh::Edge)> &edge_filter=[](Mesh::Edge){ return true;})

Public Attributes
std::vector< Eigen::Vector3d >	s

Private Member Functions
void	fillDirichletData (const int k, Eigen::Vector3d &s_value)
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
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 (Tetra::Tet &tet, std::vector< double > &BE)

Private Attributes
std::vector< double >	valDirichlet
std::vector< int >	idxDirichlet
std::vector< double >	V
const int	precision = 8

Additional Inherited Members
Protected Member Functions inherited from solver< DIM_PB_SPIN_ACC >
algebra::MatrixShape	build_shape (const std::function< bool(Mesh::Edge)> &edge_filter) const
void	buildMat (std::array< int, N > &ind, Eigen::Matrix< double, DIM_PROBLEM *N, DIM_PROBLEM *N > &Ke)
void	buildVect (std::array< int, N > &ind, std::vector< double > &Le)
Protected Attributes inherited from solver< DIM_PB_SPIN_ACC >
Mesh::mesh *	msh
const int	NOD
const std::vector< Tetra::prm > &	paramTet
const std::vector< Triangle::prm > &	paramTri
const bool	verbose
algebra::iteration< double >	iter
algebra::SparseMatrix	K
std::vector< double >	L_rhs
Static Protected Attributes inherited from solver< DIM_PB_SPIN_ACC >
static const int	DIM_PB

Detailed Description

container for Spin Accumulation constants, solver and related datas. The model obeys the diffusive spin equation. The boundary conditions are Dirichlet type. User has to provide through mesh and settings a surface S0 where spin diffusion vector s is a constant (zero recommended), and another surface S1 where current density and spin polarization is defined. This surface S1 must be the same as the one given to the potential solver for its own boundary conditions.

Constructor & Destructor Documentation

spinAcc()

spinAcc::spinAcc	(	const Settings &	mySettings,
		Mesh::mesh &	_msh,
		const double	_tol,
		const int	max_iter
	)

spin accumulation constructor if mySettings.spin_acc bool is false it is a do nothing constructor

Parameters

[in]	mySettings
[in]	_msh	ref to the mesh
[in]	_tol	tolerance for bicg_dir solver
[in]	max_iter	maximum number of iterations

Member Function Documentation

boundaryConditions()

void spinAcc::boundaryConditions

(

void

)

boundary conditions: a surface with a fixed s, and another surface with fixed normal current density J and polarization vector P set valDirichlet values and fill vector of indices idxDirichlet

checkBoundaryConditions()

void spinAcc::checkBoundaryConditions ( void  ) const

overridevirtual

check boundary conditions: mesh and settings have to define a single surface with constant normal current density J, a vector polarization P and another single surface where spin diffusion = 0

Implements solver< DIM_PB_SPIN_ACC >.

compute()

bool spinAcc::compute

(

void

)

call solver and update spin diffusion solution, returns true if solver succeeded

fillDirichletData()

void spinAcc::fillDirichletData ( const int  k, Eigen::Vector3d &  s_value  )

private

fill valDirichlet and idxDirichlet vectors with k, a node index from tri.ind and the corresponding spin diffusion s value

getDiffusionCst()

double spinAcc::getDiffusionCst ( const Tetra::Tet &  tet ) const

private

diffusion constant, units: s^-1 m^2

getLsd()

double spinAcc::getLsd ( const Tetra::Tet &  tet ) const

private

length s-d : only in magnetic material

getLsf()

double spinAcc::getLsf ( const Tetra::Tet &  tet ) const

private

spin flip length : exists in both non-magnetic and magnetic metals

getMs()

double spinAcc::getMs ( const Tetra::Tet &  tet ) const

private

returns Ms

getPolarizationRate()

double spinAcc::getPolarizationRate ( const Tetra::Tet &  tet ) const

private

\( P \) is polarization rate of the current density

getSigma()

double spinAcc::getSigma ( const Tetra::Tet &  tet ) const

private

returns sigma of the tetraedron, (conductivity in (Ohm.m)^-1

integrales() [1/2]

void spinAcc::integrales ( const Tetra::Tet &  tet, Eigen::Matrix< double, DIM_PB *Tetra::N, DIM_PB *Tetra::N > &  AE  ) const

private

computes all contributions to matrix AE from tetrahedron tet (LHS) all = non magnetic metal + magnetic metal

integrales() [2/2]

void spinAcc::integrales ( Tetra::Tet &  tet, std::vector< double > &  BE  )

private

computes magnetic metal contributions to spin diffusion from tetrahedron tet (RHS)

prepareExtraField()

void spinAcc::prepareExtraField ( void  ) const

private

affect extraField member function of all tetrahedrons extraField is computing the contribution from the spin diffusion s to llg

solve()

bool spinAcc::solve ( void  )

private

solver, using biconjugate stabilized gradient, with diagonal preconditionner and Dirichlet boundary conditions

Member Data Documentation

idxDirichlet

std::vector<int> spinAcc::idxDirichlet

private

list of the indices for Dirichlet boundary conditions, it contains the indices of the nodes where s is given by the user to be a constant on a surface

precision

const int spinAcc::precision = 8

private

number of digits in the optional output file

s

std::vector<Eigen::Vector3d> spinAcc::s

solution of the spin diffusion vector s over the nodes

V

std::vector<double> spinAcc::V

private

electrostatic potential V on the nodes unit [V] = Volt = kg A^-1 s^-3

valDirichlet

std::vector<double> spinAcc::valDirichlet

private

Dirichlet values of the components of s on the nodes, it is zero if the node is not in idxDirichlet

---

The documentation for this class was generated from the following files:

• src/spinAccumulationSolver.h
• src/spinAccumulationSolver.cpp