MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM > Class Template Reference

#include <MonodomainPurkinjeCableAssembler.hpp>

Inheritance diagram for MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM >:

Collaboration diagram for MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM >:

Public Member Functions
	MonodomainPurkinjeCableAssembler (MixedDimensionMesh< ELEMENT_DIM, SPACE_DIM > *pMesh)
Public Member Functions inherited from AbstractFeCableIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2, false, true, NORMAL >
	AbstractFeCableIntegralAssembler (MixedDimensionMesh< ELEMENT_DIM, SPACE_DIM > *pMesh)
virtual	~AbstractFeCableIntegralAssembler ()
Public Member Functions inherited from AbstractFeAssemblerCommon< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX, INTERPOLATION_LEVEL >
	AbstractFeAssemblerCommon ()
void	SetCurrentSolution (Vec currentSolution)
virtual	~AbstractFeAssemblerCommon ()
Public Member Functions inherited from AbstractFeAssemblerInterface< CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX >
	AbstractFeAssemblerInterface ()
void	SetMatrixToAssemble (Mat &rMatToAssemble, bool zeroMatrixBeforeAssembly=true)
void	SetVectorToAssemble (Vec &rVecToAssemble, bool zeroVectorBeforeAssembly)
void	Assemble ()
void	AssembleMatrix ()
void	AssembleVector ()
virtual	~AbstractFeAssemblerInterface ()

Private Member Functions
c_matrix< double, PROBLEM_DIM *2, PROBLEM_DIM *2 >	ComputeCableMatrixTerm (c_vector< double, 2 > &rPhi, c_matrix< double, ELEMENT_DIM, 2 > &rGradPhi, ChastePoint< SPACE_DIM > &rX, c_vector< double, PROBLEM_DIM > &rU, c_matrix< double, PROBLEM_DIM, SPACE_DIM > &rGradU, Element< 1, SPACE_DIM > *pElement)

Static Private Attributes
static const unsigned	PROBLEM_DIM =2

Additional Inherited Members
Protected Types inherited from AbstractFeCableIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2, false, true, NORMAL >
typedef LinearBasisFunction< 1 >	CableBasisFunction
Protected Member Functions inherited from AbstractFeCableIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2, false, true, NORMAL >
void	ComputeTransformedBasisFunctionDerivatives (const ChastePoint< CABLE_ELEMENT_DIM > &rPoint, const c_matrix< double, CABLE_ELEMENT_DIM, SPACE_DIM > &rInverseJacobian, c_matrix< double, SPACE_DIM, NUM_CABLE_ELEMENT_NODES > &rReturnValue)
void	DoAssemble ()
virtual c_matrix< double, PROBLEM_DIM *NUM_CABLE_ELEMENT_NODES, PROBLEM_DIM *NUM_CABLE_ELEMENT_NODES >	ComputeCableMatrixTerm (c_vector< double, NUM_CABLE_ELEMENT_NODES > &rPhi, c_matrix< double, SPACE_DIM, NUM_CABLE_ELEMENT_NODES > &rGradPhi, ChastePoint< SPACE_DIM > &rX, c_vector< double, PROBLEM_DIM > &rU, c_matrix< double, PROBLEM_DIM, SPACE_DIM > &rGradU, Element< CABLE_ELEMENT_DIM, SPACE_DIM > *pElement)
virtual c_vector< double, PROBLEM_DIM *NUM_CABLE_ELEMENT_NODES >	ComputeCableVectorTerm (c_vector< double, NUM_CABLE_ELEMENT_NODES > &rPhi, c_matrix< double, SPACE_DIM, NUM_CABLE_ELEMENT_NODES > &rGradPhi, ChastePoint< SPACE_DIM > &rX, c_vector< double, PROBLEM_DIM > &rU, c_matrix< double, PROBLEM_DIM, SPACE_DIM > &rGradU, Element< CABLE_ELEMENT_DIM, SPACE_DIM > *pElement)
virtual void	AssembleOnCableElement (Element< CABLE_ELEMENT_DIM, SPACE_DIM > &rElement, c_matrix< double, PROBLEM_DIM *NUM_CABLE_ELEMENT_NODES, PROBLEM_DIM *NUM_CABLE_ELEMENT_NODES > &rAElem, c_vector< double, PROBLEM_DIM *NUM_CABLE_ELEMENT_NODES > &rBElem)
virtual bool	ElementAssemblyCriterion (Element< CABLE_ELEMENT_DIM, SPACE_DIM > &rElement)
Protected Member Functions inherited from AbstractFeAssemblerCommon< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX, INTERPOLATION_LEVEL >
virtual double	GetCurrentSolutionOrGuessValue (unsigned nodeIndex, unsigned indexOfUnknown)
virtual void	ResetInterpolatedQuantities ()
virtual void	IncrementInterpolatedQuantities (double phiI, const Node< SPACE_DIM > *pNode)
virtual void	IncrementInterpolatedGradientQuantities (const c_matrix< double, SPACE_DIM, ELEMENT_DIM+1 > &rGradPhi, unsigned phiIndex, const Node< SPACE_DIM > *pNode)
Protected Attributes inherited from AbstractFeCableIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2, false, true, NORMAL >
MixedDimensionMesh< ELEMENT_DIM, SPACE_DIM > *	mpMesh
GaussianQuadratureRule< 1 > *	mpCableQuadRule
Protected Attributes inherited from AbstractFeAssemblerCommon< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX, INTERPOLATION_LEVEL >
ReplicatableVector	mCurrentSolutionOrGuessReplicated
Protected Attributes inherited from AbstractFeAssemblerInterface< CAN_ASSEMBLE_VECTOR, CAN_ASSEMBLE_MATRIX >
Vec	mVectorToAssemble
Mat	mMatrixToAssemble
bool	mAssembleMatrix
bool	mAssembleVector
bool	mZeroMatrixBeforeAssembly
bool	mZeroVectorBeforeAssembly
PetscInt	mOwnershipRangeLo
PetscInt	mOwnershipRangeHi
Static Protected Attributes inherited from AbstractFeCableIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2, false, true, NORMAL >
static const unsigned	CABLE_ELEMENT_DIM
static const unsigned	NUM_CABLE_ELEMENT_NODES

