From a41bc835186e4c9795daf023b65646657957eb7f Mon Sep 17 00:00:00 2001 From: Nils Margenberg Date: Sun, 15 Jan 2023 22:09:17 +0100 Subject: [PATCH] Replace ucd_to_deal with a local vertex numbering in Triangulation --- source/grid/tria.cc | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) 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::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 -- 2.39.5