Inheritance diagram for ContinuumMechanicsNeumannBcsAssembler< DIM >:

Collaboration diagram for ContinuumMechanicsNeumannBcsAssembler< DIM >:

Public Member Functions
	ContinuumMechanicsNeumannBcsAssembler (AbstractTetrahedralMesh< DIM, DIM > pMesh, ContinuumMechanicsProblemDefinition< DIM > pProblemDefinition)

virtual	~ContinuumMechanicsNeumannBcsAssembler ()

Public Member Functions inherited from AbstractFeAssemblerInterface< true, false >
	AbstractFeAssemblerInterface ()

void	SetMatrixToAssemble (Mat &rMatToAssemble, bool zeroMatrixBeforeAssembly=true)

void	SetVectorToAssemble (Vec &rVecToAssemble, bool zeroVectorBeforeAssembly)

void	Assemble ()

void	AssembleMatrix ()

void	AssembleVector ()

virtual	~AbstractFeAssemblerInterface ()

Protected Member Functions
void	DoAssemble ()

void	AssembleOnBoundaryElement (BoundaryElement< DIM-1, DIM > &rElement, c_vector< double, STENCIL_SIZE > &rBElem, unsigned boundaryConditionIndex)

Protected Attributes
AbstractTetrahedralMesh< DIM, DIM > *	mpMesh

ContinuumMechanicsProblemDefinition< DIM > *	mpProblemDefinition

GaussianQuadratureRule< DIM-1 > *	mpQuadRule

Protected Attributes inherited from AbstractFeAssemblerInterface< true, false >
Vec	mVectorToAssemble

Mat	mMatrixToAssemble

bool	mAssembleMatrix

bool	mAssembleVector

bool	mZeroMatrixBeforeAssembly

bool	mZeroVectorBeforeAssembly

PetscInt	mOwnershipRangeLo

PetscInt	mOwnershipRangeHi

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

template<unsigned DIM>
class ContinuumMechanicsNeumannBcsAssembler< DIM >

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	(	AbstractTetrahedralMesh< DIM, DIM > *	pMesh,
		ContinuumMechanicsProblemDefinition< DIM > *	pProblemDefinition
	)

inline

Constructor

Parameters

pMesh	Pointer to the mesh
pProblemDefinition	Pointer to the problem definition object

template<unsigned DIM>

virtual ContinuumMechanicsNeumannBcsAssembler< DIM >::~ContinuumMechanicsNeumannBcsAssembler ( )

inlinevirtual

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