Chaste Release::3.1
|
#include <ExtendedBidomainMassMatrixAssembler.hpp>
Public Member Functions | |
ExtendedBidomainMassMatrixAssembler (AbstractTetrahedralMesh< DIM, DIM > *pMesh) | |
~ExtendedBidomainMassMatrixAssembler () | |
Protected Member Functions | |
virtual c_matrix< double, 3 *(DIM+1), 3 *(DIM+1)> | ComputeMatrixTerm (c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 3 > &rU, c_matrix< double, 3, DIM > &rGradU, Element< DIM, DIM > *pElement) |
Constructs a matrix with the mass matrix in the voltage-voltage block.
Ie. IF the extended bidomain unknowns were ordered [phi1_1,..,phi1_n, phi2_1, ..., phi2_n, phie_1,..,phie_n], the matrix would be, in block form
[ M 0 0] [ 0 M 0] [ 0 0 M]
where M is the standard nxn mass matrix.
Since the bidomain ordering is not [phi1_1,..,phi1_n,phi2_1,..,phi2_n, phie_1,...phie_n] but [phi1_1,phi2_1,phie_1,..,phi1_n,phi2_n,phie_n], the matrix has a different form.
Definition at line 64 of file ExtendedBidomainMassMatrixAssembler.hpp.
ExtendedBidomainMassMatrixAssembler< DIM >::ExtendedBidomainMassMatrixAssembler | ( | AbstractTetrahedralMesh< DIM, DIM > * | pMesh | ) | [inline] |
Constructor
pMesh | pointer to the mesh |
Definition at line 96 of file ExtendedBidomainMassMatrixAssembler.hpp.
ExtendedBidomainMassMatrixAssembler< DIM >::~ExtendedBidomainMassMatrixAssembler | ( | ) | [inline] |
Destructor.
Definition at line 104 of file ExtendedBidomainMassMatrixAssembler.hpp.
c_matrix< double, 3 *(DIM+1), 3 *(DIM+1)> ExtendedBidomainMassMatrixAssembler< DIM >::ComputeMatrixTerm | ( | c_vector< double, DIM+1 > & | rPhi, |
c_matrix< double, DIM, DIM+1 > & | rGradPhi, | ||
ChastePoint< DIM > & | rX, | ||
c_vector< double, 3 > & | rU, | ||
c_matrix< double, 3, DIM > & | rGradU, | ||
Element< DIM, DIM > * | pElement | ||
) | [protected, virtual] |
This method is called by AssembleOnElement() and tells the assembler the contribution to add to the element stiffness matrix.
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 |
Definition at line 41 of file ExtendedBidomainMassMatrixAssembler.cpp.
References AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetUnsignedAttribute(), and HeartRegionCode::IsRegionBath().