Defines an abstract cell of an unstructured volume mesh. More...

#include <MeshVizXLM/mesh/cell/MiVolumeCell.h>

Inheritance diagram for MiVolumeCell:

Public Member Functions
virtual size_t	getNumFacets () const =0
	Gets the number of facets.

virtual size_t	appendNodesIndexOfFacet (size_t facet, std::vector< size_t > &facetNodes) const =0
	Appends the array of node indices belonging to a given facet to facetNode.

virtual size_t	getSubTriangleNodesIndex (size_t facet, std::vector< size_t > &triangleNodeIds) const
	Get a triangulation of the given facet of a non linear volume cell.

virtual size_t	getSubTetrahedronNodesIndex (std::vector< size_t > &tetrahedronNodeIds) const
	Get a decomposition in sub-tetrahedra of this volume cell.

virtual size_t	appendLocalNodesIndexOfFacet (size_t facet, std::vector< size_t > &facetNodes) const
	Appends the array of local node indices belonging to a given facet to facetNodes.

virtual void	getIsosurfTopology (unsigned char caseId, std::vector< std::pair< size_t, size_t > > &edgeList) const =0
	Gets the list of topological polygons defining the part of the isosurface topology in this cell.

virtual unsigned char	getMarchingCaseId (std::vector< bool > &nodesSign, size_t beginNodeId=0) const
	Gets the isosurface "marching cube" entry according to the sign of each node in the cell.

virtual MbVec3d	getFacetCenter (size_t facet, const MiGeometryI &geometry) const
	Gets the center of the specified cell facet.

Public Member Functions inherited from MiCell
virtual	~MiCell ()
	Virtual destructor.

virtual size_t	getNumNodes () const =0
	Gets the number of nodes.

virtual size_t	getNodeIndex (size_t node) const =0
	Gets the index (in the mesh) of a node.

virtual double	getRelativeSize (const MiGeometryI *meshGeometry) const
	Gets the relative size of the cell.

virtual void	getWeight (const MbVec3d &ipcoord, std::vector< double > &weight) const
	Gets the weights of a point defined by its iso parametric coordinates.

virtual MbVec3d	getIsoParametricCoord (size_t nodeIndex) const
	Gets the iso parametric coordinates of a cell's node.

virtual MbVec3d	getCenter (const MiGeometryI &geometry) const
	Gets the center of the cell.

virtual bool	isPointInsideCell (const MiGeometryI &meshGeometry, const MbVec3d &point, std::vector< double > &weights) const
	Checks if a point is inside or outside a cell.

Protected Member Functions
virtual std::ostream &	toStream (std::ostream &s) const
	Performs an output operation on a stream which is returned afterward.

Detailed Description

Defines an abstract cell of an unstructured volume mesh.

A volume cell is a polyhedron element that connects adjacent nodes of a volume mesh.

This class defines additional methods beyond MiCell to gives a complete description of a volume cell.

Definition of the facet numbering.
See getNumFacets(), appendNodesIndexOfFacet(), appendLocalNodesIndexOfFacet().
Definition of the isosurface topology in this cell.
See getIsosurfTopology() and getMarchingCaseId().
Definition of a sub triangulation for each facet of the cell.
See getSubTriangleNodesIndex().
Definition of sub tetrahedra of the cell.
See getSubTetrahedronNodesIndex().

The methods getSubTriangleNodesIndex() and getSubTetrahedronNodesIndex() are only useful for non linear volume cells and generic polyhedron cells. For the most usual basic volume cells (tetrahedron, hexahedron, wedge, pyramid), the application only needs to implement facet numbering, the isosurface topology, and the methods defined in MiCell.

When the result of getSubTetrahedronNodesIndex() is > 0, the cell is considered to be a non linear cell, or a generic polyhedron cell. In that case, the given sub triangles and sub tetrahedra are used by various extractors (for instance MiSkinExtractUnstructured) and mesh representation classes (for instance MoMeshSkin, MoMeshCellShape). Furthermore if an MiTessellator is given, each sub triangle is tessellated by the method MiTessellator::tessellateSurfaceCell().

The utility classes MxXXXCellExtract can be used to help the application implement the various methods of MiCell and MiVolumeCell. These classes provide static implementations of the necessary methods for the most common linear and quadratic volume cells.

The implementation of all the methods of this class must be thread-safe because it might be called concurrently by multiple threads. If the implementation can't be thread-safe, MeshViz XLM must be forced to use only 1 thread. It can be done by calling MiMeshViz::init(1) or MoMeshViz::init(1) or setting the parallel argument to false for each extractor instantiation. For instance by calling MiSkinExtractUnstructured::getNewInstance(&mesh, false);

