#include <Cylindrical2dMesh.hpp>

Inheritance diagram for Cylindrical2dMesh:

Collaboration diagram for Cylindrical2dMesh:

Public Member Functions
	Cylindrical2dMesh (double width)

	Cylindrical2dMesh (double width, std::vector< Node< 2 > * > nodes)

	~Cylindrical2dMesh ()

void	ReMesh (NodeMap &rMap)

c_vector< double, 2 >	GetVectorFromAtoB (const c_vector< double, 2 > &rLocation1, const c_vector< double, 2 > &rLocation2)

void	SetNode (unsigned index, ChastePoint< 2 > point, bool concreteMove)

double	GetWidth (const unsigned &rDimension) const

unsigned	AddNode (Node< 2 > *pNewNode)

bool	GetInstanceOfMismatchedBoundaryNodes ()

Public Member Functions inherited from MutableMesh< 2, 2 >
bool	CheckIsVoronoi (double maxPenetration=0.0)

	MutableMesh ()

	MutableMesh (std::vector< Node< SPACE_DIM > * > nodes)

virtual	~MutableMesh ()

void	Clear ()

unsigned	GetNumNodes () const

unsigned	GetNumElements () const

unsigned	GetNumBoundaryElements () const

void	RescaleMeshFromBoundaryNode (ChastePoint< 1 > updatedPoint, unsigned boundaryNodeIndex)

virtual unsigned	AddNode (Node< SPACE_DIM > *pNewNode)

unsigned	AddElement (Element< ELEMENT_DIM, SPACE_DIM > *pNewElement)

virtual void	SetNode (unsigned index, ChastePoint< SPACE_DIM > point, bool concreteMove=true)

void	MoveMergeNode (unsigned index, unsigned targetIndex, bool concreteMove=true)

virtual void	DeleteNode (unsigned index)

virtual void	DeleteElement (unsigned index)

void	DeleteNodePriorToReMesh (unsigned index)

unsigned	RefineElement (Element< ELEMENT_DIM, SPACE_DIM > *pElement, ChastePoint< SPACE_DIM > point)

void	DeleteBoundaryNodeAt (unsigned index)

void	ReIndex (NodeMap &map)

void	ReMesh ()

std::vector< c_vector< unsigned, 5 > >	SplitLongEdges (double cutoffLength)

c_vector< unsigned, 3 >	SplitEdge (Node< SPACE_DIM > pNodeA, Node< SPACE_DIM > pNodeB)

Public Member Functions inherited from TetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
	TetrahedralMesh ()

void	ConstructFromMeshReader (AbstractMeshReader< ELEMENT_DIM, SPACE_DIM > &rMeshReader)

void	ReadNodesPerProcessorFile (const std::string &rNodesPerProcessorFile)

bool	CheckIsConforming ()

double	GetVolume ()

double	GetSurfaceArea ()

void	RefreshMesh ()

void	PermuteNodes ()

void	PermuteNodes (const std::vector< unsigned > &perm)

unsigned	GetContainingElementIndexWithInitialGuess (const ChastePoint< SPACE_DIM > &rTestPoint, unsigned startingElementGuess, bool strict=false)

unsigned	GetNearestElementIndex (const ChastePoint< SPACE_DIM > &rTestPoint)

std::vector< unsigned >	GetContainingElementIndices (const ChastePoint< SPACE_DIM > &rTestPoint)

double	GetAngleBetweenNodes (unsigned indexA, unsigned indexB)

virtual void	RefreshJacobianCachedData ()

virtual void	GetJacobianForElement (unsigned elementIndex, c_matrix< double, SPACE_DIM, SPACE_DIM > &rJacobian, double &rJacobianDeterminant) const

virtual void	GetInverseJacobianForElement (unsigned elementIndex, c_matrix< double, SPACE_DIM, ELEMENT_DIM > &rJacobian, double &rJacobianDeterminant, c_matrix< double, ELEMENT_DIM, SPACE_DIM > &rInverseJacobian) const

virtual void	GetWeightedDirectionForElement (unsigned elementIndex, c_vector< double, SPACE_DIM > &rWeightedDirection, double &rJacobianDeterminant) const

virtual void	GetWeightedDirectionForBoundaryElement (unsigned elementIndex, c_vector< double, SPACE_DIM > &rWeightedDirection, double &rJacobianDeterminant) const

