#include <mesh.h>

Public Member Functions
	mesh (Settings &mySets)

	mesh (const mesh &)=delete

mesh &	operator= (const mesh &)=delete

int	getNbNodes (void) const

int	getNbFacs (void) const

int	getNbTets (void) const

const Eigen::Vector3d	getNode_p (const int i) const

const Eigen::Vector3d	getNode_u (const int i) const

const Eigen::Vector3d	getNode_v (const int i) const

double	getProj_ep (const int i) const

double	getProj_eq (const int i) const

void	set_node_u0 (const int i, Eigen::Vector3d const &val)

void	set_node_zero_v (const int i)

void	infos (void) const

void	setBasis (const double r)

void	updateNode (int i, double vp, double vq, const double dt)

void	evolution (void)

double	readSol (bool VERBOSE, const std::string fileName)

void	init_distrib (Settings const &mySets)

double	avg (std::function< double(Nodes::Node, Nodes::index)> getter, Nodes::index d, int region=-1) const

double	max_angle () const

void	savesol (const int precision, const std::string fileName, std::string const &metadata, bool withSpinAcc, std::vector< Eigen::Vector3d > &s) const

void	set (const int i, std::function< void(Nodes::Node &, const double)> what_to_set, const double val)

int	getNodeIndex (const int i) const

bool	isMagnetic (const Tetra::Tet &t)

bool	isMagnetic (const Facette::Fac &f)

Public Attributes
Eigen::Vector3d	c

Eigen::Vector3d	l

double	vol

double	totalMagVol

std::vector< Facette::Fac >	fac

std::vector< Tetra::Tet >	tet

std::vector< Tetra::prm > &	paramTetra

std::vector< std::pair< int, int > >	edges

std::vector< bool >	magNode

std::vector< int >	magTet

std::vector< int >	magFac

Private Member Functions
void	checkMeshFile (Settings const &mySets)

void	readNodes (Settings const &mySets)

void	readTetraedrons (Settings const &mySets)

void	readTriangles (Settings const &mySets)

void	readMesh (Settings const &mySets)

double	doOnNodes (const double init_val, const Nodes::index coord, std::function< bool(double, double)> whatToDo) const

double	minNodes (const Nodes::index coord) const

double	maxNodes (const Nodes::index coord) const

void	indexReorder ()

void	sortNodes (Nodes::index long_axis)

double	surface (std::vector< int > &facIndices)

bool	isInMagList (std::vector< int > &idxMagList, Facette::Fac &f)

Private Attributes
std::vector< Nodes::Node >	node

std::vector< int >	node_index

std::vector< std::vector< int > >	volumeRegions

Detailed Description

class for storing the mesh, including mesh geometry values, containers for the nodes, triangular faces and tetrahedrons. nodes data are private. They are accessible only through getter and setter.

Constructor & Destructor Documentation

◆ mesh() [1/2]

Mesh::mesh::mesh ( Settings & mySets )

inline

constructor : read mesh file, reorder indices and computes some values related to the mesh : center and length along coordinates,full volume

Parameters

[in] mySets

◆ mesh() [2/2]

Mesh::mesh::mesh ( const mesh & )

delete

operator copy deleted

Member Function Documentation

◆ avg()

double mesh::avg	(	std::function< double(Nodes::Node, Nodes::index)>	getter,
		Nodes::index	d,
		int	region = `-1`
	)		const

average component of either u or v through getter on the whole set of tetetrahedron

Parameters

[in]	getter
[in]	d
	region	region index, or -1 for all magnetic regions

◆ checkMeshFile()

void Mesh::mesh::checkMeshFile ( Settings const & mySets )

private

test if mesh file contains surfaces and regions mentionned in yaml settings and their dimensions

Parameters

[in] mySets

◆ doOnNodes()

double mesh::doOnNodes	(	const double	init_val,
		const Nodes::index	coord,
		std::function< bool(double, double)>	whatToDo
	)		const

private

loop on nodes to apply predicate 'whatTodo'

Parameters

[in]	init_val
[in]	coord
[in]	whatToDo

◆ evolution()

void Mesh::mesh::evolution ( void )

inline

call evolution for all the nodes

◆ getNbFacs()

int Mesh::mesh::getNbFacs ( void ) const

inline

return number of triangular fac

◆ getNbNodes()

int Mesh::mesh::getNbNodes ( void ) const

inline

return number of nodes

◆ getNbTets()

int Mesh::mesh::getNbTets ( void ) const

inline

return number of tetrahedrons

◆ getNode_p()

const Eigen::Vector3d Mesh::mesh::getNode_p ( const int i ) const

inline

getter : return node.p

◆ getNode_u()

