BidomainMassMatrixAssembler< DIM > Class Template Reference

#include <BidomainMassMatrixAssembler.hpp>

Inheritance diagram for BidomainMassMatrixAssembler< DIM >:

Inheritance graph
[legend]
Collaboration diagram for BidomainMassMatrixAssembler< DIM >:

Collaboration graph
[legend]

List of all members.

Public Member Functions

 BidomainMassMatrixAssembler (AbstractTetrahedralMesh< DIM, DIM > *pMesh)
 ~BidomainMassMatrixAssembler ()

Protected Member Functions

c_matrix< double, 2 *(DIM+1),
2 *(DIM+1)> 
ComputeMatrixTerm (c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 2 > &rU, c_matrix< double, 2, DIM > &rGradU, Element< DIM, DIM > *pElement)


Detailed Description

template<unsigned DIM>
class BidomainMassMatrixAssembler< DIM >

Constructs a matrix with the mass matrix in the voltage-voltage block.

Ie. IF the bidomain unknowns were ordered [V1,..,Vn,phie_1,..,phie_n], the matrix would be, in block form

[ M 0 ] [ 0 0 ]

where M is the standard nxn mass matrix.

Since the bidomain ordering is not [V1,..,Vn,phie_1,..,phie_n] but [V1,phie1,..,Vn,phie_n], the matrix has a different form.

WORKS FOR BATH PROBLEMS AS WELL AS NON-BATH PROBLEMS (sets zeros in the voltage-voltage block for bath nodes, ie [ M 0 0 0 ] [ 0 0 0 0 ] [ 0 0 0 0 ] [ 0 0 0 0 ] where the ordering is Vtissue, Vbath, phi_tissue, phi_bath

Definition at line 63 of file BidomainMassMatrixAssembler.hpp.


Constructor & Destructor Documentation

template<unsigned DIM>
BidomainMassMatrixAssembler< DIM >::BidomainMassMatrixAssembler ( AbstractTetrahedralMesh< DIM, DIM > *  pMesh  )  [inline]

Constructor

Parameters:
pMesh pointer to the mesh

Definition at line 95 of file BidomainMassMatrixAssembler.hpp.

template<unsigned DIM>
BidomainMassMatrixAssembler< DIM >::~BidomainMassMatrixAssembler (  )  [inline]

Destructor.

Definition at line 103 of file BidomainMassMatrixAssembler.hpp.


Member Function Documentation

template<unsigned DIM>
c_matrix< double, 2 *(DIM+1), 2 *(DIM+1)> BidomainMassMatrixAssembler< DIM >::ComputeMatrixTerm ( c_vector< double, DIM+1 > &  rPhi,
c_matrix< double, DIM, DIM+1 > &  rGradPhi,
ChastePoint< DIM > &  rX,
c_vector< double, 2 > &  rU,
c_matrix< double, 2, DIM > &  rGradU,
Element< DIM, DIM > *  pElement 
) [inline, protected]

ComputeMatrixTerm()

This method is called by AssembleOnElement() and tells the assembler the contribution to add to the element stiffness 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

Definition at line 35 of file BidomainMassMatrixAssembler.cpp.

References HeartRegionCode::BATH, and AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetRegion().


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

Generated on Mon Nov 1 12:36:05 2010 for Chaste by  doxygen 1.5.5