QuadraticMesh< DIM > Class Template Reference

#include <QuadraticMesh.hpp>

Collaboration diagram for QuadraticMesh< DIM >:

Public Member Functions
	QuadraticMesh ()
	QuadraticMesh (double spaceStep, double width, double height=0, double depth=0)
void	ConstructFromMeshReader (AbstractMeshReader< DIM, DIM > &rMeshReader)
void	ConstructFromLinearMeshReader (AbstractMeshReader< DIM, DIM > &rMeshReader)
void	WriteBoundaryElementFile (std::string directory, std::string fileName)
unsigned	GetNumVertices ()
Protected Member Functions
void	CountAndCheckVertices ()
void	AddNodesToBoundaryElements (TrianglesMeshReader< DIM, DIM > *pMeshReader)
void	AddNodeToBoundaryElement (BoundaryElement< DIM-1, DIM > pBoundaryElement, Element< DIM, DIM > pElement, unsigned internalNode)
void	AddExtraBoundaryNodes (BoundaryElement< DIM-1, DIM > pBoundaryElement, Element< DIM, DIM > pElement, unsigned nodeIndexOppositeToFace)
void	HelperMethod1 (unsigned boundaryElemNode0, unsigned boundaryElemNode1, Element< DIM, DIM > *pElement, unsigned node0, unsigned node1, unsigned node2, unsigned &rOffset, bool &rReverse)
void	HelperMethod2 (BoundaryElement< DIM-1, DIM > pBoundaryElement, Element< DIM, DIM > pElement, unsigned internalNode0, unsigned internalNode1, unsigned internalNode2, unsigned offset, bool reverse)
template<class Archive >
void	serialize (Archive &archive, const unsigned int version)
void	ConstructLinearMesh (unsigned numElemX)
void	ConstructRectangularMesh (unsigned numElemX, unsigned numElemY, bool unused=true)
void	ConstructCuboid (unsigned numElemX, unsigned numElemY, unsigned numElemZ)
Protected Attributes
unsigned	mNumVertices
Friends
class	boost::serialization::access

Detailed Description

template<unsigned DIM>
class QuadraticMesh< DIM >

A concrete quadratic mesh class that inherits from TetrahedralMesh.

Definition at line 44 of file QuadraticMesh.hpp.

Constructor & Destructor Documentation

template<unsigned DIM>

QuadraticMesh< DIM >::QuadraticMesh ( ) [inline]

Constructor

Definition at line 200 of file QuadraticMesh.hpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mMeshIsLinear.

template<unsigned DIM>

QuadraticMesh< DIM >::QuadraticMesh	(	double	spaceStep,
		double	width,
		double	height = `0`,
		double	depth = `0`
	)			`[inline]`

Create a quadratic mesh on a slab (on a line in 1D, rectangle in 2d, cuboid in 3D), with the given widths and given spacestep. In 1D height and depth need to passed in as 0 (the default value), in 2D depth must be zero

Parameters:

	spaceStep	The spatial stepsize
	width	the width of the cuboid
	height	the height of the cuboid
	depth	the depth of the cuboid

Definition at line 61 of file QuadraticMesh.cpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::ConstructRegularSlabMesh().

Member Function Documentation

template<unsigned DIM>

void QuadraticMesh< DIM >::AddExtraBoundaryNodes	(	BoundaryElement< DIM-1, DIM > *	pBoundaryElement,
		Element< DIM, DIM > *	pElement,
		unsigned	nodeIndexOppositeToFace
	)			`[inline, protected]`

Given a face in an element (defined by giving an element and the opposite node number to the face) that corresponds to a given boundary element, this method adds in the face's internal nodes to the boundary element (in the correct order).

Parameters:

	pBoundaryElement	pointer to a boundary element in the mesh
	pElement	pointer to an element in the mesh
	nodeIndexOppositeToFace	index of a node in the mesh

Definition at line 486 of file QuadraticMesh.cpp.

References QuadraticMesh< DIM >::AddNodeToBoundaryElement(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), QuadraticMesh< DIM >::HelperMethod1(), and QuadraticMesh< DIM >::HelperMethod2().

Referenced by QuadraticMesh< DIM >::AddNodesToBoundaryElements().

template<unsigned DIM>

void QuadraticMesh< DIM >::AddNodesToBoundaryElements ( TrianglesMeshReader< DIM, DIM > * pMeshReader ) [inline, protected]

Top level method for making 2D edges have 3 nodes not 2 and making 3D faces have 6 nodes not 3 (ie linear to quadratic).

Parameters:

pMeshReader

Pointer to the reader. Only used if boundaryElemFileHasContainElementInfo==true (can be null if not).

Definition at line 358 of file QuadraticMesh.cpp.