Detailed Description

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
class MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM >

An assembler for the purkinje part of the left hand side matrix of the monodomain-purkinje linear system. The full LHS matrix will look like (written in block form, although the code uses striping) [ chi C M/dt + K 0 ] [ 0 chi' C' M'/dt + K' ] where the top left block is the standard LHS matrix in a monodomain problem, and the bottom right block is the equivalent from integrals over cable elements.

This class assembles the matrix [ 0 0 ] [ 0 chi' C' M'/dt + K' ] The entries in this right-hand block are only non-zero on Purkinje nodes

Definition at line 59 of file MonodomainPurkinjeCableAssembler.hpp.

Constructor & Destructor Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM >::MonodomainPurkinjeCableAssembler ( MixedDimensionMesh< ELEMENT_DIM, SPACE_DIM > *  pMesh )

inline

Constructor

Parameters

pMesh

a pointer to a MixedDimensionMesh

Definition at line 120 of file MonodomainPurkinjeCableAssembler.hpp.

References EXCEPTION, MixedDimensionMesh< ELEMENT_DIM, SPACE_DIM >::GetCableElementIteratorBegin(), and MixedDimensionMesh< ELEMENT_DIM, SPACE_DIM >::GetCableElementIteratorEnd().

Member Function Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

c_matrix<double,PROBLEM_DIM*2,PROBLEM_DIM*2 > MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM >::ComputeCableMatrixTerm ( c_vector< double, 2 > &  rPhi, c_matrix< double, ELEMENT_DIM, 2 > &  rGradPhi, ChastePoint< SPACE_DIM > &  rX, c_vector< double, PROBLEM_DIM > &  rU, c_matrix< double, PROBLEM_DIM, SPACE_DIM > &  rGradU, Element< 1, SPACE_DIM > *  pElement  )

inlineprivate

Compute the cable element contribution to the matrix

Parameters

rPhi	The basis functions, rPhi(i) = phi_i, i=1..numBases.
rGradPhi	Basis gradients, rGradPhi(i,j) = d(phi_j)/d(X_i).
rX	The point in space.
rU	The unknown as a vector, u(i) = u_i.
rGradU	The gradient of the unknown as a matrix, rGradU(i,j) = d(u_i)/d(X_j).
pElement	Pointer to the element.

Returns

the stencil matrix

Definition at line 78 of file MonodomainPurkinjeCableAssembler.hpp.

References AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetAttribute(), PdeSimulationTime::GetPdeTimeStepInverse(), HeartConfig::GetPurkinjeCapacitance(), HeartConfig::GetPurkinjeConductivity(), HeartConfig::GetPurkinjeSurfaceAreaToVolumeRatio(), and HeartConfig::Instance().

Member Data Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

const unsigned MonodomainPurkinjeCableAssembler< ELEMENT_DIM, SPACE_DIM >::PROBLEM_DIM =2

staticprivate

Myocardium and purkinje potentials exist for each node, so PROBLEM_DIM=2.

Definition at line 65 of file MonodomainPurkinjeCableAssembler.hpp.

---

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

• heart/src/solver/electrics/monodomain/MonodomainPurkinjeCableAssembler.hpp