Replace ucd_to_deal with a local vertex numbering in Triangulation

diff --git a/source/grid/tria.cc b/source/grid/tria.cc
index 78ca6cecbb3edb0e8ed149a3bce7ca15c5f9749b..8d234f529f8bff4a14a331b0cbad538148b48c41 100644 (file)
--- a/source/grid/tria.cc
+++ b/source/grid/tria.cc
@@ -447,13 +447,16 @@ namespace
                         SubCellData &             subcelldata)
   {
     unsigned int tmp[GeometryInfo<3>::vertices_per_cell];
+    static constexpr std::array<unsigned int,
+                                GeometryInfo<3>::vertices_per_cell>
+      local_vertex_numbering{{0, 1, 5, 4, 2, 3, 7, 6}};
     for (auto &cell : cells)
       if (cell.vertices.size() == GeometryInfo<3>::vertices_per_cell)
         {
           for (const unsigned int i : GeometryInfo<3>::vertex_indices())
             tmp[i] = cell.vertices[i];
           for (const unsigned int i : GeometryInfo<3>::vertex_indices())
-            cell.vertices[GeometryInfo<3>::ucd_to_deal[i]] = tmp[i];
+            cell.vertices[local_vertex_numbering[i]] = tmp[i];
         }
 
     // now points in boundary quads