diff --git a/test/basic/CMakeLists.txt b/test/basic/CMakeLists.txt
index 66c4b33..725bba6 100644
--- a/test/basic/CMakeLists.txt
+++ b/test/basic/CMakeLists.txt
@@ -1,53 +1,54 @@
cmake_minimum_required(VERSION 3.0)
if(CMAKE_CURRENT_SOURCE_DIR STREQUAL CMAKE_SOURCE_DIR)
project(GooseFEM-test-basic)
find_package(GooseFEM REQUIRED CONFIG)
endif()
set(ASSERT ON)
set(DEBUG ON)
option(XSIMD "Use xsimd optimisations" ON)
set(test_name "test-basic")
find_package(Catch2 REQUIRED)
find_package(xtensor REQUIRED)
add_executable(${test_name}
main.cpp
version.cpp
Allocate.cpp
Element.cpp
ElementHex8.cpp
ElementQuad4.cpp
ElementQuad4Planar.cpp
Iterate.cpp
Matrix.cpp
MatrixDiagonal.cpp
Mesh.cpp
MeshQuad4.cpp
+ MeshHex8.cpp
Vector.cpp
VectorPartitioned.cpp)
target_link_libraries(${test_name} PRIVATE
Catch2::Catch2
GooseFEM
GooseFEM::compiler_warnings
xtensor::optimize)
if(ASSERT)
target_link_libraries(${test_name} PRIVATE GooseFEM::assert)
endif()
if(DEBUG)
target_link_libraries(${test_name} PRIVATE GooseFEM::debug)
endif()
if(SIMD)
find_package(xsimd REQUIRED)
target_link_libraries(${test_name} PRIVATE xtensor::use_xsimd)
endif()
add_test(NAME ${test_name} COMMAND ${test_name})
diff --git a/test/basic/MeshHex8.cpp b/test/basic/MeshHex8.cpp
new file mode 100644
index 0000000..9a97136
--- /dev/null
+++ b/test/basic/MeshHex8.cpp
@@ -0,0 +1,763 @@
+#include <catch2/catch.hpp>
+#include <GooseFEM/MeshHex8.h>
+
+TEST_CASE("GooseFEM::MeshHex8", "MeshHex8.h")
+{
+ SECTION("Regular")
+ {
+ xt::xtensor<double, 2> coor_ =
+ {{0., 0., 0.},
+ {1., 0., 0.},
+ {2., 0., 0.},
+ {3., 0., 0.},
+ {0., 1., 0.},
+ {1., 1., 0.},
+ {2., 1., 0.},
+ {3., 1., 0.},
+ {0., 2., 0.},
+ {1., 2., 0.},
+ {2., 2., 0.},
+ {3., 2., 0.},
+ {0., 3., 0.},
+ {1., 3., 0.},
+ {2., 3., 0.},
+ {3., 3., 0.},
+ {0., 0., 1.},
+ {1., 0., 1.},
+ {2., 0., 1.},
+ {3., 0., 1.},
+ {0., 1., 1.},
+ {1., 1., 1.},
+ {2., 1., 1.},
+ {3., 1., 1.},
+ {0., 2., 1.},
+ {1., 2., 1.},
+ {2., 2., 1.},
+ {3., 2., 1.},
+ {0., 3., 1.},
+ {1., 3., 1.},
+ {2., 3., 1.},
+ {3., 3., 1.},
+ {0., 0., 2.},
+ {1., 0., 2.},
+ {2., 0., 2.},
+ {3., 0., 2.},
+ {0., 1., 2.},
+ {1., 1., 2.},
+ {2., 1., 2.},
+ {3., 1., 2.},
+ {0., 2., 2.},
+ {1., 2., 2.},
+ {2., 2., 2.},
+ {3., 2., 2.},
+ {0., 3., 2.},
+ {1., 3., 2.},
+ {2., 3., 2.},
+ {3., 3., 2.},
+ {0., 0., 3.},
+ {1., 0., 3.},
+ {2., 0., 3.},
+ {3., 0., 3.},
+ {0., 1., 3.},
+ {1., 1., 3.},
+ {2., 1., 3.},
+ {3., 1., 3.},
+ {0., 2., 3.},
+ {1., 2., 3.},
+ {2., 2., 3.},
+ {3., 2., 3.},
+ {0., 3., 3.},
+ {1., 3., 3.},
+ {2., 3., 3.},
+ {3., 3., 3.}};
+
+ xt::xtensor<double, 2> conn_ =
+ {{ 0, 1, 5, 4, 16, 17, 21, 20},
+ { 1, 2, 6, 5, 17, 18, 22, 21},
+ { 2, 3, 7, 6, 18, 19, 23, 22},
+ { 4, 5, 9, 8, 20, 21, 25, 24},
+ { 5, 6, 10, 9, 21, 22, 26, 25},
+ { 6, 7, 11, 10, 22, 23, 27, 26},
+ { 8, 9, 13, 12, 24, 25, 29, 28},
+ { 9, 10, 14, 13, 25, 26, 30, 29},
+ {10, 11, 15, 14, 26, 27, 31, 30},
+ {16, 17, 21, 20, 32, 33, 37, 36},
+ {17, 18, 22, 21, 33, 34, 38, 37},
+ {18, 19, 23, 22, 34, 35, 39, 38},
+ {20, 21, 25, 24, 36, 37, 41, 40},
+ {21, 22, 26, 25, 37, 38, 42, 41},
+ {22, 23, 27, 26, 38, 39, 43, 42},
+ {24, 25, 29, 28, 40, 41, 45, 44},
+ {25, 26, 30, 29, 41, 42, 46, 45},
+ {26, 27, 31, 30, 42, 43, 47, 46},
+ {32, 33, 37, 36, 48, 49, 53, 52},
+ {33, 34, 38, 37, 49, 50, 54, 53},
+ {34, 35, 39, 38, 50, 51, 55, 54},
+ {36, 37, 41, 40, 52, 53, 57, 56},
+ {37, 38, 42, 41, 53, 54, 58, 57},
+ {38, 39, 43, 42, 54, 55, 59, 58},
+ {40, 41, 45, 44, 56, 57, 61, 60},
+ {41, 42, 46, 45, 57, 58, 62, 61},
+ {42, 43, 47, 46, 58, 59, 63, 62}};
+
+ xt::xtensor<size_t, 1> nodesBottom_ =
+ { 0, 1, 2, 3, 16, 17, 18, 19, 32, 33, 34, 35, 48, 49, 50, 51};
+
+ xt::xtensor<size_t, 1> nodesTop_ =
+ {12, 13, 14, 15, 28, 29, 30, 31, 44, 45, 46, 47, 60, 61, 62, 63};
+
+ xt::xtensor<size_t, 1> nodesLeft_ =
+ { 0, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60};
+
+ xt::xtensor<size_t, 1> nodesRight_ =
+ { 3, 7, 11, 15, 19, 23, 27, 31, 35, 39, 43, 47, 51, 55, 59, 63};
+
+ xt::xtensor<size_t, 1> nodesFront_ =
+ { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15};
+
+ xt::xtensor<size_t, 1> nodesBack_ =
+ {48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63};
+
+ xt::xtensor<size_t, 1> nodesFrontTopEdge_ = {12, 13, 14, 15};
+ xt::xtensor<size_t, 1> nodesFrontLeftEdge_ = { 0, 4, 8, 12};
+ xt::xtensor<size_t, 1> nodesFrontRightEdge_ = { 3, 7, 11, 15};
+ xt::xtensor<size_t, 1> nodesBackBottomEdge_ = {48, 49, 50, 51};
+ xt::xtensor<size_t, 1> nodesBackTopEdge_ = {60, 61, 62, 63};
+ xt::xtensor<size_t, 1> nodesBackLeftEdge_ = {48, 52, 56, 60};
+ xt::xtensor<size_t, 1> nodesBackRightEdge_ = {51, 55, 59, 63};
+ xt::xtensor<size_t, 1> nodesBottomLeftEdge_ = { 0, 16, 32, 48};
+ xt::xtensor<size_t, 1> nodesBottomRightEdge_ = { 3, 19, 35, 51};
+ xt::xtensor<size_t, 1> nodesTopLeftEdge_ = {12, 28, 44, 60};
+ xt::xtensor<size_t, 1> nodesTopRightEdge_ = {15, 31, 47, 63};
+ xt::xtensor<size_t, 1> nodesFrontFace_ = { 5, 6, 9, 10};
+ xt::xtensor<size_t, 1> nodesBackFace_ = {53, 54, 57, 58};
+ xt::xtensor<size_t, 1> nodesLeftFace_ = {20, 24, 36, 40};
+ xt::xtensor<size_t, 1> nodesRightFace_ = {23, 27, 39, 43};
+ xt::xtensor<size_t, 1> nodesBottomFace_ = {17, 18, 33, 34};
+ xt::xtensor<size_t, 1> nodesTopFace_ = {29, 30, 45, 46};
+ xt::xtensor<size_t, 1> nodesFrontBottomOpenEdge_ = {1, 2};
+ xt::xtensor<size_t, 1> nodesFrontTopOpenEdge_ = {13, 14};
+ xt::xtensor<size_t, 1> nodesFrontLeftOpenEdge_ = {4, 8};
+ xt::xtensor<size_t, 1> nodesFrontRightOpenEdge_ = { 7, 11};
+ xt::xtensor<size_t, 1> nodesBackBottomOpenEdge_ = {49, 50};
+ xt::xtensor<size_t, 1> nodesBackTopOpenEdge_ = {61, 62};
+ xt::xtensor<size_t, 1> nodesBackLeftOpenEdge_ = {52, 56};
+ xt::xtensor<size_t, 1> nodesBackRightOpenEdge_ = {55, 59};
+ xt::xtensor<size_t, 1> nodesBottomLeftOpenEdge_ = {16, 32};
+ xt::xtensor<size_t, 1> nodesBottomRightOpenEdge_ = {19, 35};
+ xt::xtensor<size_t, 1> nodesTopLeftOpenEdge_ = {28, 44};
+ xt::xtensor<size_t, 1> nodesTopRightOpenEdge_ = {31, 47};
+ size_t nodesFrontBottomLeftCorner_ = 0;
+ size_t nodesFrontBottomRightCorner_ = 3;
+ size_t nodesFrontTopLeftCorner_ = 12;
+ size_t nodesFrontTopRightCorner_ = 15;
+ size_t nodesBackBottomLeftCorner_ = 48;
+ size_t nodesBackBottomRightCorner_ = 51;
+ size_t nodesBackTopLeftCorner_ = 60;
+ size_t nodesBackTopRightCorner_ = 63;
+
+ xt::xtensor<size_t, 2> dofs_ =
+ {{ 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 9, 10, 11},
+ { 12, 13, 14},
+ { 15, 16, 17},
+ { 18, 19, 20},
+ { 21, 22, 23},
+ { 24, 25, 26},
+ { 27, 28, 29},
+ { 30, 31, 32},
+ { 33, 34, 35},
+ { 36, 37, 38},
+ { 39, 40, 41},
+ { 42, 43, 44},
+ { 45, 46, 47},
+ { 48, 49, 50},
+ { 51, 52, 53},
+ { 54, 55, 56},
+ { 57, 58, 59},
+ { 60, 61, 62},
+ { 63, 64, 65},
+ { 66, 67, 68},
+ { 69, 70, 71},
+ { 72, 73, 74},
+ { 75, 76, 77},
+ { 78, 79, 80},
+ { 81, 82, 83},
+ { 84, 85, 86},
+ { 87, 88, 89},
+ { 90, 91, 92},
+ { 93, 94, 95},
+ { 96, 97, 98},
+ { 99, 100, 101},
+ {102, 103, 104},
+ {105, 106, 107},
+ {108, 109, 110},
+ {111, 112, 113},
+ {114, 115, 116},
+ {117, 118, 119},
+ {120, 121, 122},
+ {123, 124, 125},
+ {126, 127, 128},
+ {129, 130, 131},
+ {132, 133, 134},
+ {135, 136, 137},
+ {138, 139, 140},
+ {141, 142, 143},
+ {144, 145, 146},
+ {147, 148, 149},
+ {150, 151, 152},
+ {153, 154, 155},
+ {156, 157, 158},
+ {159, 160, 161},
+ {162, 163, 164},
+ {165, 166, 167},
+ {168, 169, 170},
+ {171, 172, 173},
+ {174, 175, 176},
+ {177, 178, 179},
+ {180, 181, 182},
+ {183, 184, 185},
+ {186, 187, 188},
+ {189, 190, 191}};
+
+ xt::xtensor<size_t, 2> dofsPeriodic_ =
+ {{ 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2},
+ { 9, 10, 11},
+ {12, 13, 14},
+ {15, 16, 17},
+ { 9, 10, 11},
+ {18, 19, 20},
+ {21, 22, 23},
+ {24, 25, 26},
+ {18, 19, 20},
+ { 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2},
+ {27, 28, 29},
+ {30, 31, 32},
+ {33, 34, 35},
+ {27, 28, 29},
+ {36, 37, 38},
+ {39, 40, 41},
+ {42, 43, 44},
+ {36, 37, 38},
+ {45, 46, 47},
+ {48, 49, 50},
+ {51, 52, 53},
+ {45, 46, 47},
+ {27, 28, 29},
+ {30, 31, 32},
+ {33, 34, 35},
+ {27, 28, 29},
+ {54, 55, 56},
+ {57, 58, 59},
+ {60, 61, 62},
+ {54, 55, 56},
+ {63, 64, 65},
+ {66, 67, 68},
+ {69, 70, 71},
+ {63, 64, 65},
+ {72, 73, 74},
+ {75, 76, 77},
+ {78, 79, 80},
+ {72, 73, 74},
+ {54, 55, 56},
+ {57, 58, 59},
+ {60, 61, 62},
+ {54, 55, 56},
+ { 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2},
+ { 9, 10, 11},
+ {12, 13, 14},
+ {15, 16, 17},
+ { 9, 10, 11},
+ {18, 19, 20},
+ {21, 22, 23},
+ {24, 25, 26},
+ {18, 19, 20},
+ { 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2}};
+
+ xt::xtensor<size_t, 2> nodesPeriodic_ =
+ {{ 0, 3},
+ { 0, 51},
+ { 0, 48},
+ { 0, 12},
+ { 0, 15},
+ { 0, 63},
+ { 0, 60},
+ { 1, 49},
+ { 2, 50},
+ { 1, 61},
+ { 2, 62},
+ { 1, 13},
+ { 2, 14},
+ {16, 19},
+ {32, 35},
+ {16, 31},
+ {32, 47},
+ {16, 28},
+ {32, 44},
+ { 4, 7},
+ { 8, 11},
+ { 4, 55},
+ { 8, 59},
+ { 4, 52},
+ { 8, 56},
+ { 5, 53},
+ { 6, 54},
+ { 9, 57},
+ {10, 58},
+ {20, 23},
+ {24, 27},
+ {36, 39},
+ {40, 43},
+ {17, 29},
+ {18, 30},
+ {33, 45},
+ {34, 46}};
+
+ size_t nodesOrigin_ = 0;
+
+ GooseFEM::Mesh::Hex8::Regular mesh(3, 3, 3);
+
+ REQUIRE(xt::allclose(mesh.coor(), coor_));
+ REQUIRE(xt::all(xt::equal(mesh.conn(), conn_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottom(), nodesBottom_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTop(), nodesTop_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesLeft(), nodesLeft_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesRight(), nodesRight_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFront(), nodesFront_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBack(), nodesBack_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontTopEdge(), nodesFrontTopEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontLeftEdge(), nodesFrontLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontRightEdge(), nodesFrontRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackBottomEdge(), nodesBackBottomEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackTopEdge(), nodesBackTopEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackLeftEdge(), nodesBackLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackRightEdge(), nodesBackRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomLeftEdge(), nodesBottomLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomRightEdge(), nodesBottomRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopLeftEdge(), nodesTopLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopRightEdge(), nodesTopRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontFace(), nodesFrontFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackFace(), nodesBackFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesLeftFace(), nodesLeftFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesRightFace(), nodesRightFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomFace(), nodesBottomFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopFace(), nodesTopFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontBottomOpenEdge(), nodesFrontBottomOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontTopOpenEdge(), nodesFrontTopOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontLeftOpenEdge(), nodesFrontLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontRightOpenEdge(), nodesFrontRightOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackBottomOpenEdge(), nodesBackBottomOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackTopOpenEdge(), nodesBackTopOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackLeftOpenEdge(), nodesBackLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackRightOpenEdge(), nodesBackRightOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomLeftOpenEdge(), nodesBottomLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomRightOpenEdge(), nodesBottomRightOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopLeftOpenEdge(), nodesTopLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopRightOpenEdge(), nodesTopRightOpenEdge_)));
+ REQUIRE(mesh.nodesFrontBottomLeftCorner() == nodesFrontBottomLeftCorner_);
+ REQUIRE(mesh.nodesFrontBottomRightCorner() == nodesFrontBottomRightCorner_);
+ REQUIRE(mesh.nodesFrontTopLeftCorner() == nodesFrontTopLeftCorner_);
+ REQUIRE(mesh.nodesFrontTopRightCorner() == nodesFrontTopRightCorner_);
+ REQUIRE(mesh.nodesBackBottomLeftCorner() == nodesBackBottomLeftCorner_);
+ REQUIRE(mesh.nodesBackBottomRightCorner() == nodesBackBottomRightCorner_);
+ REQUIRE(mesh.nodesBackTopLeftCorner() == nodesBackTopLeftCorner_);
+ REQUIRE(mesh.nodesBackTopRightCorner() == nodesBackTopRightCorner_);
+ REQUIRE(xt::all(xt::equal(mesh.dofs(), dofs_)));
+ REQUIRE(xt::all(xt::equal(mesh.dofsPeriodic(), dofsPeriodic_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesPeriodic(), nodesPeriodic_)));
+ REQUIRE(mesh.nodesOrigin() == nodesOrigin_);
+ }
+
+ SECTION("FineLayer - basic")
+ {
+ GooseFEM::Mesh::Hex8::FineLayer mesh(3, 3, 3);
+
+ xt::xtensor<double, 2> coor_ =
+ {{ 0., 0., 0.},
+ { 1., 0., 0.},
+ { 2., 0., 0.},
+ { 3., 0., 0.},
+ { 0., 0., 1.},
+ { 1., 0., 1.},
+ { 2., 0., 1.},
+ { 3., 0., 1.},
+ { 0., 0., 2.},
+ { 1., 0., 2.},
+ { 2., 0., 2.},
+ { 3., 0., 2.},
+ { 0., 0., 3.},
+ { 1., 0., 3.},
+ { 2., 0., 3.},
+ { 3., 0., 3.},
+ { 0., 1., 0.},
+ { 1., 1., 0.},
+ { 2., 1., 0.},
+ { 3., 1., 0.},
+ { 0., 1., 1.},
+ { 1., 1., 1.},
+ { 2., 1., 1.},
+ { 3., 1., 1.},
+ { 0., 1., 2.},
+ { 1., 1., 2.},
+ { 2., 1., 2.},
+ { 3., 1., 2.},
+ { 0., 1., 3.},
+ { 1., 1., 3.},
+ { 2., 1., 3.},
+ { 3., 1., 3.},
+ { 0., 2., 0.},
+ { 1., 2., 0.},
+ { 2., 2., 0.},
+ { 3., 2., 0.},
+ { 0., 2., 1.},
+ { 1., 2., 1.},
+ { 2., 2., 1.},
+ { 3., 2., 1.},
+ { 0., 2., 2.},
+ { 1., 2., 2.},
+ { 2., 2., 2.},
+ { 3., 2., 2.},
+ { 0., 2., 3.},
+ { 1., 2., 3.},
+ { 2., 2., 3.},
+ { 3., 2., 3.},
+ { 0., 3., 0.},
+ { 1., 3., 0.},
+ { 2., 3., 0.},
+ { 3., 3., 0.},
+ { 0., 3., 1.},
+ { 1., 3., 1.},
+ { 2., 3., 1.},
+ { 3., 3., 1.},
+ { 0., 3., 2.},
+ { 1., 3., 2.},
+ { 2., 3., 2.},
+ { 3., 3., 2.},
+ { 0., 3., 3.},
+ { 1., 3., 3.},
+ { 2., 3., 3.},
+ { 3., 3., 3.}};
+
+ xt::xtensor<double, 2> conn_ =
+ {{ 0, 1, 17, 16, 4, 5, 21, 20},
+ { 1, 2, 18, 17, 5, 6, 22, 21},
+ { 2, 3, 19, 18, 6, 7, 23, 22},
+ { 4, 5, 21, 20, 8, 9, 25, 24},
+ { 5, 6, 22, 21, 9, 10, 26, 25},
+ { 6, 7, 23, 22, 10, 11, 27, 26},
+ { 8, 9, 25, 24, 12, 13, 29, 28},
+ { 9, 10, 26, 25, 13, 14, 30, 29},
+ {10, 11, 27, 26, 14, 15, 31, 30},
+ {16, 17, 33, 32, 20, 21, 37, 36},
+ {17, 18, 34, 33, 21, 22, 38, 37},
+ {18, 19, 35, 34, 22, 23, 39, 38},
+ {20, 21, 37, 36, 24, 25, 41, 40},
+ {21, 22, 38, 37, 25, 26, 42, 41},
+ {22, 23, 39, 38, 26, 27, 43, 42},
+ {24, 25, 41, 40, 28, 29, 45, 44},
+ {25, 26, 42, 41, 29, 30, 46, 45},
+ {26, 27, 43, 42, 30, 31, 47, 46},
+ {32, 33, 49, 48, 36, 37, 53, 52},
+ {33, 34, 50, 49, 37, 38, 54, 53},
+ {34, 35, 51, 50, 38, 39, 55, 54},
+ {36, 37, 53, 52, 40, 41, 57, 56},
+ {37, 38, 54, 53, 41, 42, 58, 57},
+ {38, 39, 55, 54, 42, 43, 59, 58},
+ {40, 41, 57, 56, 44, 45, 61, 60},
+ {41, 42, 58, 57, 45, 46, 62, 61},
+ {42, 43, 59, 58, 46, 47, 63, 62}};
+
+ xt::xtensor<size_t, 1> nodesBottom_ =
+ { 0, 4, 8, 12, 1, 5, 9, 13, 2, 6, 10, 14, 3, 7, 11, 15};
+
+ xt::xtensor<size_t, 1> nodesTop_ =
+ {48, 52, 56, 60, 49, 53, 57, 61, 50, 54, 58, 62, 51, 55, 59, 63};
+
+ xt::xtensor<size_t, 1> nodesLeft_ =
+ { 0, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60};
+
+ xt::xtensor<size_t, 1> nodesRight_ =
+ { 3, 7, 11, 15, 19, 23, 27, 31, 35, 39, 43, 47, 51, 55, 59, 63};
+
+ xt::xtensor<size_t, 1> nodesFront_ =
+ { 0, 1, 2, 3, 16, 17, 18, 19, 32, 33, 34, 35, 48, 49, 50, 51};
+
+ xt::xtensor<size_t, 1> nodesBack_ =
+ {12, 13, 14, 15, 28, 29, 30, 31, 44, 45, 46, 47, 60, 61, 62, 63};
+
+ xt::xtensor<size_t, 1> nodesFrontTopEdge_ = {48, 49, 50, 51};
+ xt::xtensor<size_t, 1> nodesFrontLeftEdge_ = { 0, 16, 32, 48};
+ xt::xtensor<size_t, 1> nodesFrontRightEdge_ = { 3, 19, 35, 51};
+ xt::xtensor<size_t, 1> nodesBackBottomEdge_ = {12, 13, 14, 15};
+ xt::xtensor<size_t, 1> nodesBackTopEdge_ = {60, 61, 62, 63};
+ xt::xtensor<size_t, 1> nodesBackLeftEdge_ = {12, 28, 44, 60};
+ xt::xtensor<size_t, 1> nodesBackRightEdge_ = {15, 31, 47, 63};
+ xt::xtensor<size_t, 1> nodesBottomLeftEdge_ = { 0, 4, 8, 12};
+ xt::xtensor<size_t, 1> nodesBottomRightEdge_ = { 3, 7, 11, 15};
+ xt::xtensor<size_t, 1> nodesTopLeftEdge_ = {48, 52, 56, 60};
+ xt::xtensor<size_t, 1> nodesTopRightEdge_ = {51, 55, 59, 63};
+ xt::xtensor<size_t, 1> nodesFrontFace_ = {17, 18, 33, 34};
+ xt::xtensor<size_t, 1> nodesBackFace_ = {29, 30, 45, 46};
+ xt::xtensor<size_t, 1> nodesLeftFace_ = {20, 24, 36, 40};
+ xt::xtensor<size_t, 1> nodesRightFace_ = {23, 27, 39, 43};
+ xt::xtensor<size_t, 1> nodesBottomFace_ = { 5, 9, 6, 10};
+ xt::xtensor<size_t, 1> nodesTopFace_ = {53, 57, 54, 58};
+ xt::xtensor<size_t, 1> nodesFrontBottomOpenEdge_ = {1, 2};
+ xt::xtensor<size_t, 1> nodesFrontTopOpenEdge_ = {49, 50};
+ xt::xtensor<size_t, 1> nodesFrontLeftOpenEdge_ = {16, 32};
+ xt::xtensor<size_t, 1> nodesFrontRightOpenEdge_ = {19, 35};
+ xt::xtensor<size_t, 1> nodesBackBottomOpenEdge_ = {13, 14};
+ xt::xtensor<size_t, 1> nodesBackTopOpenEdge_ = {61, 62};
+ xt::xtensor<size_t, 1> nodesBackLeftOpenEdge_ = {28, 44};
+ xt::xtensor<size_t, 1> nodesBackRightOpenEdge_ = {31, 47};
+ xt::xtensor<size_t, 1> nodesBottomLeftOpenEdge_ = {4, 8};
+ xt::xtensor<size_t, 1> nodesBottomRightOpenEdge_ = { 7, 11};
+ xt::xtensor<size_t, 1> nodesTopLeftOpenEdge_ = {52, 56};
+ xt::xtensor<size_t, 1> nodesTopRightOpenEdge_ = {55, 59};
+ size_t nodesFrontBottomLeftCorner_ = 0;
+ size_t nodesFrontBottomRightCorner_ = 3;
+ size_t nodesFrontTopLeftCorner_ = 48;
+ size_t nodesFrontTopRightCorner_ = 51;
+ size_t nodesBackBottomLeftCorner_ = 12;
+ size_t nodesBackBottomRightCorner_ = 15;
+ size_t nodesBackTopLeftCorner_ = 60;
+ size_t nodesBackTopRightCorner_ = 63;
+
+ xt::xtensor<size_t, 2> dofs_ =
+ {{ 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 9, 10, 11},
+ { 12, 13, 14},
+ { 15, 16, 17},
+ { 18, 19, 20},
+ { 21, 22, 23},
+ { 24, 25, 26},
+ { 27, 28, 29},
+ { 30, 31, 32},
+ { 33, 34, 35},
+ { 36, 37, 38},
+ { 39, 40, 41},
+ { 42, 43, 44},
+ { 45, 46, 47},
+ { 48, 49, 50},
+ { 51, 52, 53},
+ { 54, 55, 56},
+ { 57, 58, 59},
+ { 60, 61, 62},
+ { 63, 64, 65},
+ { 66, 67, 68},
+ { 69, 70, 71},
+ { 72, 73, 74},
+ { 75, 76, 77},
+ { 78, 79, 80},
+ { 81, 82, 83},
+ { 84, 85, 86},
+ { 87, 88, 89},
+ { 90, 91, 92},
+ { 93, 94, 95},
+ { 96, 97, 98},
+ { 99, 100, 101},
+ {102, 103, 104},
+ {105, 106, 107},
+ {108, 109, 110},
+ {111, 112, 113},
+ {114, 115, 116},
+ {117, 118, 119},
+ {120, 121, 122},
+ {123, 124, 125},
+ {126, 127, 128},
+ {129, 130, 131},
+ {132, 133, 134},
+ {135, 136, 137},
+ {138, 139, 140},
+ {141, 142, 143},
+ {144, 145, 146},
+ {147, 148, 149},
+ {150, 151, 152},
+ {153, 154, 155},
+ {156, 157, 158},
+ {159, 160, 161},
+ {162, 163, 164},
+ {165, 166, 167},
+ {168, 169, 170},
+ {171, 172, 173},
+ {174, 175, 176},
+ {177, 178, 179},
+ {180, 181, 182},
+ {183, 184, 185},
+ {186, 187, 188},
+ {189, 190, 191}};
+
+ xt::xtensor<size_t, 2> dofsPeriodic_ =
+ {{ 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2},
+ { 9, 10, 11},
+ {12, 13, 14},
+ {15, 16, 17},
+ { 9, 10, 11},
+ {18, 19, 20},
+ {21, 22, 23},
+ {24, 25, 26},
+ {18, 19, 20},
+ { 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2},
+ {27, 28, 29},
+ {30, 31, 32},
+ {33, 34, 35},
+ {27, 28, 29},
+ {36, 37, 38},
+ {39, 40, 41},
+ {42, 43, 44},
+ {36, 37, 38},
+ {45, 46, 47},
+ {48, 49, 50},
+ {51, 52, 53},
+ {45, 46, 47},
+ {27, 28, 29},
+ {30, 31, 32},
+ {33, 34, 35},
+ {27, 28, 29},
+ {54, 55, 56},
+ {57, 58, 59},
+ {60, 61, 62},
+ {54, 55, 56},
+ {63, 64, 65},
+ {66, 67, 68},
+ {69, 70, 71},
+ {63, 64, 65},
+ {72, 73, 74},
+ {75, 76, 77},
+ {78, 79, 80},
+ {72, 73, 74},
+ {54, 55, 56},
+ {57, 58, 59},
+ {60, 61, 62},
+ {54, 55, 56},
+ { 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2},
+ { 9, 10, 11},
+ {12, 13, 14},
+ {15, 16, 17},
+ { 9, 10, 11},
+ {18, 19, 20},
+ {21, 22, 23},
+ {24, 25, 26},
+ {18, 19, 20},
+ { 0, 1, 2},
+ { 3, 4, 5},
+ { 6, 7, 8},
+ { 0, 1, 2}};
+
+ xt::xtensor<size_t, 2> nodesPeriodic_ =
+ {{ 0, 3},
+ { 0, 15},
+ { 0, 12},
+ { 0, 48},
+ { 0, 51},
+ { 0, 63},
+ { 0, 60},
+ { 1, 13},
+ { 2, 14},
+ { 1, 61},
+ { 2, 62},
+ { 1, 49},
+ { 2, 50},
+ { 4, 7},
+ { 8, 11},
+ { 4, 55},
+ { 8, 59},
+ { 4, 52},
+ { 8, 56},
+ {16, 19},
+ {32, 35},
+ {16, 31},
+ {32, 47},
+ {16, 28},
+ {32, 44},
+ {17, 29},
+ {18, 30},
+ {33, 45},
+ {34, 46},
+ {20, 23},
+ {24, 27},
+ {36, 39},
+ {40, 43},
+ { 5, 53},
+ { 9, 57},
+ { 6, 54},
+ {10, 58}};
+
+ size_t nodesOrigin_ = 0;
+
+ REQUIRE(xt::allclose(mesh.coor(), coor_));
+ REQUIRE(xt::all(xt::equal(mesh.conn(), conn_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottom(), nodesBottom_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTop(), nodesTop_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesLeft(), nodesLeft_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesRight(), nodesRight_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFront(), nodesFront_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBack(), nodesBack_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontTopEdge(), nodesFrontTopEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontLeftEdge(), nodesFrontLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontRightEdge(), nodesFrontRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackBottomEdge(), nodesBackBottomEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackTopEdge(), nodesBackTopEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackLeftEdge(), nodesBackLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackRightEdge(), nodesBackRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomLeftEdge(), nodesBottomLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomRightEdge(), nodesBottomRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopLeftEdge(), nodesTopLeftEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopRightEdge(), nodesTopRightEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontFace(), nodesFrontFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackFace(), nodesBackFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesLeftFace(), nodesLeftFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesRightFace(), nodesRightFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomFace(), nodesBottomFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopFace(), nodesTopFace_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontBottomOpenEdge(), nodesFrontBottomOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontTopOpenEdge(), nodesFrontTopOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontLeftOpenEdge(), nodesFrontLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesFrontRightOpenEdge(), nodesFrontRightOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackBottomOpenEdge(), nodesBackBottomOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackTopOpenEdge(), nodesBackTopOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackLeftOpenEdge(), nodesBackLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBackRightOpenEdge(), nodesBackRightOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomLeftOpenEdge(), nodesBottomLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesBottomRightOpenEdge(), nodesBottomRightOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopLeftOpenEdge(), nodesTopLeftOpenEdge_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesTopRightOpenEdge(), nodesTopRightOpenEdge_)));
+ REQUIRE(mesh.nodesFrontBottomLeftCorner() == nodesFrontBottomLeftCorner_);
+ REQUIRE(mesh.nodesFrontBottomRightCorner() == nodesFrontBottomRightCorner_);
+ REQUIRE(mesh.nodesFrontTopLeftCorner() == nodesFrontTopLeftCorner_);
+ REQUIRE(mesh.nodesFrontTopRightCorner() == nodesFrontTopRightCorner_);
+ REQUIRE(mesh.nodesBackBottomLeftCorner() == nodesBackBottomLeftCorner_);
+ REQUIRE(mesh.nodesBackBottomRightCorner() == nodesBackBottomRightCorner_);
+ REQUIRE(mesh.nodesBackTopLeftCorner() == nodesBackTopLeftCorner_);
+ REQUIRE(mesh.nodesBackTopRightCorner() == nodesBackTopRightCorner_);
+ REQUIRE(xt::all(xt::equal(mesh.dofs(), dofs_)));
+ REQUIRE(xt::all(xt::equal(mesh.dofsPeriodic(), dofsPeriodic_)));
+ REQUIRE(xt::all(xt::equal(mesh.nodesPeriodic(), nodesPeriodic_)));
+ REQUIRE(mesh.nodesOrigin() == nodesOrigin_);
+ }
+}