Detailed Description

This class represents a 3d prism element. The following sketch shows the node and edge numbering.

       5
      / \
     / : \
   8/  :  \7
   /   :5  \
  /    :  6 \
 3-----------4
 |     :     |
 |     2     |
 |    . .    |
3|   .   .   |4
 | 2.     .1 |
 | .       . |
 |.         .|
 0-----------1
       0

Definition at line 44 of file PrismRule6.h.

#include <PrismRule6.h>

Inheritance diagram for MeshLib::PrismRule6:

[legend]

Collaboration diagram for MeshLib::PrismRule6:

[legend]

Public Types
using	EdgeReturn = MeshLib::LinearEdgeReturn
	Returns the i-th edge of the element. More...

Static Public Member Functions
static const Element *	getFace (const Element *e, unsigned i)
	Returns the i-th face of the element. More...

static bool	isPntInElement (Node const const nodes, MathLib::Point3d const &pnt, double eps)

static ElementErrorCode	validate (const Element *e)

static unsigned	identifyFace (Node const const , Node const *nodes[3])
	Returns the ID of a face given an array of nodes. More...

static double	computeVolume (Node const const _nodes)
	Calculates the volume of a convex hexahedron by partitioning it into six tetrahedra. More...

Static Public Member Functions inherited from MeshLib::CellRule
static bool	testElementNodeOrder (Element const &e)

Static Public Attributes
static const unsigned	n_base_nodes = 6u
	Constant: The number of base nodes for this element. More...

static const unsigned	n_all_nodes = 6u
	Constant: The number of all nodes for this element. More...

static const MeshElemType	mesh_elem_type = MeshElemType::PRISM
	Constant: The geometric type of the element. More...

static const CellType	cell_type = CellType::PRISM6
	Constant: The FEM type of the element. More...

static const unsigned	n_faces = 5
	Constant: The number of faces. More...

static const unsigned	n_edges = 9
	Constant: The number of edges. More...

static const unsigned	n_neighbors = 5
	Constant: The number of neighbors. More...

static const unsigned	face_nodes [5][4]
	Constant: Local node index table for faces. More...

static const unsigned	edge_nodes [9][2]
	Constant: Local node index table for edge. More...

static const unsigned	n_face_nodes [5] = {3, 4, 4, 4, 3}
	Constant: Table for the number of nodes for each face. More...

Static Public Attributes inherited from MeshLib::CellRule
static const unsigned	dimension = 3u
	Constant: Dimension of this mesh element. More...

Member Typedef Documentation

◆ EdgeReturn

using MeshLib::PrismRule6::EdgeReturn = MeshLib::LinearEdgeReturn

Returns the i-th edge of the element.

Definition at line 78 of file PrismRule6.h.

Member Function Documentation

◆ computeVolume()

double MeshLib::PrismRule6::computeVolume ( Node const *const * _nodes )

static

Calculates the volume of a convex hexahedron by partitioning it into six tetrahedra.

Definition at line 65 of file PrismRule6.cpp.

 {
     return MathLib::calcTetrahedronVolume(
                *_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3]) +
            MathLib::calcTetrahedronVolume(
                *_nodes[1], *_nodes[4], *_nodes[2], *_nodes[3]) +
            MathLib::calcTetrahedronVolume(
                *_nodes[2], *_nodes[4], *_nodes[5], *_nodes[3]);
 }

References MathLib::calcTetrahedronVolume().

◆ getFace()

const Element * MeshLib::PrismRule6::getFace	(	const Element *	e,
		unsigned	i
	)

static

Returns the i-th face of the element.

Definition at line 43 of file PrismRule6.cpp.

 {
     if (i < n_faces)
     {
         unsigned nFaceNodes(PrismRule6::n_face_nodes[i]);
         auto** nodes = new Node*[nFaceNodes];
         for (unsigned j = 0; j < nFaceNodes; j++)
         {
             nodes[j] = const_cast<Node*>(e->getNode(face_nodes[i][j]));
         }
  
         if (i == 0 || i == 4)
         {
             return new Tri(nodes, e->getID());
         }
  
         return new Quad(nodes);
     }
     ERR("Error in MeshLib::Element::getFace() - Index {:d} does not exist.", i);
     return nullptr;
 }

References ERR(), face_nodes, MeshLib::Element::getID(), MeshLib::Element::getNode(), n_face_nodes, and n_faces.

◆ identifyFace()

unsigned MeshLib::PrismRule6::identifyFace	(	Node const const	_nodes,
		Node const *	nodes[3]
	)

static

Returns the ID of a face given an array of nodes.

Definition at line 87 of file PrismRule6.cpp.

 {
     for (unsigned i = 0; i < 5; i++)
     {
         unsigned flag(0);
         for (unsigned j = 0; j < 4; j++)
         {
             for (unsigned k = 0; k < 3; k++)
             {
                 if (face_nodes[i][j] != 99 &&
                     _nodes[face_nodes[i][j]] == nodes[k])
                 {
                     flag++;
                 }
             }
         }
         if (flag == 3)
         {
             return i;
         }
     }
     return std::numeric_limits<unsigned>::max();
 }

References face_nodes.

