BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM > Class Template Reference

#include <BidomainMatrixBasedAssembler.hpp>

Inheritance diagram for BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >:

[legend]

Collaboration diagram for BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >:

[legend]

List of all members.

Public Member Functions
	BidomainMatrixBasedAssembler (AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > *pMesh, BidomainPde< SPACE_DIM > *pPde, BoundaryConditionsContainer< ELEMENT_DIM, SPACE_DIM, 2 > *pBcc, unsigned numQuadPoints=2)
	~BidomainMatrixBasedAssembler ()
virtual void	ConstructVectorForMatrixBasedRhsAssembly (Vec existingSolution)
Protected Attributes
BidomainRhsMatrixAssembler < SPACE_DIM > *	mpBidomainRhsMatrixAssembler

---

Detailed Description

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
class BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >

BidomainMatrixBasedAssembler

This class makes use of the functionality in its parent, AbstractDynamicAssemblerMixin, for specifying a constant matrix B and time-dependent vector z such that the finite element RHS vector b (in Ax=b) satisfies Bz=b, and therefore b can be constructed very quickly each timestep (compared to looping over elements and assembling it.

For the bidomain problem, B is the mass matrix for 2 unknowns (see BidomainRhsMatrixAssembler) ***IF*** the unknowns were ordered [V1 V2 .. V_N phi_e1 ... phi_eN ], then z would be

z = [z1] [z2] where z1 = CA V^{m}/dt - A I_ionic - Istim_intra and z2 = -Istim_extra,

where V is the vector of voltages are the last timestep, and I_ionic and I_stim are nodewise vectors of ionic currents and stimulus currents (and C is the capacitance and A surface-area-to-volume ratio).

Definition at line 194 of file BidomainMatrixBasedAssembler.hpp.

---

Constructor & Destructor Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::BidomainMatrixBasedAssembler	(	AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > *	pMesh,
		BidomainPde< SPACE_DIM > *	pPde,
		BoundaryConditionsContainer< ELEMENT_DIM, SPACE_DIM, 2 > *	pBcc,
		unsigned	numQuadPoints = 2
	)			[inline]

Constructor calls base constructor and creates and stores rhs-matrix.

Parameters:

	pMesh	pointer to the mesh
	pPde	pointer to the PDE
	pBcc	pointer to the boundary conditions
	numQuadPoints	number of quadrature points (defaults to 2)

Definition at line 140 of file BidomainMatrixBasedAssembler.cpp.

References BidomainRhsMatrixAssembler< DIM >::GetMatrix(), BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::mpBidomainRhsMatrixAssembler, AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >::mpMatrixForMatrixBasedRhsAssembly, AbstractStaticAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, NON_HEART, BidomainDg0Assembler< ELEMENT_DIM, SPACE_DIM > >::mpMesh, AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >::mUseMatrixBasedRhsAssembly, AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >::mVectorForMatrixBasedRhsAssembly, and AbstractCardiacPde< SPACE_DIM >::SetCacheReplication().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::~BidomainMatrixBasedAssembler

(

)

[inline]

Destructor.

Definition at line 161 of file BidomainMatrixBasedAssembler.cpp.

References BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::mpBidomainRhsMatrixAssembler, and AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >::mVectorForMatrixBasedRhsAssembly.

---

Member Function Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

void BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::ConstructVectorForMatrixBasedRhsAssembly

(

Vec

existingSolution

)

[inline, virtual]

This constructs the vector z such that b (in Ax=b) is given by Bz = b. See main class documentation.

Parameters:

existingSolution

the vector of ionic currents to use in the assembly.

Reimplemented from AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >.

Reimplemented in BidomainWithBathMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >.

Definition at line 169 of file BidomainMatrixBasedAssembler.cpp.

References DistributedVector::Begin(), DistributedVectorFactory::CreateDistributedVector(), DistributedVector::End(), HeartConfig::GetCapacitance(), HeartConfig::GetSurfaceAreaToVolumeRatio(), HeartConfig::Instance(), AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >::mDtInverse, BidomainDg0Assembler< ELEMENT_DIM, SPACE_DIM >::mpBidomainPde, AbstractStaticAssembler< ELEMENT_DIM, SPACE_DIM, PROBLEM_DIM, NON_HEART, BidomainDg0Assembler< ELEMENT_DIM, SPACE_DIM > >::mpMesh, AbstractDynamicAssemblerMixin< ELEMENT_DIM, SPACE_DIM, 2 >::mVectorForMatrixBasedRhsAssembly, and DistributedVector::Restore().

---

Member Data Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>

BidomainRhsMatrixAssembler<SPACE_DIM>* BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::mpBidomainRhsMatrixAssembler [protected]

Helper assembler for doing the RHS (matrix B)

Definition at line 199 of file BidomainMatrixBasedAssembler.hpp.

Referenced by BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::BidomainMatrixBasedAssembler(), and BidomainMatrixBasedAssembler< ELEMENT_DIM, SPACE_DIM >::~BidomainMatrixBasedAssembler().

---

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

• /tmp/release_1.1/heart/src/solver/BidomainMatrixBasedAssembler.hpp
• /tmp/release_1.1/heart/src/solver/BidomainMatrixBasedAssembler.cpp