See also: MeXVolumeCell; MiVolumeTopologyExplicitI; MiTessellator; MxTetrahedronCellExtract; MxHexahedronCellExtract, MxHexahedron20CellExtract, MxHexahedron27CellExtract; MxWedgeCellExtract, MxWedge18CellExtract; MxPyramidCellExtract

Definition at line 88 of file MiVolumeCell.h.

Member Function Documentation

◆ appendLocalNodesIndexOfFacet()

virtual size_t MiVolumeCell::appendLocalNodesIndexOfFacet	(	size_t	facet,
		std::vector< size_t > &	facetNodes
	)		const

inlinevirtual

Appends the array of local node indices belonging to a given facet to facetNodes.

This method is similar to appendNodesIndexOfFacet() but facetNodes must be filled only with local node ids (i.e. node ids relative to the number of nodes in the cell): 0 <= facetNodes[i] < getNumNodes() for each i.

The method must not include the additional nodes of a non linear cell. For instance an implementation of this method for a linear or quadratic tetrahedron cell like the following one

facetNodes of the pink face of these tetrahedra.

should return 3, and facetNodes should be appended by {0,2,1} if the facet argument corresponds to the pink one. (Provided that getNodeIndex() follows this tetrahedron node numbering.)

Parameters

[in]	facet	The index of the facet (where 0 <= facet < getNumFacets()).
	facetNodes	A vector in which the local node indices of the facet will be appended. The given vector should not be cleared by this method. The implementation must use the push_back() method to add elements to the vector.

Returns: the number of nodes added.

Definition at line 174 of file MiVolumeCell.h.

◆ appendNodesIndexOfFacet()

virtual size_t MiVolumeCell::appendNodesIndexOfFacet	(	size_t	facet,
		std::vector< size_t > &	facetNodes
	)		const

pure virtual

Appends the array of node indices belonging to a given facet to facetNode.

The method must not include the additional nodes of the non linear cell. For instance an implementation of this method for a linear or quadratic tetrahedron cell like the following one

facetNodes of the pink face of these tetrahedra.

should return 3, and facetNodes should be appended by = {node0,node2,node1} if the facet argument corresponds to the pink one. (Provided that getNodeIndex() follows this tetrahedron node numbering.)

Parameters

[in]	facet	The index of the facet (where 0 <= facet < getNumFacets()).
	facetNodes	A vector in which the node indices of the facet will be appended. The given vector should not be cleared by this method. The implementation must use the push_back() method to add elements to the vector.

Returns: the number of nodes appended in facetNodes.

◆ getFacetCenter()

MbVec3d MiVolumeCell::getFacetCenter	(	size_t	facet,
		const MiGeometryI &	geometry
	)		const

inlinevirtual

Gets the center of the specified cell facet.

This method is not pure virtual because it is not used by all extraction classes. The default implementation provided computes the average of the coordinates of each facet's node.

Parameters

[in]	facet	a facet id of this cell. It must be < getNumFacets().
[in]	geometry	The geometry of the mesh. Used to retrieve the node coordinates of this cell.

Returns: The center of the specified facet.

Definition at line 271 of file MiVolumeCell.h.

◆ getIsosurfTopology()

virtual void MiVolumeCell::getIsosurfTopology	(	unsigned char	caseId,
		std::vector< std::pair< size_t, size_t > > &	edgeList
	)		const

pure virtual

Gets the list of topological polygons defining the part of the isosurface topology in this cell.

This method is used during isosurface and slice extraction on an unstructured mesh (for instance MoMeshPlaneSlice, MoMeshIsosurface, MiIsosurfExtractUnstructured). Static implementations are provided for the most common basic volume cells:

For polyhedral cells and non linear cells (for which no static implementation is available), returning an empty list can produce a hole in the isosurface. To avoid such a hole, this kind of cell must provide a decomposition into sub tetrahedra by the method getSubTetrahedronNodesIndex(). In case of polyhedral cells, one can also simply provide a polyhedral tessellator (

See also

MiTessellator::getNewTessellatorPolyhedron()) to the isosurface extractor. To summarize:

For basic cell shape (hexahedron, tetrahedron, wedge, pyramid), implement this method for example by calling MxXXXXCellExtract::getIsosurfTopology() and do not override getSubTetrahedronNodesIndex().
Otherwise, implement this method by giving an empty edgeList and either implement getSubTetrahedronNodesIndex() to provide a sub tetrahedra decomposition or for polyhedral cells only, provide a polyhedral tessellator to the isosurface extractor.

Parameters

caseId the isosurf "marching cube" entry for this cell. This is the value returned by getMarchingCaseId