EdgeIterator	EdgesBegin ()

EdgeIterator	EdgesEnd ()

Public Member Functions inherited from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
ElementIterator	GetElementIteratorBegin (bool skipDeletedElements=true)

ElementIterator	GetElementIteratorEnd ()

	AbstractTetrahedralMesh ()

virtual	~AbstractTetrahedralMesh ()

virtual unsigned	GetNumLocalElements () const

virtual unsigned	GetNumLocalBoundaryElements () const

unsigned	GetNumAllElements () const

unsigned	GetNumAllBoundaryElements () const

virtual unsigned	GetNumCableElements () const

virtual unsigned	GetNumVertices () const

virtual unsigned	GetMaximumNodeIndex ()

Element< ELEMENT_DIM, SPACE_DIM > *	GetElement (unsigned index) const

BoundaryElement< ELEMENT_DIM-1, SPACE_DIM > *	GetBoundaryElement (unsigned index) const

void	ConstructFromMesh (AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > &rOtherMesh)

BoundaryElementIterator	GetBoundaryElementIteratorBegin () const

BoundaryElementIterator	GetBoundaryElementIteratorEnd () const

void	CheckOutwardNormals ()

virtual void	ConstructLinearMesh (unsigned width)

virtual void	ConstructRectangularMesh (unsigned width, unsigned height, bool stagger=true)

virtual void	ConstructCuboid (unsigned width, unsigned height, unsigned depth)

void	ConstructRegularSlabMesh (double spaceStep, double width, double height=0, double depth=0)

void	ConstructRegularSlabMeshWithDimensionSplit (unsigned dimension, double spaceStep, double width, double height=0, double depth=0)

virtual bool	CalculateDesignatedOwnershipOfBoundaryElement (unsigned faceIndex)

virtual bool	CalculateDesignatedOwnershipOfElement (unsigned elementIndex)

unsigned	CalculateMaximumNodeConnectivityPerProcess () const

virtual void	GetHaloNodeIndices (std::vector< unsigned > &rHaloIndices) const

void	CalculateNodeExchange (std::vector< std::vector< unsigned > > &rNodesToSendPerProcess, std::vector< std::vector< unsigned > > &rNodesToReceivePerProcess)

virtual c_vector< double, 2 >	CalculateMinMaxEdgeLengths ()

unsigned	GetContainingElementIndex (const ChastePoint< SPACE_DIM > &rTestPoint, bool strict=false, std::set< unsigned > testElements=std::set< unsigned >(), bool onlyTryWithTestElements=false)

unsigned	GetNearestElementIndexFromTestElements (const ChastePoint< SPACE_DIM > &rTestPoint, std::set< unsigned > testElements)

Public Member Functions inherited from AbstractMesh< ELEMENT_DIM, SPACE_DIM >
NodeIterator	GetNodeIteratorBegin (bool skipDeletedNodes=true)

NodeIterator	GetNodeIteratorEnd ()

	AbstractMesh ()

virtual	~AbstractMesh ()

unsigned	GetNumBoundaryNodes () const

virtual unsigned	GetNumAllNodes () const

unsigned	GetNumNodeAttributes () const

Node< SPACE_DIM > *	GetNode (unsigned index) const

virtual Node< SPACE_DIM > *	GetNodeOrHaloNode (unsigned index) const

Node< SPACE_DIM > *	GetNodeFromPrePermutationIndex (unsigned index) const

virtual DistributedVectorFactory *	GetDistributedVectorFactory ()

virtual void	SetDistributedVectorFactory (DistributedVectorFactory *pFactory)

BoundaryNodeIterator	GetBoundaryNodeIteratorBegin () const

BoundaryNodeIterator	GetBoundaryNodeIteratorEnd () const

std::string	GetMeshFileBaseName () const

bool	IsMeshOnDisk () const

const std::vector< unsigned > &	rGetNodePermutation () const

virtual c_vector< double, SPACE_DIM >	GetVectorFromAtoB (const c_vector< double, SPACE_DIM > &rLocationA, const c_vector< double, SPACE_DIM > &rLocationB)

double	GetDistanceBetweenNodes (unsigned indexA, unsigned indexB)

virtual ChasteCuboid< SPACE_DIM >	CalculateBoundingBox () const

