BidomainRhsMatrixAssembler< DIM > Class Template Reference

#include <BidomainMatrixBasedAssembler.hpp>

Inherits AbstractLinearAssembler< DIM, DIM, 2, false, BidomainRhsMatrixAssembler< DIM > >.

Collaboration diagram for BidomainRhsMatrixAssembler< DIM >:

Public Member Functions
virtual 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)
virtual c_vector< double, 2 *(DIM+1)>	ComputeVectorTerm (c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 2 > &u, c_matrix< double, 2, DIM > &rGradU, Element< DIM, DIM > *pElement)
virtual c_vector< double, 2 *DIM >	ComputeVectorSurfaceTerm (const BoundaryElement< DIM-1, DIM > &rSurfaceElement, c_vector< double, DIM > &rPhi, ChastePoint< DIM > &rX)
	BidomainRhsMatrixAssembler (AbstractTetrahedralMesh< DIM, DIM > *pMesh)
	~BidomainRhsMatrixAssembler ()
Mat *	GetMatrix ()
Static Public Attributes
static const unsigned	E_DIM = DIM
static const unsigned	S_DIM = DIM
static const unsigned	P_DIM = 2u

Detailed Description

template<unsigned DIM>
class BidomainRhsMatrixAssembler< DIM >

BidomainRhsMatrixAssembler

This class only exists to construct a matrix, the matrix which is used to assemble the RHS in monodomain problems. Therefore, although it inherits from the assembler hierachy, it is not an assembler for any particular PDE problem, it is just used to assemble one matrix. Therefore only ConstructMatrixTerm is properly implemented.

The matrix that is constructed is in fact the mass matrix for a 2-unknown problem. ie ***IF*** the unknowns were ordered [V1 V2 .. V_N phi_e1 ... phi_eN ], the matrix would be [M 0] [0 M] where M_ij = integral phi_i phi_j dV, where phi_k is the k-th basis function

Definition at line 66 of file BidomainMatrixBasedAssembler.hpp.

Constructor & Destructor Documentation

template<unsigned DIM>

BidomainRhsMatrixAssembler< DIM >::BidomainRhsMatrixAssembler ( AbstractTetrahedralMesh< DIM, DIM > * pMesh ) [inline]

Constructor takes in a mesh and calls AssembleSystem to construct the matrix

Parameters:

pMesh

pointer to mesh

Definition at line 102 of file BidomainMatrixBasedAssembler.cpp.

References AbstractStaticAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, NON_HEART, BidomainRhsMatrixAssembler< DIM > >::AssembleSystem(), AbstractAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM >::mpBoundaryConditions, AbstractStaticAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, NON_HEART, BidomainRhsMatrixAssembler< DIM > >::mpLinearSystem, and AbstractStaticAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, NON_HEART, BidomainRhsMatrixAssembler< DIM > >::mpMesh.

template<unsigned DIM>

BidomainRhsMatrixAssembler< DIM >::~BidomainRhsMatrixAssembler ( ) [inline]

Destructor.

Definition at line 121 of file BidomainMatrixBasedAssembler.cpp.

References AbstractAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM >::mpBoundaryConditions.

Member Function Documentation

template<unsigned DIM>

c_matrix< double, 2 (DIM+1), 2 (DIM+1)> BidomainRhsMatrixAssembler< 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, virtual]`

Integrand in matrix definition integral (see class documentation).

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 42 of file BidomainMatrixBasedAssembler.cpp.

template<unsigned DIM>

c_vector< double, 2 *DIM > BidomainRhsMatrixAssembler< DIM >::ComputeVectorSurfaceTerm	(	const BoundaryElement< DIM-1, DIM > &	rSurfaceElement,
		c_vector< double, DIM > &	rPhi,
		ChastePoint< DIM > &	rX
	)			`[inline, virtual]`

The term arising from boundary conditions to be added to the element stiffness vector - except this class is only used fpr constructing a matrix so this is never called.

Parameters:

	rSurfaceElement	the element which is being considered.
	rPhi	The basis functions, rPhi(i) = phi_i, i=1..numBases
	rX	The point in space

Definition at line 89 of file BidomainMatrixBasedAssembler.cpp.

template<unsigned DIM>

c_vector< double, 2 *(DIM+1)> BidomainRhsMatrixAssembler< DIM >::ComputeVectorTerm	(	c_vector< double, DIM+1 > &	rPhi,
		c_matrix< double, DIM, DIM+1 > &	rGradPhi,
		ChastePoint< DIM > &	rX,
		c_vector< double, 2 > &	u,
		c_matrix< double, 2, DIM > &	rGradU,
		Element< DIM, DIM > *	pElement
	)			`[inline, virtual]`

The term to be added to the element stiffness vector - except this class is only used for constructing a matrix so this is never called.

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
	u	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