+static int Hexa_oppF[] = { 1, 0, 4, 5, 2, 3 }; // oppopsite facet indices
+
+/*
+// N8 +------+ N9
+// / \
+// / \
+// N7 + + N10
+// \ /
+// \ /
+// N6 +------+ N11
+// HEXAGONAL PRISM
+// N2 +------+ N3
+// / \
+// / \
+// N1 + + N4
+// \ /
+// \ /
+// N0 +------+ N5
+*/
+static int HexPrism_F [8][7] = { // FORWARD
+ { 0, 1, 2, 3, 4, 5, 0 },
+ { 6,11,10, 9, 8, 7, 6 },
+ { 0, 6, 7, 1, 0, 0, 0 },
+ { 1, 7, 8, 2, 1, 1, 1 },
+ { 2, 8, 9, 3, 2, 2, 2 },
+ { 3, 9,10, 4, 3, 3, 3 },
+ { 4,10,11, 5, 4, 4, 4 },
+ { 5,11, 6, 0, 5, 5, 5 }};
+static int HexPrism_RE [8][7] = { // REVERSED -> EXTERNAL
+ { 0, 5, 4, 3, 2, 1, 0 },
+ { 6,11,10, 9, 8, 7, 6 },
+ { 0, 6, 7, 1, 0, 0, 0 },
+ { 1, 7, 8, 2, 1, 1, 1 },
+ { 2, 8, 9, 3, 2, 2, 2 },
+ { 3, 9,10, 4, 3, 3, 3 },
+ { 4,10,11, 5, 4, 4, 4 },
+ { 5,11, 6, 0, 5, 5, 5 }};
+static int HexPrism_nbN [] = { 6, 6, 4, 4, 4, 4, 4, 4 };