virtual unsigned	GetNearestNodeIndex (const ChastePoint< SPACE_DIM > &rTestPoint)

virtual void	Scale (const double xFactor=1.0, const double yFactor=1.0, const double zFactor=1.0)

virtual void	Translate (const c_vector< double, SPACE_DIM > &rDisplacement)

void	Translate (const double xMovement=0.0, const double yMovement=0.0, const double zMovement=0.0)

virtual void	Rotate (c_matrix< double, SPACE_DIM, SPACE_DIM > rotationMatrix)

void	Rotate (c_vector< double, 3 > axis, double angle)

void	RotateX (const double theta)

void	RotateY (const double theta)

void	RotateZ (const double theta)

void	Rotate (double theta)

bool	IsMeshChanging () const

unsigned	CalculateMaximumContainingElementsPerProcess () const

void	SetMeshHasChangedSinceLoading ()

Private Member Functions
void	UpdateTopAndBottom ()

void	CreateHaloNodes ()

void	CreateMirrorNodes ()

void	ReconstructCylindricalMesh ()

void	DeleteHaloNodes ()

void	CorrectNonPeriodicMesh ()

void	GenerateVectorsOfElementsStraddlingPeriodicBoundaries ()

unsigned	GetCorrespondingNodeIndex (unsigned nodeIndex)

void	UseTheseElementsToDecideMeshing (std::set< unsigned > &rMainSideElements)

template<class Archive >
void	serialize (Archive &archive, const unsigned int version)

Private Attributes
double	mWidth

double	mTop

double	mBottom

std::vector< unsigned >	mLeftOriginals

std::vector< unsigned >	mLeftImages

std::map< unsigned, unsigned >	mImageToLeftOriginalNodeMap

std::vector< unsigned >	mRightOriginals

std::vector< unsigned >	mRightImages

std::map< unsigned, unsigned >	mImageToRightOriginalNodeMap

std::set< unsigned >	mLeftPeriodicBoundaryElementIndices

std::set< unsigned >	mRightPeriodicBoundaryElementIndices

std::vector< unsigned >	mTopHaloNodes

std::vector< unsigned >	mBottomHaloNodes

bool	mMismatchedBoundaryElements

Friends
class	TestCylindrical2dMesh

class	boost::serialization::access

Additional Inherited Members
Public Types inherited from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
typedef std::vector< BoundaryElement< ELEMENT_DIM-1, SPACE_DIM > * >::const_iterator	BoundaryElementIterator

Public Types inherited from AbstractMesh< ELEMENT_DIM, SPACE_DIM >
typedef std::vector< Node< SPACE_DIM > * >::const_iterator	BoundaryNodeIterator

Protected Member Functions inherited from TetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
unsigned	SolveNodeMapping (unsigned index) const

unsigned	SolveElementMapping (unsigned index) const

unsigned	SolveBoundaryElementMapping (unsigned index) const

template<class MESHER_IO >
void	ExportToMesher (NodeMap &map, MESHER_IO &mesherInput, int *elementList=nullptr)

template<class MESHER_IO >
void	ImportFromMesher (MESHER_IO &mesherOutput, unsigned numberOfElements, int elementList, unsigned numberOfFaces, int faceList, int *edgeMarkerList)

void	InitialiseTriangulateIo (triangulateio &mesherIo)

void	FreeTriangulateIo (triangulateio &mesherIo)

Protected Member Functions inherited from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
void	SetElementOwnerships ()

Protected Member Functions inherited from AbstractMesh< ELEMENT_DIM, SPACE_DIM >
ChasteCuboid< SPACE_DIM >	CalculateBoundingBox (const std::vector< Node< SPACE_DIM > * > &rNodes) const

Protected Attributes inherited from MutableMesh< 2, 2 >
std::vector< unsigned >	mDeletedElementIndices

std::vector< unsigned >	mDeletedBoundaryElementIndices

std::vector< unsigned >	mDeletedNodeIndices

bool	mAddedNodes

Protected Attributes inherited from TetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
std::vector< c_vector< double, SPACE_DIM > >	mElementWeightedDirections

std::vector< c_matrix< double, SPACE_DIM, ELEMENT_DIM > >	mElementJacobians

std::vector< c_matrix< double, ELEMENT_DIM, SPACE_DIM > >	mElementInverseJacobians