◆ isPntInElement()

bool MeshLib::PrismRule6::isPntInElement	(	Node const const	nodes,
		MathLib::Point3d const &	pnt,
		double	eps
	)

static

Checks if a point is inside the element.

Parameters

pnt	a 3D MathLib::Point3d object
eps	tolerance for numerical algorithm used or computing the property

Returns: true if the point is not outside the element, false otherwise

Parameters

nodes the nodes of the element.

Definition at line 75 of file PrismRule6.cpp.

 {
     return (MathLib::isPointInTetrahedron(
                 pnt, *nodes[0], *nodes[1], *nodes[2], *nodes[3], eps) ||
             MathLib::isPointInTetrahedron(
                 pnt, *nodes[1], *nodes[4], *nodes[2], *nodes[3], eps) ||
             MathLib::isPointInTetrahedron(
                 pnt, *nodes[2], *nodes[4], *nodes[5], *nodes[3], eps));
 }

References MathLib::isPointInTetrahedron().

◆ validate()

ElementErrorCode MeshLib::PrismRule6::validate ( const Element * e )

static

Tests if the element is geometrically valid.

Definition at line 112 of file PrismRule6.cpp.

 {
     ElementErrorCode error_code;
     error_code[ElementErrorFlag::ZeroVolume] = hasZeroVolume(*e);
  
     for (unsigned i = 1; i < 4; ++i)
     {
         const auto* quad(dynamic_cast<const MeshLib::Quad*>(e->getFace(i)));
         if (quad)
         {
             error_code |= quad->validate();
         }
         else
         {
             error_code.set(ElementErrorFlag::NodeOrder);
         }
         delete quad;
     }
     error_code[ElementErrorFlag::NodeOrder] = !e->testElementNodeOrder();
     return error_code;
 }

References MeshLib::Element::getFace(), MeshLib::hasZeroVolume(), NodeOrder, ElementErrorCode::set(), MeshLib::Element::testElementNodeOrder(), and ZeroVolume.

Member Data Documentation

◆ cell_type

const CellType MeshLib::PrismRule6::cell_type = CellType::PRISM6

static

Constant: The FEM type of the element.

Definition at line 57 of file PrismRule6.h.

◆ edge_nodes

const unsigned MeshLib::PrismRule6::edge_nodes

static

Initial value:

= {
    {0, 1},  
    {1, 2},  
    {0, 2},  
    {0, 3},  
    {1, 4},  
    {2, 5},  
    {3, 4},  
    {4, 5},  
    {3, 5}   
}

Constant: Local node index table for edge.

Definition at line 72 of file PrismRule6.h.

◆ face_nodes

const unsigned MeshLib::PrismRule6::face_nodes

static

Initial value:

= {
    {0, 2, 1, 99},  
    {0, 1, 4, 3},   
    {1, 2, 5, 4},   
    {2, 0, 3, 5},   
    {3, 4, 5, 99}   
}

Constant: Local node index table for faces.

Definition at line 69 of file PrismRule6.h.

Referenced by getFace(), and identifyFace().

◆ mesh_elem_type

const MeshElemType MeshLib::PrismRule6::mesh_elem_type = MeshElemType::PRISM

static

Constant: The geometric type of the element.

Definition at line 54 of file PrismRule6.h.

◆ n_all_nodes

const unsigned MeshLib::PrismRule6::n_all_nodes = 6u

static

Constant: The number of all nodes for this element.

Definition at line 51 of file PrismRule6.h.

◆ n_base_nodes

const unsigned MeshLib::PrismRule6::n_base_nodes = 6u

static

Constant: The number of base nodes for this element.

Definition at line 48 of file PrismRule6.h.

◆ n_edges

const unsigned MeshLib::PrismRule6::n_edges = 9

static

Constant: The number of edges.

Definition at line 63 of file PrismRule6.h.

◆ n_face_nodes

const unsigned MeshLib::PrismRule6::n_face_nodes = {3, 4, 4, 4, 3}

static

Constant: Table for the number of nodes for each face.

Definition at line 75 of file PrismRule6.h.

Referenced by getFace().

◆ n_faces

const unsigned MeshLib::PrismRule6::n_faces = 5

static

Constant: The number of faces.

Definition at line 60 of file PrismRule6.h.

Referenced by MeshLib::PrismRule15::getFace(), and getFace().

◆ n_neighbors

const unsigned MeshLib::PrismRule6::n_neighbors = 5

static

Constant: The number of neighbors.

Definition at line 66 of file PrismRule6.h.

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

MeshLib/Elements/PrismRule6.h
MeshLib/Elements/PrismRule6.cpp

Detailed Description

Public Types

Static Public Member Functions

Static Public Attributes

Member Typedef Documentation

◆ EdgeReturn

Member Function Documentation

◆ computeVolume()

◆ getFace()

◆ identifyFace()

◆ isPntInElement()

◆ validate()

Member Data Documentation

◆ cell_type

◆ edge_nodes

◆ face_nodes

◆ mesh_elem_type

◆ n_all_nodes

◆ n_base_nodes

◆ n_edges

◆ n_face_nodes

◆ n_faces

◆ n_neighbors