#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

◆ MonodomainPurkinjeCableAssembler()

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

◆ ComputeCableMatrixTerm()

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.