const Eigen::Vector3d Mesh::mesh::getNode_u ( const int i ) const

inline

getter : return node.u

◆ getNode_v()

const Eigen::Vector3d Mesh::mesh::getNode_v ( const int i ) const

inline

getter : return node.v

◆ getNodeIndex()

int Mesh::mesh::getNodeIndex ( const int i ) const

inline

getter for the node_index value at position i

◆ getProj_ep()

double Mesh::mesh::getProj_ep ( const int i ) const

inline

return projection of speed at node i along ep

◆ getProj_eq()

double Mesh::mesh::getProj_eq ( const int i ) const

inline

return projection of speed at node i along eq

◆ indexReorder()

void mesh::indexReorder ( )

private

redefine orientation of triangular faces in accordance with the tetrahedron reorientation of the tetrahedrons if needed; definition of Ms on facette elements Indices and orientation convention :

                v
              .
            ,/
           /
        2(ic)                                 2
      ,/|`\                                 ,/|`\
    ,/  |  `\                             ,/  |  `\
  ,/    '.   `\                         ,6    '.   `5
,/       |     `\                     ,/       8     `\

,/ | \ ,/ |\ 0(ia)----—'.-----—1(ib) --> u 0-----—4–'.-----—1 \. | ,/. | ,/ \. | ,/. | ,9 ‘. ’. ,/ ‘7. ’. ,/ \. |/. |/ 3(id)3 \. w

◆ infos()

void mesh::infos ( void ) const

basic informations on the mesh

◆ init_distrib()

void Mesh::mesh::init_distrib ( Settings const & mySets )

inline

computes an analytical initial magnetization distribution as a starting point for the simulation. If the node is not magnetic then it is set to NAN.

Parameters

[in] mySets

◆ isInMagList()

bool Mesh::mesh::isInMagList	(	std::vector< int > &	idxMagList,
		Facette::Fac &	f
	)

inlineprivate

return true if facette f is already indexed in the list idxMagList. Uses operator== for Fac.

◆ isMagnetic() [1/2]

bool Mesh::mesh::isMagnetic ( const Facette::Fac & f )

inline

return true if this facet is magnetic

◆ isMagnetic() [2/2]

bool Mesh::mesh::isMagnetic ( const Tetra::Tet & t )

inline

return true if this tetraedron is magnetic

◆ max_angle()

double Mesh::mesh::max_angle ( ) const

inline

Compute the maximum angle of the magnetization between two adjacent nodes.

◆ maxNodes()

double Mesh::mesh::maxNodes ( const Nodes::index coord ) const

inlineprivate

return the maximum of all nodes coordinate along coord axis

Parameters

[in] coord

◆ minNodes()

double Mesh::mesh::minNodes ( const Nodes::index coord ) const

inlineprivate

return the minimum of all nodes coordinate along coord axis

Parameters

[in] coord

◆ operator=()

mesh& Mesh::mesh::operator= ( const mesh & )

delete

operator= deleted

◆ readMesh()

void Mesh::mesh::readMesh ( Settings const & mySets )

private

reading mesh format 2.2 text file function

Parameters

[in] mySets

◆ readNodes()

void Mesh::mesh::readNodes ( Settings const & mySets )

private

read Nodes from mesh file

Parameters

[in] mySets

◆ readSol()

double Mesh::mesh::readSol	(	bool	VERBOSE,
		const std::string	fileName
	)

read a solution from a file (tsv formated) and initialize fem struct to restart computation from that distribution, return time

Parameters

[in]	VERBOSE
[in]	fileName	input .sol text file

◆ readTetraedrons()

void Mesh::mesh::readTetraedrons ( Settings const & mySets )

private

read tetraedrons of the settings volume regions

Parameters

[in] mySets

◆ readTriangles()

void Mesh::mesh::readTriangles ( Settings const & mySets )

private

read facettes of the settings surface regions

Parameters

[in] mySets

◆ savesol()

void Mesh::mesh::savesol	(	const int	precision,
		const std::string	fileName,
		std::string const &	metadata,
		bool	withSpinAcc,
		std::vector< Eigen::Vector3d > &	s
	)		const

text file (tsv) writing function for a solution, node indices are zero based first column is index, then magnetization components, then scalar magnetic potential. If spin accumulation was involved, there are three extra columns for the componants of the spin diffusion vector. Since some volume region might be non magnetic, magnetization is undefined on those nodes and nan is used

Parameters

[in]	precision	numeric precision in .sol output text file
[in]	fileName
[in]	metadata
[in]	withSpinAcc
[in]	s	spin accumulation (might be empty)

◆ set()

void Mesh::mesh::set	(	const int	i,
		std::function< void(Nodes::Node &, const double)>	what_to_set,
		const double	val
	)

inline

setter for node[i]; what_to_set will fix what is the part of the node struct to set (usefull for fmm_demag.h)

Parameters

[in]	i
[in]	what_to_set
[in]	val

◆ set_node_u0()

void Mesh::mesh::set_node_u0	(	const int	i,
		Eigen::Vector3d const &	val
	)

inline

setter for u0

◆ set_node_zero_v()

void Mesh::mesh::set_node_zero_v ( const int i )

inline

fix to zero node[i].v

◆ setBasis()

void Mesh::mesh::setBasis ( const double r )

inline

call setBasis for all nodes, and update P matrix for all elements

◆ sortNodes()

void mesh::sortNodes ( Nodes::index long_axis )

private

Sort the nodes along the longest axis of the sample. This should reduce the bandwidth of the matrix we will have to solve for.

Parameters

[in] long_axis

◆ surface()

double mesh::surface ( std::vector< int > & facIndices )

private

returns the surface defined by the set of facets of indices in facIndices each elementary surface triangle defined by points p0,p1,p2 is computed using norm(cross(p0p1,p0p2))/2, it is always positive

◆ updateNode()

void Mesh::mesh::updateNode	(	int	i,
		double	vp,
		double	vq,
		const double	dt
	)

inline

make_evol on i^th node

Member Data Documentation

◆ c

Eigen::Vector3d Mesh::mesh::c

isobarycenter

◆ edges

std::vector<std::pair<int, int> > Mesh::mesh::edges

Edges of all the tetrahedrons, i.e. all the unique pairs of adjacent node indices. Each pair is sorted: first < second. The list is sorted lexicographically, as per std::pair::operator<().

◆ fac

std::vector<Facette::Fac> Mesh::mesh::fac

face container

◆ l

Eigen::Vector3d Mesh::mesh::l

lengths along x,y,z axis

◆ magFac

std::vector<int> Mesh::mesh::magFac

list of the indices of all magnetic facets from all surface regions

◆ magNode

std::vector<bool> Mesh::mesh::magNode

list of the magnetic nodes, using inner indices (non gmsh indices). If true it is magnetic. a node is magnetic if it belongs to a magnetic tetrahedron. Consequently any node on a magnetic/non magnetic interface is set to true in magNode.

◆ magTet

std::vector<int> Mesh::mesh::magTet

list of the indices of all magnetic tetrahedrons from all volume regions

◆ node

std::vector<Nodes::Node> Mesh::mesh::node

private

node container: not initialized by constructor, but later while reading the mesh by member function init_node

◆ node_index

std::vector<int> Mesh::mesh::node_index

private

Index of a node in the node vector. This vector is itself indexed by the node position in the *.msh and *.sol files. In other words, the node found at file_idx in a file is stored as node[node_index[file_idx]].

This is the inverse of the permutation we applied when sorting the nodes.

◆ paramTetra

std::vector<Tetra::prm>& Mesh::mesh::paramTetra

Reference to the volume regions in Settings.

◆ tet

std::vector<Tetra::Tet> Mesh::mesh::tet

tetrahedron container

◆ totalMagVol

double Mesh::mesh::totalMagVol

total magnetic volume of the mesh

◆ vol

double Mesh::mesh::vol

total volume of the mesh

◆ volumeRegions

std::vector<std::vector<int> > Mesh::mesh::volumeRegions

private

List of tetrahedrons making each volume region.

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

mesh.h
mesh.cpp
read.cpp
save.cpp

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ mesh() [1/2]

◆ mesh() [2/2]

Member Function Documentation

◆ avg()

◆ checkMeshFile()

◆ doOnNodes()

◆ evolution()

◆ getNbFacs()

◆ getNbNodes()

◆ getNbTets()

◆ getNode_p()

◆ getNode_u()

◆ getNode_v()

◆ getNodeIndex()

◆ getProj_ep()

◆ getProj_eq()

◆ indexReorder()

◆ infos()

◆ init_distrib()

◆ isInMagList()

◆ isMagnetic() [1/2]

◆ isMagnetic() [2/2]

◆ max_angle()

◆ maxNodes()

◆ minNodes()

◆ operator=()

◆ readMesh()

◆ readNodes()

◆ readSol()

◆ readTetraedrons()

◆ readTriangles()

◆ savesol()

◆ set()

◆ set_node_u0()

◆ set_node_zero_v()

◆ setBasis()

◆ sortNodes()

◆ surface()

◆ updateNode()

Member Data Documentation

◆ c

◆ edges

◆ fac

◆ l

◆ magFac

◆ magNode

◆ magTet

◆ node

◆ node_index

◆ paramTetra

◆ tet

◆ totalMagVol

◆ vol

◆ volumeRegions