References QuadraticMesh< DIM >::AddExtraBoundaryNodes(), ElementData::ContainingElement, Node< SPACE_DIM >::ContainingElementsBegin(), Node< SPACE_DIM >::ContainingElementsEnd(), EXCEPTION, AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetBoundaryElementIteratorBegin(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetBoundaryElementIteratorEnd(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElement(), TrianglesMeshReader< ELEMENT_DIM, SPACE_DIM >::GetNextFaceData(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), TrianglesMeshReader< ELEMENT_DIM, SPACE_DIM >::GetReadContainingElementOfBoundaryElement(), and TrianglesMeshReader< ELEMENT_DIM, SPACE_DIM >::Reset().

Referenced by QuadraticMesh< DIM >::ConstructCuboid(), QuadraticMesh< DIM >::ConstructFromLinearMeshReader(), QuadraticMesh< DIM >::ConstructFromMeshReader(), and QuadraticMesh< DIM >::ConstructRectangularMesh().

template<unsigned DIM>

void QuadraticMesh< DIM >::AddNodeToBoundaryElement	(	BoundaryElement< DIM-1, DIM > *	pBoundaryElement,
		Element< DIM, DIM > *	pElement,
		unsigned	internalNode
	)			`[inline, protected]`

This method adds the given node (defined by an element and a node index) to the given boundary element, and also sets the node as a boundary element and adds it to the std::vector of boundary elements.

Parameters:

	pBoundaryElement	pointer to a boundary element in the mesh
	pElement	pointer to an element in the mesh
	internalNode	index of a node in the mesh

Definition at line 465 of file QuadraticMesh.cpp.

References AbstractElement< ELEMENT_DIM, SPACE_DIM >::AddNode(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNode(), Node< SPACE_DIM >::IsBoundaryNode(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::mBoundaryNodes, and Node< SPACE_DIM >::SetAsBoundaryNode().

Referenced by QuadraticMesh< DIM >::AddExtraBoundaryNodes(), and QuadraticMesh< DIM >::HelperMethod2().

template<unsigned DIM>

void QuadraticMesh< DIM >::ConstructCuboid	(	unsigned	numElemX,
		unsigned	numElemY,
		unsigned	numElemZ
	)			`[inline, protected, virtual]`

Create a quadratic mesh on a cuobid from (0,0,0) to (numElemX,numElemY,numElemZ) with that number of elements in each direction. This writes a temporary node file and uses tetgen to mesh this nodefile. The method overloads the equivalent method in AbstractTetrahedralMesh. This is private, users should call ConstructRegularSlabMesh();

Parameters:

	numElemX	Number of elements in x-direction (also, the width of the final mesh)
	numElemY	Number of elements in y-direction (also, the height of the final mesh)
	numElemZ	Number of elements in y-direction (also, the depth of the final mesh)

Reimplemented from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >.

Definition at line 165 of file QuadraticMesh.cpp.

References QuadraticMesh< DIM >::AddNodesToBoundaryElements(), QuadraticMesh< DIM >::CountAndCheckVertices(), TetrahedralMesh< DIM, DIM >::ImportFromMesher(), and AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mMeshIsLinear.

template<unsigned DIM>

void QuadraticMesh< DIM >::ConstructFromLinearMeshReader ( AbstractMeshReader< DIM, DIM > & rMeshReader ) [inline]

Load a quadratic mesh from a linear mesh file.

Constructs as linear mesh, then exports to triangle/tetgen

Parameters:

rMeshReader

the mesh reader

Definition at line 217 of file QuadraticMesh.cpp.

References QuadraticMesh< DIM >::AddNodesToBoundaryElements(), QuadraticMesh< DIM >::ConstructFromMeshReader(), QuadraticMesh< DIM >::CountAndCheckVertices(), TetrahedralMesh< DIM, DIM >::ExportToMesher(), TetrahedralMesh< DIM, DIM >::FreeTriangulateIo(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetNumElements(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNumNodes(), TetrahedralMesh< DIM, DIM >::ImportFromMesher(), and TetrahedralMesh< DIM, DIM >::InitialiseTriangulateIo().

Referenced by QuadraticMesh< DIM >::ConstructFromMeshReader().

template<unsigned DIM>

void QuadraticMesh< DIM >::ConstructFromMeshReader ( AbstractMeshReader< DIM, DIM > & rMeshReader ) [inline]

Load a quadratic mesh from a file.

Parameters:

rMeshReader

the mesh reader

Definition at line 267 of file QuadraticMesh.cpp.

Referenced by QuadraticMesh< DIM >::ConstructFromLinearMeshReader().

template<unsigned DIM>

void QuadraticMesh< DIM >::ConstructLinearMesh ( unsigned numElemX ) [inline, protected, virtual]

Create a quadratic mesh on the interval [0,numElemX] with numElemX elements in each direction. This is private, users should call ConstructRegularSlabMesh();

Very badly named. This creates a QUADRATIC mesh, the linear just refers to the fact the mesh is a line in 1D. The name is inherited from the parent class.

Parameters:

numElemX

Number of elements in x-direction (also, the width of the final mesh)

Reimplemented from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >.

Definition at line 67 of file QuadraticMesh.cpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mBoundaryElements, AbstractMesh< ELEMENT_DIM, SPACE_DIM >::mBoundaryNodes, AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mElements, AbstractMesh< ELEMENT_DIM, SPACE_DIM >::mNodes, QuadraticMesh< DIM >::mNumVertices, and TetrahedralMesh< DIM, DIM >::RefreshMesh().

template<unsigned DIM>

void QuadraticMesh< DIM >::ConstructRectangularMesh	(	unsigned	numElemX,
		unsigned	numElemY,
		bool	unused = `true`
	)			`[inline, protected, virtual]`

Create a quadratic mesh on a rectangle from (0,0) to (numElemX,numElemY) with that number of elements in each direction. This writes a temporary node file and uses triangle to mesh this nodefile. The method overloads the equivalent method in AbstractTetrahedralMesh. This is private, users should call ConstructRegularSlabMesh();

Parameters:

	numElemX	Number of elements in x-direction (also, the width of the final mesh)
	numElemY	Number of elements in y-direction (also, the height of the final mesh)
	unused	(defaults to true; must always be true) is for compatibility of the interface of this method with same name in AbstractTetrahedralMesh.

Reimplemented from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >.

Definition at line 113 of file QuadraticMesh.cpp.

References QuadraticMesh< DIM >::AddNodesToBoundaryElements(), QuadraticMesh< DIM >::CountAndCheckVertices(), TetrahedralMesh< DIM, DIM >::FreeTriangulateIo(), TetrahedralMesh< DIM, DIM >::ImportFromMesher(), TetrahedralMesh< DIM, DIM >::InitialiseTriangulateIo(), and AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mMeshIsLinear.

template<unsigned DIM>

void QuadraticMesh< DIM >::CountAndCheckVertices ( ) [inline, protected]

Count nodes which are vertices (not marked as internal) Check that internals appear after vertices in the list

Definition at line 43 of file QuadraticMesh.cpp.

References AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNode(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNumNodes(), and QuadraticMesh< DIM >::mNumVertices.

Referenced by QuadraticMesh< DIM >::ConstructCuboid(), QuadraticMesh< DIM >::ConstructFromLinearMeshReader(), QuadraticMesh< DIM >::ConstructFromMeshReader(), and QuadraticMesh< DIM >::ConstructRectangularMesh().

template<unsigned DIM>

unsigned QuadraticMesh< DIM >::GetNumVertices ( ) [inline]

Get the number of vertices, ie non-internal (non-quadratic), nodes.

Definition at line 210 of file QuadraticMesh.cpp.

References QuadraticMesh< DIM >::mNumVertices.

template<unsigned DIM>

void QuadraticMesh< DIM >::HelperMethod1	(	unsigned	boundaryElemNode0,
		unsigned	boundaryElemNode1,
		Element< DIM, DIM > *	pElement,
		unsigned	node0,
		unsigned	node1,
		unsigned	node2,
		unsigned &	rOffset,
		bool &	rReverse
	)			`[inline, protected]`

Nasty helper method for AddNodeToBoundaryElement() in 3D.

This method takes in the three vertices of a face which match the given boundary element, and figure out if the order of the nodes in the face is reversed in the boundary element (returned in the bool 'rReverse'). Also, the offset between the first node in the face (as given to this method) and the first node in the boundary element is computed (returned in the variable 'rOffset'). Offset should then be applied before reverse to match the face nodes to the boundary element nodes.

Todo:: document these parameters

Parameters:

	boundaryElemNode0
	boundaryElemNode1
	pElement
	node0
	node1
	node2
	rOffset
	rReverse

Definition at line 588 of file QuadraticMesh.cpp.

References AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex().

Referenced by QuadraticMesh< DIM >::AddExtraBoundaryNodes().

template<unsigned DIM>

void QuadraticMesh< DIM >::HelperMethod2	(	BoundaryElement< DIM-1, DIM > *	pBoundaryElement,
		Element< DIM, DIM > *	pElement,
		unsigned	internalNode0,
		unsigned	internalNode1,
		unsigned	internalNode2,
		unsigned	offset,
		bool	reverse
	)			`[inline, protected]`

Nasty helper method for AddNodeToBoundaryElement() in 3D.

This method takes the three internal nodes for some face in some element, applies the given offset and reverse (see HelperMethod1) to them, to get the ordered internal nodes which should given to the boundary element. It then calls AddNodeToBoundaryElement with each of the three internal nodes.