#include <PapillaryFibreCalculator.hpp>

Collaboration diagram for PapillaryFibreCalculator:

Public Member Functions
	PapillaryFibreCalculator (TetrahedralMesh< 3, 3 > &rMesh)
std::vector< c_vector< double, 3 > >	CalculateFibreOrientations ()
Private Member Functions
c_vector< double, 3 >	GetRadiusVectorForOneElement (unsigned elementIndex)
void	GetRadiusVectors ()
void	ConstructStructureTensors ()
void	SmoothStructureTensors ()
Private Attributes
TetrahedralMesh< 3, 3 > &	mrMesh
std::vector< c_vector< double, 3 > >	mRadiusVectors
std::vector< c_matrix< double, 3, 3 > >	mStructureTensors
std::vector< c_matrix< double, 3, 3 > >	mSmoothedStructureTensors
Friends
class	TestPapillaryFibreCalculator

Detailed Description

Assigns fibre orientation vectors for papillary muscle structures. Vectors are assigned to be parallel to the direction of the axis of the muscle. To do this we use the "Structure Tensor" method.

Definition at line 46 of file PapillaryFibreCalculator.hpp.

Constructor & Destructor Documentation

PapillaryFibreCalculator::PapillaryFibreCalculator ( TetrahedralMesh< 3, 3 > & rMesh )

Constructor saves mesh and allocates memory

Parameters:

rMesh The mesh to calculate fibres on

Definition at line 113 of file PapillaryFibreCalculator.hpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetNumElements(), mRadiusVectors, mrMesh, mSmoothedStructureTensors, and mStructureTensors.

Member Function Documentation

std::vector< c_vector< double, 3 > > PapillaryFibreCalculator::CalculateFibreOrientations ( )

Main method - calculate the fibre orientations

Returns:: A fibre vector for each element

Definition at line 121 of file PapillaryFibreCalculator.hpp.

References ConstructStructureTensors(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetNumElements(), GetRadiusVectors(), mrMesh, mSmoothedStructureTensors, and SmoothStructureTensors().

void PapillaryFibreCalculator::ConstructStructureTensors ( ) [private]

This generates structure tensors from the radial vectors by taking

T = r.r'

Definition at line 183 of file PapillaryFibreCalculator.hpp.

References mRadiusVectors, and mStructureTensors.

Referenced by CalculateFibreOrientations().

c_vector< double, 3 > PapillaryFibreCalculator::GetRadiusVectorForOneElement ( unsigned elementIndex ) [private]

This method calculates the vector from the centroid of an element to all of the boundary nodes. It returns the shortest of the vectors.

Parameters:

elementIndex The index of the element we are calculating radial vectors for

Returns:: The shortest radial vector

Definition at line 140 of file PapillaryFibreCalculator.hpp.

References AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetBoundaryNodeIteratorBegin(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetBoundaryNodeIteratorEnd(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElement(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNode(), and mrMesh.

Referenced by GetRadiusVectors().

void PapillaryFibreCalculator::GetRadiusVectors ( ) [private]

This method calls GetRadiusVectorForOneElement() for each of the elements sets the radial vector for each element of the mesh.

Definition at line 171 of file PapillaryFibreCalculator.hpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorBegin(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorEnd(), GetRadiusVectorForOneElement(), mRadiusVectors, and mrMesh.

Referenced by CalculateFibreOrientations().

void PapillaryFibreCalculator::SmoothStructureTensors ( ) [private]

Smoothes the structure tensor components for each papillary element by looping over all other papillary elements, calculating distance geometric distance between the two elements; if it is within a certain limit, include this in the Gaussian kernel

Here mStructureTensors[i] is the 'rough' tensor for each element and mSmoothedStructureTensors[i] is the smoothed tensor for each element

Definition at line 191 of file PapillaryFibreCalculator.hpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorBegin(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorEnd(), mrMesh, mSmoothedStructureTensors, and mStructureTensors.

Referenced by CalculateFibreOrientations().