From c25eaaead5941a69a08bb6d43bd5879915ac36e1 Mon Sep 17 00:00:00 2001 From: Nils Margenberg Date: Sun, 15 Jan 2023 22:09:30 +0100 Subject: [PATCH] Replace ucd_to_deal with a local vertex numbering in grid_reordering.cc --- source/grid/grid_reordering.cc | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/source/grid/grid_reordering.cc b/source/grid/grid_reordering.cc index d031215988..7e37ecd9f1 100644 --- a/source/grid/grid_reordering.cc +++ b/source/grid/grid_reordering.cc @@ -40,6 +40,9 @@ namespace * do the reordering of their * arguments in-place. */ + constexpr std::array local_vertex_numbering{ + {0, 1, 5, 4, 2, 3, 7, 6}}; + void reorder_new_to_old_style(std::vector> &) {} @@ -62,7 +65,7 @@ namespace 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[i] = tmp[GeometryInfo<3>::ucd_to_deal[i]]; + cell.vertices[i] = tmp[local_vertex_numbering[i]]; } } @@ -93,7 +96,7 @@ namespace 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]; } } } // namespace -- 2.39.5