From: Nils Margenberg <margenbn@hsu-hh.de>
Date: Sun, 15 Jan 2023 21:09:17 +0000 (+0100)
Subject: Replace ucd_to_deal with a local vertex numbering in Triangulation
X-Git-Tag: v9.5.0-rc1~642^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a41bc835186e4c9795daf023b65646657957eb7f;p=dealii.git

Replace ucd_to_deal with a local vertex numbering in Triangulation
---

diff --git a/source/grid/tria.cc b/source/grid/tria.cc
index 78ca6cecbb..8d234f529f 100644
--- 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