Inheritance diagram for ContinuumMechanicsNeumannBcsAssembler< DIM >:

Collaboration diagram for ContinuumMechanicsNeumannBcsAssembler< DIM >:

Public Member Functions
	ContinuumMechanicsNeumannBcsAssembler (QuadraticMesh< DIM > pMesh, ContinuumMechanicsProblemDefinition< DIM > pProblemDefinition, unsigned numQuadPoints=2)
virtual	~ContinuumMechanicsNeumannBcsAssembler ()
Protected Member Functions
void	DoAssemble ()
void	AssembleOnBoundaryElement (BoundaryElement< DIM-1, DIM > &rElement, c_vector< double, STENCIL_SIZE > &rBElem, unsigned boundaryConditionIndex)
Protected Attributes
QuadraticMesh< DIM > *	mpMesh
ContinuumMechanicsProblemDefinition < DIM > *	mpProblemDefinition
GaussianQuadratureRule< DIM-1 > *	mpQuadRule
Static Private Attributes
static const unsigned	NUM_VERTICES_PER_ELEMENT = DIM
static const unsigned	NUM_NODES_PER_ELEMENT = DIM*(DIM+1)/2
static const unsigned	SPATIAL_BLOCK_SIZE_ELEMENTAL = DIM*NUM_NODES_PER_ELEMENT
static const unsigned	PRESSURE_BLOCK_SIZE_ELEMENTAL = NUM_VERTICES_PER_ELEMENT
static const unsigned	STENCIL_SIZE = DIM*NUM_NODES_PER_ELEMENT + NUM_VERTICES_PER_ELEMENT

Detailed Description

Abstract class for assembling surface-integral parts of vectors in continuum mechanics problems.

For such problems, the matrix has the form [A B1] [B2^T C ] (where often B1=B2 and C=0) and the vector has the form [b1] [b2]

This class adds surface integral components, arising from natural Neumann boundary conditions, to the spatial part of the RHS vector, ie to b1.

Constructor & Destructor Documentation

template<unsigned DIM>

ContinuumMechanicsNeumannBcsAssembler< DIM >::ContinuumMechanicsNeumannBcsAssembler	(	QuadraticMesh< DIM > *	pMesh,
		ContinuumMechanicsProblemDefinition< DIM > *	pProblemDefinition,
		unsigned	numQuadPoints = `2`
	)		`[inline]`

Constructor

Parameters:

pMesh	Pointer to the mesh
pProblemDefinition	Pointer to the problem definition object
numQuadPoints	Number of quadrature points in each direction, defaults to 2

template<unsigned DIM>

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

Destructor.

template<unsigned DIM>

void ContinuumMechanicsNeumannBcsAssembler< DIM >::AssembleOnBoundaryElement	(	BoundaryElement< DIM-1, DIM > &	rElement,
		c_vector< double, STENCIL_SIZE > &	rBElem,
		unsigned	boundaryConditionIndex
	)		`[protected]`

Calculate the contribution of a single element to the linear system.

Parameters:

rElement	The element to assemble on.
rBElem	The element's contribution to the RHS vector is returned in this vector of length n, the no. of nodes in this element. There is no need to zero this vector before calling.
boundaryConditionIndex	into the traction BCs structures in the problem definition

template<unsigned DIM>

void ContinuumMechanicsNeumannBcsAssembler< DIM >::DoAssemble ( ) [protected, virtual]

The main assembly method. Protected, should only be called through Assemble(), or AssembleVector() in the parent class.