Loading [MathJax]/extensions/MathMenu.js
OGS
HexRule.cpp
Go to the documentation of this file.
1
10#include "HexRule.h"
11
12#include "MathLib/GeometricBasics.h"
13#include "MeshLib/Node.h"
14#include "Quad.h"
15
16namespace MeshLib
17{
18double HexRule::computeVolume(Node const* const* element_nodes)
19{
20 return MathLib::calcTetrahedronVolume(*element_nodes[4],
21 *element_nodes[7],
22 *element_nodes[5],
23 *element_nodes[0]) +
24 MathLib::calcTetrahedronVolume(*element_nodes[5],
25 *element_nodes[3],
26 *element_nodes[1],
27 *element_nodes[0]) +
28 MathLib::calcTetrahedronVolume(*element_nodes[5],
29 *element_nodes[7],
30 *element_nodes[3],
31 *element_nodes[0]) +
32 MathLib::calcTetrahedronVolume(*element_nodes[5],
33 *element_nodes[7],
34 *element_nodes[6],
35 *element_nodes[2]) +
36 MathLib::calcTetrahedronVolume(*element_nodes[1],
37 *element_nodes[3],
38 *element_nodes[5],
39 *element_nodes[2]) +
40 MathLib::calcTetrahedronVolume(*element_nodes[3],
41 *element_nodes[7],
42 *element_nodes[5],
43 *element_nodes[2]);
44}
45
46bool HexRule::isPntInElement(Node const* const* nodes,
47 MathLib::Point3d const& pnt,
48 double eps)
49{
50 return (MathLib::isPointInTetrahedron(
51 pnt, *nodes[4], *nodes[7], *nodes[5], *nodes[0], eps) ||
52 MathLib::isPointInTetrahedron(
53 pnt, *nodes[5], *nodes[3], *nodes[1], *nodes[0], eps) ||
54 MathLib::isPointInTetrahedron(
55 pnt, *nodes[5], *nodes[7], *nodes[3], *nodes[0], eps) ||
56 MathLib::isPointInTetrahedron(
57 pnt, *nodes[5], *nodes[7], *nodes[6], *nodes[2], eps) ||
58 MathLib::isPointInTetrahedron(
59 pnt, *nodes[1], *nodes[3], *nodes[5], *nodes[2], eps) ||
60 MathLib::isPointInTetrahedron(
61 pnt, *nodes[3], *nodes[7], *nodes[5], *nodes[2], eps));
62}
63
64ElementErrorCode HexRule::validate(const Element* e)
65{
66 ElementErrorCode error_code;
67 error_code[ElementErrorFlag::ZeroVolume] = hasZeroVolume(*e);
68
69 for (unsigned i = 0; i < 6; ++i)
70 {
71 if (error_code.all())
72 {
73 break;
74 }
75
76 const MeshLib::Element* quad(e->getFace(i));
77 error_code |= quad->validate();
78 delete quad;
79 }
80 error_code[ElementErrorFlag::NodeOrder] = !e->testElementNodeOrder();
81 return error_code;
82}
83
84} // end namespace MeshLib
Node.h
Definition of the Node class.
Quad.h
Definition of the Quad class.
ElementErrorCode
Collects error flags for mesh elements.
Definition ElementErrorCode.h:34
MeshLib::Element::validate
virtual ElementErrorCode validate() const =0
MeshLib::Element::getFace
virtual const Element * getFace(unsigned i) const =0
Returns the i-th face of the element.
MeshLib::Element::testElementNodeOrder
virtual bool testElementNodeOrder() const =0
MeshLib::HexRule::validate
static ElementErrorCode validate(const Element *e)
Definition HexRule.cpp:64
MeshLib::HexRule::isPntInElement
static bool isPntInElement(Node const *const *nodes, MathLib::Point3d const &pnt, double eps)
Definition HexRule.cpp:46
MeshLib::HexRule::computeVolume
static double computeVolume(Node const *const *element_nodes)
Definition HexRule.cpp:18
MathLib::calcTetrahedronVolume
double calcTetrahedronVolume(MathLib::Point3d const &a, MathLib::Point3d const &b, MathLib::Point3d const &c, MathLib::Point3d const &d)
Definition GeometricBasics.cpp:30
MathLib::isPointInTetrahedron
bool isPointInTetrahedron(MathLib::Point3d const &p, MathLib::Point3d const &a, MathLib::Point3d const &b, MathLib::Point3d const &c, MathLib::Point3d const &d, double eps)
Definition GeometricBasics.cpp:50
MeshLib::hasZeroVolume
bool hasZeroVolume(MeshLib::Element const &element)
Returns true if the element has zero length/area/volume.
Definition Element.cpp:119