OGS
ComputeElementVolumeNumerically.cpp
Go to the documentation of this file.
1
11
12#include "ComputeElementVolumeNumerically.h"
13
14#include <typeinfo>
15
16#include "MeshLib/Elements/Element.h"
17#include "MeshLib/Elements/Hex.h"
18#include "MeshLib/Elements/Line.h"
19#include "MeshLib/Elements/Prism.h"
20#include "MeshLib/Elements/Pyramid.h"
21#include "MeshLib/Elements/Quad.h"
22#include "MeshLib/Elements/Tet.h"
23#include "MeshLib/Elements/Tri.h"
24#include "MeshLib/MeshEnums.h"
25#include "NumLib/Fem/InitShapeMatrices.h"
26#include "NumLib/Fem/Integration/IntegrationMethodRegistry.h"
27#include "NumLib/Fem/ShapeFunction/ShapeHex20.h"
28#include "NumLib/Fem/ShapeFunction/ShapeHex8.h"
29#include "NumLib/Fem/ShapeFunction/ShapeLine2.h"
30#include "NumLib/Fem/ShapeFunction/ShapeLine3.h"
31#include "NumLib/Fem/ShapeFunction/ShapePrism15.h"
32#include "NumLib/Fem/ShapeFunction/ShapePrism6.h"
33#include "NumLib/Fem/ShapeFunction/ShapePyra13.h"
34#include "NumLib/Fem/ShapeFunction/ShapePyra5.h"
35#include "NumLib/Fem/ShapeFunction/ShapeQuad4.h"
36#include "NumLib/Fem/ShapeFunction/ShapeQuad8.h"
37#include "NumLib/Fem/ShapeFunction/ShapeQuad9.h"
38#include "NumLib/Fem/ShapeFunction/ShapeTet10.h"
39#include "NumLib/Fem/ShapeFunction/ShapeTet4.h"
40#include "NumLib/Fem/ShapeFunction/ShapeTri3.h"
41#include "NumLib/Fem/ShapeFunction/ShapeTri6.h"
42#include "NumLib/Fem/ShapeMatrixPolicy.h"
43
44namespace MeshToolsLib
45{
46template <typename ShapeFunction>
47double computeElementVolumeNumerically(MeshLib::Element const& e)
48{
49 // Space dimension is set to 3 in case that 1D or 2D element is inclined.
50 constexpr int space_dim = 3;
51 using ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, space_dim>;
52
53 // Integration order is set to 3:
54 auto const& integration_method =
55 NumLib::IntegrationMethodRegistry::template getIntegrationMethod<
56 typename ShapeFunction::MeshElement>(NumLib::IntegrationOrder{3});
57
58 auto const shape_function_data =
59 NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType, space_dim>(
60 e, false /*is_axially_symmetric*/, integration_method);
61
62 auto const n_integration_points = integration_method.getNumberOfPoints();
63 double volume = 0.0;
64 for (unsigned ip = 0; ip < n_integration_points; ++ip)
65 {
66 auto const weight = integration_method.getWeightedPoint(ip).getWeight();
67 volume += shape_function_data[ip].detJ * weight;
68 }
69
70 return volume;
71}
72
73double computeElementVolumeNumerically(MeshLib::Element const& e)
74{
75 switch (e.getCellType())
76 {
77 case MeshLib::CellType::LINE2:
78 return computeElementVolumeNumerically<NumLib::ShapeLine2>(e);
79 case MeshLib::CellType::LINE3:
80 return computeElementVolumeNumerically<NumLib::ShapeLine3>(e);
81 case MeshLib::CellType::TRI3:
82 return computeElementVolumeNumerically<NumLib::ShapeTri3>(e);
83 case MeshLib::CellType::TRI6:
84 return computeElementVolumeNumerically<NumLib::ShapeTri6>(e);
85 case MeshLib::CellType::QUAD4:
86 return computeElementVolumeNumerically<NumLib::ShapeQuad4>(e);
87 case MeshLib::CellType::QUAD8:
88 return computeElementVolumeNumerically<NumLib::ShapeQuad8>(e);
89 case MeshLib::CellType::QUAD9:
90 return computeElementVolumeNumerically<NumLib::ShapeQuad9>(e);
91 case MeshLib::CellType::TET4:
92 return computeElementVolumeNumerically<NumLib::ShapeTet4>(e);
93 case MeshLib::CellType::HEX8:
94 return computeElementVolumeNumerically<NumLib::ShapeHex8>(e);
95 case MeshLib::CellType::HEX20:
96 return computeElementVolumeNumerically<NumLib::ShapeHex20>(e);
97 case MeshLib::CellType::TET10:
98 return computeElementVolumeNumerically<NumLib::ShapeTet10>(e);
99 case MeshLib::CellType::PRISM6:
100 return computeElementVolumeNumerically<NumLib::ShapePrism6>(e);
101 case MeshLib::CellType::PRISM15:
102 return computeElementVolumeNumerically<NumLib::ShapePrism15>(e);
103 case MeshLib::CellType::PYRAMID5:
104 return computeElementVolumeNumerically<NumLib::ShapePyra5>(e);
105 case MeshLib::CellType::PYRAMID13:
106 return computeElementVolumeNumerically<NumLib::ShapePyra13>(e);
107 default:
108 OGS_FATAL(
109 "Numerical volume calculation is not available for element "
110 "with type {}. ",
111 MeshLib::CellType2String(e.getCellType()));
112 }
113}
114
115} // namespace MeshToolsLib
Element.h
Definition of the Element class.
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:26
Hex.h
Definition of the Hex class.
Line.h
Definition of the Line class.
MeshEnums.h
Definition of mesh-related Enumerations.
Prism.h
Definition of the Prism class.
Pyramid.h
Definition of the Pyramid class.
Quad.h
Definition of the Quad class.
ShapeMatrixPolicyType
EigenFixedShapeMatrixPolicy< ShapeFunction, GlobalDim > ShapeMatrixPolicyType
Definition ShapeMatrixPolicy.h:141
Tet.h
Definition of the Tet class.
Tri.h
Definition of the Tri class.
MeshLib::Element::getCellType
virtual CellType getCellType() const =0
MeshLib::CellType2String
std::string CellType2String(const CellType t)
Given a MeshElemType this returns the appropriate string.
Definition MeshEnums.cpp:157
MeshToolsLib::computeElementVolumeNumerically
double computeElementVolumeNumerically(MeshLib::Element const &e)
Definition ComputeElementVolumeNumerically.cpp:47
NumLib::initShapeMatrices
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
Definition InitShapeMatrices.h:54