From 4d4d197240e87645d4fbafcef727d6567eb2afb1 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 27 Nov 2017 18:06:50 +0100
Subject: [PATCH] [NL] Add local coordinates for each element type.
---
.../NaturalNodeCoordinates.cpp | 58 ++++
.../NaturalNodeCoordinates.h | 254 ++++++++++++++++++
2 files changed, 312 insertions(+)
create mode 100644 NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.cpp
create mode 100644 NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.h
diff --git a/NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.cpp b/NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.cpp
new file mode 100644
index 00000000000..8f8bb033887
--- /dev/null
+++ b/NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.cpp
@@ -0,0 +1,58 @@
+/**
+ * \copyright
+ * Copyright (c) 2012-2018, OpenGeoSys Community (http://www.opengeosys.org)
+ * Distributed under a Modified BSD License.
+ * See accompanying file LICENSE.txt or
+ * http://www.opengeosys.org/project/license
+ *
+ */
+
+#include "NaturalNodeCoordinates.h"
+
+namespace NumLib
+{
+constexpr std::array<std::array<double, 3>, 2>
+ NaturalCoordinates<MeshLib::Line>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 3>
+ NaturalCoordinates<MeshLib::Line3>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 3>
+ NaturalCoordinates<MeshLib::Tri>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 6>
+ NaturalCoordinates<MeshLib::Tri6>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 4>
+ NaturalCoordinates<MeshLib::Quad>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 8>
+ NaturalCoordinates<MeshLib::Quad8>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 9>
+ NaturalCoordinates<MeshLib::Quad9>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 4>
+ NaturalCoordinates<MeshLib::Tet>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 10>
+ NaturalCoordinates<MeshLib::Tet10>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 6>
+ NaturalCoordinates<MeshLib::Prism>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 15>
+ NaturalCoordinates<MeshLib::Prism15>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 5>
+ NaturalCoordinates<MeshLib::Pyramid>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 13>
+ NaturalCoordinates<MeshLib::Pyramid13>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 8>
+ NaturalCoordinates<MeshLib::Hex>::coordinates;
+
+constexpr std::array<std::array<double, 3>, 20>
+ NaturalCoordinates<MeshLib::Hex20>::coordinates;
+} // namespace NumLib
diff --git a/NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.h b/NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.h
new file mode 100644
index 00000000000..a55c968f9a2
--- /dev/null
+++ b/NumLib/Fem/CoordinatesMapping/NaturalNodeCoordinates.h
@@ -0,0 +1,254 @@
+/**
+ * \copyright
+ * Copyright (c) 2012-2018, OpenGeoSys Community (http://www.opengeosys.org)
+ * Distributed under a Modified BSD License.
+ * See accompanying file LICENSE.txt or
+ * http://www.opengeosys.org/project/license
+ *
+ */
+
+#include <cassert>
+
+#include "BaseLib/Error.h"
+#include "MathLib/Point3d.h"
+
+#include "MeshLib/ElementCoordinatesMappingLocal.h"
+#include "MeshLib/Elements/HexRule20.h"
+#include "MeshLib/Elements/HexRule8.h"
+#include "MeshLib/Elements/LineRule2.h"
+#include "MeshLib/Elements/LineRule3.h"
+#include "MeshLib/Elements/PointRule1.h"
+#include "MeshLib/Elements/PrismRule15.h"
+#include "MeshLib/Elements/PrismRule6.h"
+#include "MeshLib/Elements/PyramidRule13.h"
+#include "MeshLib/Elements/PyramidRule5.h"
+#include "MeshLib/Elements/QuadRule4.h"
+#include "MeshLib/Elements/QuadRule8.h"
+#include "MeshLib/Elements/QuadRule9.h"
+#include "MeshLib/Elements/TemplateElement.h"
+#include "MeshLib/Elements/TetRule10.h"
+#include "MeshLib/Elements/TetRule4.h"
+#include "MeshLib/Elements/TriRule3.h"
+#include "MeshLib/Elements/TriRule6.h"
+
+#include "NumLib/Fem/ShapeFunction/ShapeHex20.h"
+#include "NumLib/Fem/ShapeFunction/ShapeHex8.h"
+#include "NumLib/Fem/ShapeFunction/ShapeLine2.h"
+#include "NumLib/Fem/ShapeFunction/ShapeLine3.h"
+#include "NumLib/Fem/ShapeFunction/ShapePoint1.h"
+#include "NumLib/Fem/ShapeFunction/ShapePrism15.h"
+#include "NumLib/Fem/ShapeFunction/ShapePrism6.h"
+#include "NumLib/Fem/ShapeFunction/ShapePyra13.h"
+#include "NumLib/Fem/ShapeFunction/ShapePyra5.h"
+#include "NumLib/Fem/ShapeFunction/ShapeQuad4.h"
+#include "NumLib/Fem/ShapeFunction/ShapeQuad8.h"
+#include "NumLib/Fem/ShapeFunction/ShapeQuad9.h"
+#include "NumLib/Fem/ShapeFunction/ShapeTet10.h"
+#include "NumLib/Fem/ShapeFunction/ShapeTet4.h"
+#include "NumLib/Fem/ShapeFunction/ShapeTri3.h"
+#include "NumLib/Fem/ShapeFunction/ShapeTri6.h"
+#include "NumLib/Fem/ShapeMatrixPolicy.h"
+
+#include "ShapeMatrices.h"
+
+namespace NumLib
+{
+/// A map 3D natural coordinates (-1 to 1, or 0 to 1) for each node of a
+/// MeshLib::Element.
+template <typename MeshLibElement>
+struct NaturalCoordinates;
+
+template <>
+struct NaturalCoordinates<MeshLib::Line>
+{
+ static constexpr std::array<std::array<double, 3>, 2> coordinates = {
+ {{{-1, 0, 0}}, {{1, 0, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Line3>
+{
+ static constexpr std::array<std::array<double, 3>, 3> coordinates = {
+ {{{-1, 0, 0}}, {{1, 0, 0}}, {{0, 0, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Tri>
+{
+ static constexpr std::array<std::array<double, 3>, 3> coordinates = {
+ {{{0, 0, 0}}, {{1, 0, 0}}, {{0, 1, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Tri6>
+{
+ static constexpr std::array<std::array<double, 3>, 6> coordinates = {
+ {{{0, 0, 0}},
+ {{1, 0, 0}},
+ {{0, 1, 0}},
+ {{0.5, 0, 0}},
+ {{0.5, 0.5, 0}},
+ {{0, 0.5, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Quad>
+{
+ static constexpr std::array<std::array<double, 3>, 4> coordinates = {
+ {{{1, 1, 0}}, {{-1, 1, 0}}, {{-1, -1, 0}}, {{1, -1, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Quad8>
+{
+ static constexpr std::array<std::array<double, 3>, 8> coordinates = {
+ {{{1, 1, 0}},
+ {{-1, 1, 0}},
+ {{-1, -1, 0}},
+ {{1, -1, 0}},
+ {{0, 1, 0}},
+ {{-1, 0, 0}},
+ {{0, -1, 0}},
+ {{1, 0, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Quad9>
+{
+ static constexpr std::array<std::array<double, 3>, 9> coordinates = {
+ {{{1, 1, 0}},
+ {{-1, 1, 0}},
+ {{-1, -1, 0}},
+ {{1, -1, 0}},
+ {{0, 1, 0}},
+ {{-1, 0, 0}},
+ {{0, -1, 0}},
+ {{1, 0, 0}},
+ {{0, 0, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Tet>
+{
+ static constexpr std::array<std::array<double, 3>, 4> coordinates = {
+ {{{0, 0, 0}}, {{1, 0, 0}}, {{0, 1, 0}}, {{0, 0, 1}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Tet10>
+{
+ static constexpr std::array<std::array<double, 3>, 10> coordinates = {
+ {{{0, 0, 0}},
+ {{1, 0, 0}},
+ {{0, 1, 0}},
+ {{0, 0, 1}},
+
+ {{0.5, 0, 0}},
+ {{0.5, 0.5, 0}},
+ {{0, 0.5, 0}},
+
+ {{0, 0, 0.5}},
+ {{0.5, 0, 0.5}},
+ {{0, 0.5, 0.5}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Prism>
+{
+ static constexpr std::array<std::array<double, 3>, 6> coordinates = {
+ {{{0, 0, -1}},
+ {{1, 0, -1}},
+ {{0, 1, -1}},
+ {{0, 0, 1}},
+ {{1, 0, 1}},
+ {{0, 1, 1}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Prism15>
+{
+ static constexpr std::array<std::array<double, 3>, 15> coordinates = {
+ {{{0, 0, -1}},
+ {{1, 0, -1}},
+ {{0, 1, -1}},
+ {{0, 0, 1}},
+ {{1, 0, 1}},
+ {{0, 1, 1}},
+
+ {{0.5, 0, -1}},
+ {{0.5, 0.5, -1}},
+ {{0, 0.5, -1}},
+
+ {{0.5, 0, 1}},
+ {{0.5, 0.5, 1}},
+ {{0, 0.5, 1}},
+
+ {{0, 0, 0}},
+ {{1, 0, 0}},
+ {{0, 1, 0}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Pyramid>
+{
+ static constexpr std::array<std::array<double, 3>, 5> coordinates = {
+ {{{-1, -1, -1}},
+ {{1, -1, -1}},
+ {{1, 1, -1}},
+ {{-1, 1, -1}},
+ {{0, 0, 1}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Pyramid13>
+{
+ static constexpr std::array<std::array<double, 3>, 13> coordinates = {{
+ {{-1, -1, -1}},
+ {{1, -1, -1}},
+ {{1, 1, -1}},
+ {{-1, 1, -1}},
+
+ {{0, 0, 1}},
+
+ {{0, -1, -1}},
+ {{1, 0, -1}},
+ {{0, 1, -1}},
+ {{-1, 0, -1}},
+
+ {{-1, -1, 0}},
+ {{1, -1, 0}},
+ {{1, 1, 0}},
+ {{-1, 1, 0}},
+ }};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Hex>
+{
+ static constexpr std::array<std::array<double, 3>,
+ MeshLib::Hex::n_all_nodes>
+ coordinates = {{{{-1, -1, -1}},
+ {{1, -1, -1}},
+ {{1, 1, -1}},
+ {{-1, 1, -1}},
+ {{-1, -1, 1}},
+ {{1, -1, 1}},
+ {{1, 1, 1}},
+ {{-1, 1, 1}}}};
+};
+
+template <>
+struct NaturalCoordinates<MeshLib::Hex20>
+{
+ static constexpr std::array<std::array<double, 3>,
+ MeshLib::Hex20::n_all_nodes>
+ coordinates = {
+ {{{-1, -1, -1}}, {{1, -1, -1}}, {{1, 1, -1}}, {{-1, 1, -1}},
+ {{-1, -1, 1}}, {{1, -1, 1}}, {{1, 1, 1}}, {{-1, 1, 1}},
+
+ {{0, -1, -1}}, {{1, 0, -1}}, {{0, 1, -1}}, {{-1, 0, -1}},
+ {{0, -1, 1}}, {{1, 0, 1}}, {{0, 1, 1}}, {{-1, 0, 1}},
+
+ {{-1, -1, 0}}, {{1, -1, 0}}, {{1, 1, 0}}, {{-1, 1, 0}}}};
+};
+} // namespace NumLib
--
GitLab