std::vector< double >	mElementJacobianDeterminants

std::vector< c_vector< double, SPACE_DIM > >	mBoundaryElementWeightedDirections

std::vector< double >	mBoundaryElementJacobianDeterminants

Protected Attributes inherited from AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >
bool	mMeshIsLinear

std::vector< Element< ELEMENT_DIM, SPACE_DIM > * >	mElements

std::vector< BoundaryElement< ELEMENT_DIM-1, SPACE_DIM > * >	mBoundaryElements

Protected Attributes inherited from AbstractMesh< ELEMENT_DIM, SPACE_DIM >
std::vector< Node< SPACE_DIM > * >	mNodes

std::vector< Node< SPACE_DIM > * >	mBoundaryNodes

DistributedVectorFactory *	mpDistributedVectorFactory

std::vector< unsigned >	mNodePermutation

std::string	mMeshFileBaseName

bool	mMeshChangesDuringSimulation

Detailed Description

A subclass of MutableMesh<2,2> for a rectangular mesh with periodic left and right boundaries, representing a cylindrical geometry.

The class works by overriding calls such as ReMesh() and GetVectorFromAtoB() so that simulation classes can treat this class in exactly the same way as a MutableMesh<2,2>.

Definition at line 55 of file Cylindrical2dMesh.hpp.

Constructor & Destructor Documentation

Cylindrical2dMesh::Cylindrical2dMesh ( double width )

Constructor.

Parameters

width the width of the mesh (circumference)

Definition at line 46 of file Cylindrical2dMesh.cpp.

Referenced by serialize().

Cylindrical2dMesh::Cylindrical2dMesh	(	double	width,
		std::vector< Node< 2 > * >	nodes
	)

A constructor which reads in a width and collection of nodes, then calls a ReMesh() command to create the elements of the mesh.

Parameters

width	the periodic length scale
nodes	a collection of nodes to construct the mesh with

Definition at line 57 of file Cylindrical2dMesh.cpp.

References AbstractMesh< ELEMENT_DIM, SPACE_DIM >::mNodes, MutableMesh< 2, 2 >::ReMesh(), Node< SPACE_DIM >::rGetLocation(), and UNUSED_OPT.

Cylindrical2dMesh::~Cylindrical2dMesh ( )

Destructor.

Definition at line 53 of file Cylindrical2dMesh.cpp.

Referenced by serialize().

Member Function Documentation

unsigned Cylindrical2dMesh::AddNode ( Node< 2 > * pNewNode )

Overridden AddNode() method.

Parameters

pNewNode the node to be added to the mesh

Returns: the global index of the new node

Definition at line 528 of file Cylindrical2dMesh.cpp.

References MutableMesh< ELEMENT_DIM, SPACE_DIM >::AddNode(), Node< SPACE_DIM >::GetPoint(), and SetNode().

Referenced by serialize().

void Cylindrical2dMesh::CorrectNonPeriodicMesh ( )

private

This method should only ever be called by the public ReMesh() method.

Uses mLeftPeriodicBoundaryElementIndices and mRightPeriodicBoundaryElementIndices and compares the nodes in each to ensure that both boundaries have been meshed identically. If they have not it calls UseTheseElementsToDecideMeshing() to sort out the troublesome elements which have been meshed differently on each side and uses the meshing of the elements on the right hand boundary to decide on how to mesh the left hand side.

Definition at line 539 of file Cylindrical2dMesh.cpp.

