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, unsigned numQuadPoints=2)
	~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 42 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,
		unsigned	numQuadPoints = `2`
	)			`[inline]`

Constructor.

Parameters:

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

Definition at line 35 of file CellBasedPdeSolver.cpp.

template<unsigned DIM>

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

Destructor.

Definition at line 44 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 61 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]`