OGS
NaturalCoordinatesMapping.cpp
Go to the documentation of this file.
1 
11 #include "NaturalCoordinatesMapping.h"
12 
13 #include <cassert>
14 #ifndef NDEBUG
15 #include <iostream>
16 #endif // NDEBUG
17 
18 #include "BaseLib/Error.h"
19 #include "MeshLib/ElementCoordinatesMappingLocal.h"
20 #include "MeshLib/Elements/HexRule20.h"
21 #include "MeshLib/Elements/HexRule8.h"
22 #include "MeshLib/Elements/LineRule2.h"
23 #include "MeshLib/Elements/LineRule3.h"
24 #include "MeshLib/Elements/PointRule1.h"
25 #include "MeshLib/Elements/PrismRule15.h"
26 #include "MeshLib/Elements/PrismRule6.h"
27 #include "MeshLib/Elements/PyramidRule13.h"
28 #include "MeshLib/Elements/PyramidRule5.h"
29 #include "MeshLib/Elements/QuadRule4.h"
30 #include "MeshLib/Elements/QuadRule8.h"
31 #include "MeshLib/Elements/QuadRule9.h"
32 #include "MeshLib/Elements/TemplateElement.h"
33 #include "MeshLib/Elements/TetRule10.h"
34 #include "MeshLib/Elements/TetRule4.h"
35 #include "MeshLib/Elements/TriRule3.h"
36 #include "MeshLib/Elements/TriRule6.h"
37 #include "NumLib/Fem/ShapeFunction/ShapeHex20.h"
38 #include "NumLib/Fem/ShapeFunction/ShapeHex8.h"
39 #include "NumLib/Fem/ShapeFunction/ShapeLine2.h"
40 #include "NumLib/Fem/ShapeFunction/ShapeLine3.h"
41 #include "NumLib/Fem/ShapeFunction/ShapePoint1.h"
42 #include "NumLib/Fem/ShapeFunction/ShapePrism15.h"
43 #include "NumLib/Fem/ShapeFunction/ShapePrism6.h"
44 #include "NumLib/Fem/ShapeFunction/ShapePyra13.h"
45 #include "NumLib/Fem/ShapeFunction/ShapePyra5.h"
46 #include "NumLib/Fem/ShapeFunction/ShapeQuad4.h"
47 #include "NumLib/Fem/ShapeFunction/ShapeQuad8.h"
48 #include "NumLib/Fem/ShapeFunction/ShapeQuad9.h"
49 #include "NumLib/Fem/ShapeFunction/ShapeTet10.h"
50 #include "NumLib/Fem/ShapeFunction/ShapeTet4.h"
51 #include "NumLib/Fem/ShapeFunction/ShapeTri3.h"
52 #include "NumLib/Fem/ShapeFunction/ShapeTri6.h"
53 #include "NumLib/Fem/ShapeMatrixPolicy.h"
54 #include "ShapeMatrices.h"
55 
56 namespace NumLib
57 {
58 namespace detail
59 {
60 template <ShapeMatrixType FIELD_TYPE>
61 struct FieldType
62 {
63 };
64 
65 template <class T_MESH_ELEMENT, class T_SHAPE_FUNC, class T_SHAPE_MATRICES>
66 inline void computeMappingMatrices(
67  const T_MESH_ELEMENT& /*ele*/,
68  const double* natural_pt,
69  const MeshLib::ElementCoordinatesMappingLocal& /*ele_local_coord*/,
70  T_SHAPE_MATRICES& shapemat,
71  FieldType<ShapeMatrixType::N> /*unused*/)
72 {
73  T_SHAPE_FUNC::computeShapeFunction(natural_pt, shapemat.N);
74 }
75 
76 template <class T_MESH_ELEMENT, class T_SHAPE_FUNC, class T_SHAPE_MATRICES>
77 inline void computeMappingMatrices(
78  const T_MESH_ELEMENT& /*ele*/,
79  [[maybe_unused]] const double* natural_pt,
80  const MeshLib::ElementCoordinatesMappingLocal& /*ele_local_coord*/,
81  [[maybe_unused]] T_SHAPE_MATRICES& shapemat,
82  FieldType<ShapeMatrixType::DNDR> /*unused*/)
83 {
84  if constexpr (T_SHAPE_FUNC::DIM != 0)
85  {
86  double* const dNdr = shapemat.dNdr.data();
87  T_SHAPE_FUNC::computeGradShapeFunction(natural_pt, dNdr);
88  }
89 }
90 
91 static void checkJacobianDeterminant(const double detJ,
92  MeshLib::Element const& element)
93 {
94  if (detJ > 0)
95  { // The usual case
96  return;
97  }
98 
99  if (detJ < 0)
100  {
101  ERR("det J = {:g} is negative for element {:d}.",
102  detJ,
103  element.getID());
104 #ifndef NDEBUG
105  std::cerr << element << "\n";
106 #endif // NDEBUG
107  OGS_FATAL(
108  "Please check whether the node numbering of the element is correct,"
109  "or additional elements (like boundary elements) are still present "
110  "in the mesh.");
111  }
112 
113  if (detJ == 0)
114  {
115  ERR("det J is zero for element {:d}.", element.getID());
116 #ifndef NDEBUG
117  std::cerr << element << "\n";
118 #endif // NDEBUG
119  OGS_FATAL(
120  "Please check whether:\n"
121  "\t the element nodes may have the same coordinates,\n"
122  "\t or the coordinates of all nodes are not given on the x-axis "
123  "for a 1D problem or in the xy-plane in the 2D case.\n"
124  "The first case can occur, if not all boundary elements"
125  "were removed from the bulk mesh.");
126  }
127 }
128 
129 template <class T_MESH_ELEMENT, class T_SHAPE_FUNC, class T_SHAPE_MATRICES>
130 inline void computeMappingMatrices(
131  const T_MESH_ELEMENT& ele,
132  [[maybe_unused]] const double* natural_pt,
133  const MeshLib::ElementCoordinatesMappingLocal& ele_local_coord,
134  T_SHAPE_MATRICES& shapemat,
135  FieldType<ShapeMatrixType::DNDR_J> /*unused*/)
136 {
137  if constexpr (T_SHAPE_FUNC::DIM != 0)
138  {
139  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
140  ele,
141  natural_pt,
142  ele_local_coord,
143  shapemat,
144  FieldType<ShapeMatrixType::DNDR>());
145 
146  auto const dim = T_MESH_ELEMENT::dimension;
147  auto const nnodes = T_MESH_ELEMENT::n_all_nodes;
148 
149  // jacobian: J=[dx/dr dy/dr // dx/ds dy/ds]
150  for (auto k = decltype(nnodes){0}; k < nnodes; k++)
151  {
152  const MathLib::Point3d& mapped_pt =
153  ele_local_coord.getMappedCoordinates(k);
154  // outer product of dNdr and mapped_pt for a particular node
155  for (auto i_r = decltype(dim){0}; i_r < dim; i_r++)
156  {
157  for (auto j_x = decltype(dim){0}; j_x < dim; j_x++)
158  {
159  shapemat.J(i_r, j_x) +=
160  shapemat.dNdr(i_r, k) * mapped_pt[j_x];
161  }
162  }
163  }
164 
165  shapemat.detJ = shapemat.J.determinant();
166  checkJacobianDeterminant(shapemat.detJ, ele);
167  }
168  else
169  {
170  shapemat.detJ = 1.0;
171  }
172 }
173 
174 template <class T_MESH_ELEMENT, class T_SHAPE_FUNC, class T_SHAPE_MATRICES>
175 inline void computeMappingMatrices(
176  const T_MESH_ELEMENT& ele,
177  const double* natural_pt,
178  const MeshLib::ElementCoordinatesMappingLocal& ele_local_coord,
179  T_SHAPE_MATRICES& shapemat,
180  FieldType<ShapeMatrixType::N_J> /*unused*/)
181 {
182  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
183  ele,
184  natural_pt,
185  ele_local_coord,
186  shapemat,
187  FieldType<ShapeMatrixType::N>());
188  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
189  ele,
190  natural_pt,
191  ele_local_coord,
192  shapemat,
193  FieldType<ShapeMatrixType::DNDR_J>());
194 }
195 
196 template <class T_MESH_ELEMENT, class T_SHAPE_FUNC, class T_SHAPE_MATRICES>
197 inline void computeMappingMatrices(
198  const T_MESH_ELEMENT& ele,
199  const double* natural_pt,
200  const MeshLib::ElementCoordinatesMappingLocal& ele_local_coord,
201  T_SHAPE_MATRICES& shapemat,
202  FieldType<ShapeMatrixType::DNDX> /*unused*/)
203 {
204  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
205  ele,
206  natural_pt,
207  ele_local_coord,
208  shapemat,
209  FieldType<ShapeMatrixType::DNDR_J>());
210  if constexpr (T_SHAPE_FUNC::DIM != 0)
211  {
212  checkJacobianDeterminant(shapemat.detJ, ele);
213 
214  // J^-1, dshape/dx
215  shapemat.invJ.noalias() = shapemat.J.inverse();
216 
217  assert(shapemat.dNdr.rows() == ele.getDimension());
218  shapemat.dNdx
219  .template topLeftCorner<T_SHAPE_FUNC::DIM, T_SHAPE_FUNC::NPOINTS>()
220  .noalias() = shapemat.invJ * shapemat.dNdr;
221  }
222 }
223 
224 template <class T_MESH_ELEMENT, class T_SHAPE_FUNC, class T_SHAPE_MATRICES>
225 inline void computeMappingMatrices(
226  const T_MESH_ELEMENT& ele,
227  const double* natural_pt,
228  const MeshLib::ElementCoordinatesMappingLocal& ele_local_coord,
229  T_SHAPE_MATRICES& shapemat,
230  FieldType<ShapeMatrixType::ALL> /*unused*/)
231 {
232  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
233  ele,
234  natural_pt,
235  ele_local_coord,
236  shapemat,
237  FieldType<ShapeMatrixType::N>());
238  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
239  ele,
240  natural_pt,
241  ele_local_coord,
242  shapemat,
243  FieldType<ShapeMatrixType::DNDX>());
244 }
245 
246 template <class T_MESH_ELEMENT,
247  class T_SHAPE_FUNC,
248  class T_SHAPE_MATRICES,
249  ShapeMatrixType T_SHAPE_MATRIX_TYPE>
250 void naturalCoordinatesMappingComputeShapeMatrices(const T_MESH_ELEMENT& ele,
251  const double* natural_pt,
252  T_SHAPE_MATRICES& shapemat,
253  const unsigned global_dim)
254 {
255  const MeshLib::ElementCoordinatesMappingLocal ele_local_coord(ele,
256  global_dim);
257 
258  detail::
259  computeMappingMatrices<T_MESH_ELEMENT, T_SHAPE_FUNC, T_SHAPE_MATRICES>(
260  ele,
261  natural_pt,
262  ele_local_coord,
263  shapemat,
264  detail::FieldType<T_SHAPE_MATRIX_TYPE>());
265 }
266 
267 #define OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
268  RULE, SHAPE, DIM, WHICHPART, SHAPEMATRIXPOLICY) \
269  template void naturalCoordinatesMappingComputeShapeMatrices< \
270  MeshLib::TemplateElement<MeshLib::RULE>, \
271  NumLib::SHAPE, \
272  SHAPEMATRIXPOLICY<NumLib::SHAPE, DIM>::ShapeMatrices, \
273  ShapeMatrixType::WHICHPART>( \
274  MeshLib::TemplateElement<MeshLib::RULE> const&, \
275  double const*, \
276  SHAPEMATRIXPOLICY<NumLib::SHAPE, DIM>::ShapeMatrices&, \
277  const unsigned global_dim)
278 
279 #define OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(RULE, SHAPE) \
280  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
281  RULE, SHAPE, 0, ALL, EigenDynamicShapeMatrixPolicy); \
282  /* Those instantiations are needed in unit tests only */ \
283  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
284  RULE, SHAPE, 0, N, EigenDynamicShapeMatrixPolicy); \
285  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
286  RULE, SHAPE, 0, DNDR, EigenDynamicShapeMatrixPolicy); \
287  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
288  RULE, SHAPE, 0, N_J, EigenDynamicShapeMatrixPolicy); \
289  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
290  RULE, SHAPE, 0, DNDR_J, EigenDynamicShapeMatrixPolicy); \
291  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
292  RULE, SHAPE, 0, DNDX, EigenDynamicShapeMatrixPolicy)
293 
294 #define OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(RULE, SHAPE, DIM) \
295  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
296  RULE, SHAPE, DIM, ALL, EigenFixedShapeMatrixPolicy); \
297  /* Those instantiations are needed in unit tests only */ \
298  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
299  RULE, SHAPE, DIM, N, EigenFixedShapeMatrixPolicy); \
300  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
301  RULE, SHAPE, DIM, DNDR, EigenFixedShapeMatrixPolicy); \
302  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
303  RULE, SHAPE, DIM, N_J, EigenFixedShapeMatrixPolicy); \
304  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
305  RULE, SHAPE, DIM, DNDR_J, EigenFixedShapeMatrixPolicy); \
306  OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_PART( \
307  RULE, SHAPE, DIM, DNDX, EigenFixedShapeMatrixPolicy)
308 
309 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(HexRule20, ShapeHex20);
310 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(HexRule8, ShapeHex8);
311 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(LineRule2, ShapeLine2);
312 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(LineRule3, ShapeLine3);
313 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(PointRule1, ShapePoint1);
314 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(PrismRule15, ShapePrism15);
315 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(PrismRule6, ShapePrism6);
316 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(PyramidRule13, ShapePyra13);
317 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(PyramidRule5, ShapePyra5);
318 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(QuadRule4, ShapeQuad4);
319 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(QuadRule8, ShapeQuad8);
320 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(QuadRule9, ShapeQuad9);
321 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(TetRule10, ShapeTet10);
322 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(TetRule4, ShapeTet4);
323 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(TriRule3, ShapeTri3);
324 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(TriRule6, ShapeTri6);
325 
326 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule20, ShapeHex20, 1);
327 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule8, ShapeHex8, 1);
328 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(LineRule2, ShapeLine2, 1);
329 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(LineRule3, ShapeLine3, 1);
330 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PointRule1, ShapePoint1, 1);
331 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PrismRule15, ShapePrism15, 1);
332 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PrismRule6, ShapePrism6, 1);
333 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PyramidRule13, ShapePyra13, 1);
334 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PyramidRule5, ShapePyra5, 1);
335 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule4, ShapeQuad4, 1);
336 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule8, ShapeQuad8, 1);
337 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule9, ShapeQuad9, 1);
338 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TetRule10, ShapeTet10, 1);
339 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TetRule4, ShapeTet4, 1);
340 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TriRule3, ShapeTri3, 1);
341 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TriRule6, ShapeTri6, 1);
342 
343 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule20, ShapeHex20, 2);
344 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule8, ShapeHex8, 2);
345 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(LineRule2, ShapeLine2, 2);
346 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(LineRule3, ShapeLine3, 2);
347 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PointRule1, ShapePoint1, 2);
348 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PrismRule15, ShapePrism15, 2);
349 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PrismRule6, ShapePrism6, 2);
350 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PyramidRule13, ShapePyra13, 2);
351 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PyramidRule5, ShapePyra5, 2);
352 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule4, ShapeQuad4, 2);
353 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule8, ShapeQuad8, 2);
354 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule9, ShapeQuad9, 2);
355 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TetRule10, ShapeTet10, 2);
356 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TetRule4, ShapeTet4, 2);
357 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TriRule3, ShapeTri3, 2);
358 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TriRule6, ShapeTri6, 2);
359 
360 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule20, ShapeHex20, 3);
361 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule8, ShapeHex8, 3);
362 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(LineRule2, ShapeLine2, 3);
363 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(LineRule3, ShapeLine3, 3);
364 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PointRule1, ShapePoint1, 3);
365 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PrismRule15, ShapePrism15, 3);
366 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PrismRule6, ShapePrism6, 3);
367 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PyramidRule13, ShapePyra13, 3);
368 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(PyramidRule5, ShapePyra5, 3);
369 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule4, ShapeQuad4, 3);
370 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule8, ShapeQuad8, 3);
371 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(QuadRule9, ShapeQuad9, 3);
372 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TetRule10, ShapeTet10, 3);
373 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TetRule4, ShapeTet4, 3);
374 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TriRule3, ShapeTri3, 3);
375 OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(TriRule6, ShapeTri6, 3);
376 
377 } // namespace detail
378 
379 } // namespace NumLib
OGS_FATAL
#define OGS_FATAL(...)
Definition: Error.h:26
ERR
void ERR(char const *fmt, Args const &... args)
Definition: Logging.h:42
MeshLib::ElementCoordinatesMappingLocal::getMappedCoordinates
MathLib::Point3d const & getMappedCoordinates(std::size_t node_id) const
return mapped coordinates of the node
Definition: ElementCoordinatesMappingLocal.h:42
MeshLib::Element::getID
virtual std::size_t getID() const final
Returns the ID of the element.
Definition: Element.h:82
NumLib::ShapePoint1
Shape function for a point element in natural coordinates.
Definition: ShapePoint1.h:19
NumLib::detail::naturalCoordinatesMappingComputeShapeMatrices
void naturalCoordinatesMappingComputeShapeMatrices(const T_MESH_ELEMENT &ele, const double *natural_pt, T_SHAPE_MATRICES &shapemat, const unsigned global_dim)
Used to explicitly instantiate the NaturalCoordinatesMapping class template.
Definition: NaturalCoordinatesMapping.cpp:250
NumLib::detail::computeMappingMatrices
void computeMappingMatrices(const T_MESH_ELEMENT &, const double *natural_pt, const MeshLib::ElementCoordinatesMappingLocal &, T_SHAPE_MATRICES &shapemat, FieldType< ShapeMatrixType::N >)
Definition: NaturalCoordinatesMapping.cpp:66
NumLib::detail::checkJacobianDeterminant
static void checkJacobianDeterminant(const double detJ, MeshLib::Element const &element)
Definition: NaturalCoordinatesMapping.cpp:91
NumLib::detail::OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX
OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_FIX(HexRule20, ShapeHex20, 1)
NumLib::detail::OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN
OGS_INSTANTIATE_NATURAL_COORDINATES_MAPPING_DYN(HexRule20, ShapeHex20)
NumLib::ShapeMatrixType
ShapeMatrixType
Shape matrix type to be calculated.
Definition: ShapeMatrices.h:24