Commit 78d5af37 authored by jonathan.froehlich's avatar jonathan.froehlich
Browse files

[195-fix-tetrahedron-p4-shape-derivative] Added higher order tests

parent 541c0358
......@@ -152,7 +152,21 @@ void LagrangeNodalPointsTetrahedron(const cell &c, vector<PointT<T, sDim, tDim>>
LagrangeNodalPointsTetrahedron(*c, z, degree);
}
// Used to identify NodalPointsOnFace
// Used to identify NodalPointsOnFace * Numbering of nodal points on tetrahedron (example P4)
// *
// * front: rear:
// *
// * 34
// * | -
// * 31 32 33
// * | - | -
// * 25 26 27 28 29 30
// * | - | - | -
// * 15 16 17 18 19 20 21 22 23 24
// * | / - | / - | / - | / -
// * 0 -- 1 -- 2 -- 3 -- 4 5 -- 6 -- 7 -- 8 9 -- 10 -- 11 12 -- 13 14
// * / \ / \
// */
static std::vector<std::vector<std::vector<int>>> tetrahdronNodalMap{
{
{0, 2, 1},
......@@ -173,10 +187,10 @@ static std::vector<std::vector<std::vector<int>>> tetrahdronNodalMap{
{0, 1, 2, 3, 10, 11, 12, 16, 17, 19}
},
{
{0, 5, 9, 12, 14, 1, 6, 10, 13, 2, 7, 11, 3, 8, 4},
{4, 8, 11, 13, 14, 18, 21, 23, 24, 27, 29, 30, 32, 33, 34},
{0, 15, 25, 31, 34, 5, 19, 28, 33, 9, 22, 30, 12, 24, 14},
{0, 1, 2, 3, 4, 15, 16, 17, 18, 25, 26, 27, 31, 32, 34}
{0, 5, 9, 12, 14, 1, 6, 10, 13, 2, 7, 11, 3, 8, 4}, // OK
{4, 8, 11, 13, 14, 18, 21, 23, 24, 27, 29, 30, 32, 33, 34}, // OK
{0, 15, 25, 31, 34, 5, 19, 28, 33, 9, 22, 30, 12, 24, 14}, // OK
{0, 1, 2, 3, 4, 15, 16, 17, 18, 25, 26, 27, 31, 32, 34} // OK
}
};
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment