CellBasedPdeSolver< DIM > Class Template Reference

#include <CellBasedPdeSolver.hpp>

Inherits SimpleLinearEllipticSolver< DIM, DIM >.

Collaboration diagram for CellBasedPdeSolver< DIM >:

Public Member Functions
	CellBasedPdeSolver (TetrahedralMesh< DIM, DIM > pMesh, AbstractLinearEllipticPde< DIM, DIM > pPde, BoundaryConditionsContainer< DIM, DIM, 1 > *pBoundaryConditions)
virtual	~CellBasedPdeSolver ()
Protected Member Functions
virtual c_vector< double, 1 *(DIM+1)>	ComputeVectorTerm (c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 1 > &rU, c_matrix< double, 1, DIM > &rGradU, Element< DIM, DIM > *pElement)
virtual c_matrix< double, 1 (DIM+1), 1 (DIM+1)>	ComputeMatrixTerm (c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 1 > &rU, c_matrix< double, 1, DIM > &rGradU, Element< DIM, DIM > *pElement)
void	ResetInterpolatedQuantities ()
void	IncrementInterpolatedQuantities (double phiI, const Node< DIM > *pNode)
void	InitialiseForSolve (Vec initialSolution)
Private Attributes
double	mConstantInUSourceTerm
double	mLinearInUCoeffInSourceTerm

Detailed Description

template<unsigned DIM>
class CellBasedPdeSolver< DIM >

A purpose-made elliptic solver that interpolates the source terms from node onto Gauss points, as for a cell-based simulation with PDEs the source will only be known at the cells (nodes), not the Gauss points.

Definition at line 48 of file CellBasedPdeSolver.hpp.

Constructor & Destructor Documentation

template<unsigned DIM>

CellBasedPdeSolver< DIM >::CellBasedPdeSolver	(	TetrahedralMesh< DIM, DIM > *	pMesh,
		AbstractLinearEllipticPde< DIM, DIM > *	pPde,
		BoundaryConditionsContainer< DIM, DIM, 1 > *	pBoundaryConditions
	)			`[inline]`

Constructor.

Parameters:

	pMesh	pointer to the mesh
	pPde	pointer to the PDE
	pBoundaryConditions	pointer to the boundary conditions

Definition at line 41 of file CellBasedPdeSolver.cpp.

template<unsigned DIM>

CellBasedPdeSolver< DIM >::~CellBasedPdeSolver ( ) [inline, virtual]

Destructor.

Definition at line 49 of file CellBasedPdeSolver.cpp.

Member Function Documentation

template<unsigned DIM>

c_matrix< double, 1 (DIM+1), 1 (DIM+1)> CellBasedPdeSolver< DIM >::ComputeMatrixTerm	(	c_vector< double, DIM+1 > &	rPhi,
		c_matrix< double, DIM, DIM+1 > &	rGradPhi,
		ChastePoint< DIM > &	rX,
		c_vector< double, 1 > &	rU,
		c_matrix< double, 1, DIM > &	rGradU,
		Element< DIM, DIM > *	pElement
	)			`[inline, protected, virtual]`

The SimpleLinearEllipticSolver version of this method is overloaded using the interpolated source term.

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

Returns:: The matrix term for the stiffness matrix

Definition at line 66 of file CellBasedPdeSolver.cpp.

References CellBasedPdeSolver< DIM >::mLinearInUCoeffInSourceTerm, and SimpleLinearEllipticSolver< DIM, DIM >::mpEllipticPde.

template<unsigned DIM>

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