References GenerateVectorsOfElementsStraddlingPeriodicBoundaries(), GetCorrespondingNodeIndex(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElement(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), mLeftPeriodicBoundaryElementIndices, mRightPeriodicBoundaryElementIndices, NEVER_REACHED, and UseTheseElementsToDecideMeshing().

Referenced by ReMesh().

void Cylindrical2dMesh::CreateHaloNodes ( )

private

This method creates a compressed row of nodes, just above and below the main mesh at a constant height (a 'halo'). These will mesh to a known configuration (each one on the boundary), this avoids boundary elements of lengths over half the mesh width, which prevent the process of cylindrical meshing.

The nodes which are created are later removed by DeleteHaloNodes().

Definition at line 155 of file Cylindrical2dMesh.cpp.

References MutableMesh< ELEMENT_DIM, SPACE_DIM >::AddNode(), mBottom, mBottomHaloNodes, mTop, mTopHaloNodes, mWidth, and UpdateTopAndBottom().

Referenced by ReMesh().

void Cylindrical2dMesh::CreateMirrorNodes ( )

private

Creates a set of mirrored nodes for a cylindrical re-mesh. Updates mRightImages and mLeftImages. All mesh points should be 0 < x < mWidth.

This method should only ever be called by the public ReMesh() method.

Definition at line 88 of file Cylindrical2dMesh.cpp.

References MutableMesh< ELEMENT_DIM, SPACE_DIM >::AddNode(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNodeIteratorBegin(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNodeIteratorEnd(), mImageToLeftOriginalNodeMap, mImageToRightOriginalNodeMap, mLeftImages, mLeftOriginals, mLeftPeriodicBoundaryElementIndices, AbstractMesh< ELEMENT_DIM, SPACE_DIM >::mNodes, mRightImages, mRightOriginals, mRightPeriodicBoundaryElementIndices, and mWidth.

Referenced by ReMesh().

void Cylindrical2dMesh::DeleteHaloNodes ( )

private

This method should only ever be called by the public ReMesh method.

This method removes the nodes which were added by CreateHaloNodes() before the remeshing algorithm was called.

Definition at line 463 of file Cylindrical2dMesh.cpp.

References MutableMesh< 2, 2 >::DeleteBoundaryNodeAt(), mBottomHaloNodes, and mTopHaloNodes.

Referenced by ReMesh().

void Cylindrical2dMesh::GenerateVectorsOfElementsStraddlingPeriodicBoundaries ( )

private

This method should only ever be called by the public ReMesh method.

The elements which straddle the periodic boundaries need to be identified in order to compare the list on the right with the list on the left and reconstruct a cylindrical mesh.

Empties and repopulates the member variables mLeftPeriodicBoundaryElementIndices and mRightPeriodicBoundaryElementIndices

Definition at line 731 of file Cylindrical2dMesh.cpp.

References AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorBegin(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorEnd(), mImageToLeftOriginalNodeMap, mImageToRightOriginalNodeMap, mLeftPeriodicBoundaryElementIndices, and mRightPeriodicBoundaryElementIndices.

Referenced by CorrectNonPeriodicMesh().

unsigned Cylindrical2dMesh::GetCorrespondingNodeIndex ( unsigned nodeIndex )

private

This method should only ever be called by the public ReMesh() method.

Parameters

nodeIndex The index of an original/mirrored node

Returns: the index of the corresponding mirror image of that node (can be either an original or mirror node)

Definition at line 804 of file Cylindrical2dMesh.cpp.

References CHASTE_CLASS_EXPORT, mLeftImages, mLeftOriginals, mRightImages, and mRightOriginals.

Referenced by CorrectNonPeriodicMesh(), and UseTheseElementsToDecideMeshing().

bool Cylindrical2dMesh::GetInstanceOfMismatchedBoundaryNodes ( )

Returns: whether you have mismatched numbers of left and right boundary nodes

Referenced by serialize().

c_vector< double, 2 > Cylindrical2dMesh::GetVectorFromAtoB	(	const c_vector< double, 2 > &	rLocation1,
		const c_vector< double, 2 > &	rLocation2
	)

Overridden GetVectorFromAtoB() method.

Evaluates the (surface) distance between two points in a 2D cylindrical geometry.

Parameters

rLocation1	the x and y co-ordinates of point 1
rLocation2	the x and y co-ordinates of point 2

Returns: the vector from location1 to location2

Definition at line 473 of file Cylindrical2dMesh.cpp.

References mWidth.

Referenced by serialize().

double Cylindrical2dMesh::GetWidth ( const unsigned & rDimension ) const

virtual

Overridden GetWidth() method.

Calculate the 'width' of any dimension of the mesh, taking periodicity into account.

Parameters

rDimension a dimension (0 or 1)

Returns: The maximum distance between any nodes in this dimension.

Reimplemented from AbstractMesh< ELEMENT_DIM, SPACE_DIM >.

Definition at line 513 of file Cylindrical2dMesh.cpp.

References AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetWidth(), and mWidth.

Referenced by serialize().

void Cylindrical2dMesh::ReconstructCylindricalMesh ( )

private

After any corrections have been made to the boundary elements (see UseTheseElementsToDecideMeshing()) this method deletes the mirror image nodes, elements and boundary elements created for a cylindrical remesh by cycling through the elements and changing elements with partly real and partly imaginary elements to be real with periodic real nodes instead of mirror image nodes. We end up with very strangely shaped elements which cross the whole mesh but specify the correct connections between nodes.

This method should only ever be called by the public ReMesh() method.

Definition at line 329 of file Cylindrical2dMesh.cpp.

Referenced by ReMesh().

void Cylindrical2dMesh::ReMesh ( NodeMap & rMap )

virtual

Overridden ReMesh() method.

Conduct a cylindrical remesh by calling CreateMirrorNodes() to create mirror image nodes, then calling ReMesh() on the parent class, then mapping the new node indices and calling ReconstructCylindricalMesh() to remove surplus nodes, leaving a fully periodic mesh.

Parameters

rMap	a reference to a nodemap which should be created with the required number of nodes.

Reimplemented from MutableMesh< 2, 2 >.

Definition at line 184 of file Cylindrical2dMesh.cpp.

template<class Archive >

void Cylindrical2dMesh::serialize	(	Archive &	archive,
		const unsigned int	version
	)

inlineprivate

Archives the member variables of the Cylindrical2dMesh class which have to be preserved during the lifetime of the mesh.

The remaining member variables are re-initialised before being used by each ReMesh() call so they do not need to be archived.

Parameters

archive	the archive
version	the current version of this class the current version of this class

Definition at line 209 of file Cylindrical2dMesh.hpp.

References AddNode(), CHASTE_CLASS_EXPORT, Cylindrical2dMesh(), GetInstanceOfMismatchedBoundaryNodes(), GetVectorFromAtoB(), GetWidth(), mBottom, mTop, mWidth, MutableMesh< 2, 2 >::ReMesh(), SetNode(), and ~Cylindrical2dMesh().

void Cylindrical2dMesh::SetNode	(	unsigned	index,
		ChastePoint< 2 >	point,
		bool	concreteMove
	)

Overridden SetNode() method.

If the location should be set outside a cylindrical boundary, it is moved back onto the cylinder.

Parameters

index	is the index of the node to be moved
point	is the new target location of the node
concreteMove	is set to false if we want to skip the signed area tests

Definition at line 495 of file Cylindrical2dMesh.cpp.

References mWidth, ChastePoint< DIM >::rGetLocation(), ChastePoint< DIM >::SetCoordinate(), and MutableMesh< ELEMENT_DIM, SPACE_DIM >::SetNode().

Referenced by AddNode(), and serialize().

void Cylindrical2dMesh::UpdateTopAndBottom ( )

private

Calls AbstractMesh<2,2>::CalculateBoundingBox() to calculate mTop and mBottom for the cylindrical mesh.

This method should only ever be called by the public ReMesh() method.

Definition at line 81 of file Cylindrical2dMesh.cpp.

References AbstractMesh< ELEMENT_DIM, SPACE_DIM >::CalculateBoundingBox(), mBottom, mTop, ChasteCuboid< SPACE_DIM >::rGetLowerCorner(), and ChasteCuboid< SPACE_DIM >::rGetUpperCorner().

Referenced by CreateHaloNodes().

void Cylindrical2dMesh::UseTheseElementsToDecideMeshing ( std::set< unsigned > & rMainSideElements )

private

This method takes in two elements which are not meshed in the same way on the opposite boundary. It deletes the corresponding two elements (connecting the same four nodes) and makes two new elements which are connected in the same way. We should then be able to reconstruct the cylindrical mesh properly.

Parameters

rMainSideElements two elements (usually in a square) which have been meshed differently on the opposite boundary

Definition at line 649 of file Cylindrical2dMesh.cpp.

References GetCorrespondingNodeIndex(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElement(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorBegin(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetElementIteratorEnd(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNode(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetNumAllElements(), AbstractMesh< ELEMENT_DIM, SPACE_DIM >::GetNumAllNodes(), MutableMesh< 2, 2 >::mDeletedElementIndices, TetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mElementInverseJacobians, TetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mElementJacobianDeterminants, TetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mElementJacobians, AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::mElements, and MutableMesh< 2, 2 >::ReIndex().