edgeList the list of pairs of nodes defining the isosurface topology in this cell. The given vector should not be cleared by this method. New pairs of node ids are appended to this vector. Each pair of nodes must correspond to an edge of the cell which is cut by the isosurface. An edge is cut by the isosurface when the two nodes of the edge have different signs. The pair of nodes must be correctly ordered to define a non manifold polygon. When the isosurface part is made up of several polygons, extra polygons must be separated by specific pairs (for instance <-1,-1> for mandatory polygons and <-1,0> for optional polygons).

◆ getMarchingCaseId()

unsigned char MiVolumeCell::getMarchingCaseId	(	std::vector< bool > &	nodesSign,
		size_t	beginNodeId = `0`
	)		const

inlinevirtual

Gets the isosurface "marching cube" entry according to the sign of each node in the cell.

The return entry corresponds to the first argument of getIsosurfTopology(). This method is used during isosurface and slice extraction on an unstructured mesh (for instance MoMeshPlaneSlice, MoMeshIsosurface, MiIsosurfExtractUnstructured).

Parameters

nodesSign	A vector containing as many booleans as the number of nodes in the mesh. nodesSign[i] is relative to the sign of the (isovalue-i-th node's value) .
beginNodeId	is the first node id used by the mesh topology.

See also: MiVolumeTopologyExplicitI::getBeginNodeId.

Returns: a binary word containing as many bits as the number of nodes in the cell. the i-th bit in the returned word equals nodesSign[getNodeIndex(i)].

Note: MeshViz XLM assumes that the cell does not contain any part of an isosurf if and only if the "marching cube" entry equals 0. Thus if all the bits in this word are 1, the method returns 0.

Definition at line 254 of file MiVolumeCell.h.

◆ getNumFacets()

virtual size_t MiVolumeCell::getNumFacets ( ) const

pure virtual

Gets the number of facets.

◆ getSubTetrahedronNodesIndex()

virtual size_t MiVolumeCell::getSubTetrahedronNodesIndex ( std::vector< size_t > & tetrahedronNodeIds ) const

inlinevirtual

Get a decomposition in sub-tetrahedra of this volume cell.

This method must be overridden to correctly handle
non-linear cells. By default it returns 0 and the argument tetrahedronNodeIds is not modified. This method must not be overridden for linear cells.

Parameters

[out] tetrahedronNodeIds sub tetrahedra node id (where: 0 <= tetrahedronNodeIds[i] < getNumNodes() for each i). If the cell is not linear, tetrahedronNodeIds must contain 4 times the number of sub-tetrahedra.

Returns: the number of sub-tetrahedra in this cell. 0 by default.

Definition at line 148 of file MiVolumeCell.h.

◆ getSubTriangleNodesIndex()

virtual size_t MiVolumeCell::getSubTriangleNodesIndex	(	size_t	facet,
		std::vector< size_t > &	triangleNodeIds
	)		const

inlinevirtual

Get a triangulation of the given facet of a non linear volume cell.

This method must be overridden to correctly handle non-linear cells. For instance an implementation of this method for a quadratic tetrahedron like the following one

Triangulation of the pink face of this quadratic tetrahedron.

should return 4, and triangleNodeIds = {3,7,9, 7,0,6, 9,7,6, 9,6,2}. (provided that getNodeIndex() follows this tetrahedron node numbering).

By default the method returns 0 and the argument triangleNodeIds is not modified. It must not be overridden for linear cell.

Parameters

[in]	facet	number of the given facet. (where: 0 <= face < getNumFacets() ).
[out]	triangleNodeIds	sub triangle node id (where: 0 <= triangleNodeIds[i] < getNumNodes() for each i).

Returns: the number of sub triangles in the facet of the cell.

Definition at line 133 of file MiVolumeCell.h.

◆ toStream()

std::ostream & MiVolumeCell::toStream ( std::ostream & s ) const

inlineprotectedvirtual

Performs an output operation on a stream which is returned afterward.

This method is protected since only called by the output operator <<. This method is not pure virtual and output by default the geometry and topology of a mesh. Whenever a different output is needed, this method can be implemented in the subclass to get a different output format. No std::ostream& operator << needs to be defined in any subclass. This allows you to benefit from the virtual inheritance when using the output operator on an abstract class.

Reimplemented from MiCell.

Definition at line 283 of file MiVolumeCell.h.

The documentation for this class was generated from the following file:

MeshVizXLM/mesh/cell/MiVolumeCell.h

Open Inventor Release 2023.2.3

Public Member Functions

Protected Member Functions

Detailed Description

Member Function Documentation

◆ appendLocalNodesIndexOfFacet()

◆ appendNodesIndexOfFacet()

◆ getFacetCenter()

◆ getIsosurfTopology()

◆ getMarchingCaseId()

◆ getNumFacets()

◆ getSubTetrahedronNodesIndex()

◆ getSubTriangleNodesIndex()

◆ toStream()