Chaste  Release::2017.1
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_NODESComputeCableMatrixTerm (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_NODESComputeCableVectorTerm (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

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
rPhiThe basis functions, rPhi(i) = phi_i, i=1..numBases.
rGradPhiBasis gradients, rGradPhi(i,j) = d(phi_j)/d(X_i).
rXThe point in space.
rUThe unknown as a vector, u(i) = u_i.
rGradUThe gradient of the unknown as a matrix, rGradU(i,j) = d(u_i)/d(X_j).
pElementPointer 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: