From 9fee04c195a2e8cf25ec4f9903ea771f22717360 Mon Sep 17 00:00:00 2001 From: David Wells Date: Sat, 13 Apr 2024 18:11:44 -0400 Subject: [PATCH] Explicitly use the combined orientation. --- include/deal.II/grid/tria_accessor.templates.h | 4 +++- source/grid/tria.cc | 18 +++++++++++++++--- 2 files changed, 18 insertions(+), 4 deletions(-) diff --git a/include/deal.II/grid/tria_accessor.templates.h b/include/deal.II/grid/tria_accessor.templates.h index 6fd99cee44..79932777cb 100644 --- a/include/deal.II/grid/tria_accessor.templates.h +++ b/include/deal.II/grid/tria_accessor.templates.h @@ -973,7 +973,9 @@ namespace internal corner); const auto vertex_within_line_index = accessor.reference_cell().standard_to_real_face_vertex( - vertex_index, line_index, accessor.line_orientation(line_index)); + vertex_index, + line_index, + accessor.combined_face_orientation(line_index)); return accessor.line(line_index) ->vertex_index(vertex_within_line_index); diff --git a/source/grid/tria.cc b/source/grid/tria.cc index de3681962d..02b724666e 100644 --- a/source/grid/tria.cc +++ b/source/grid/tria.cc @@ -6425,7 +6425,11 @@ namespace internal s.insert(i); #endif - new_quad->set_line_orientation(f, orientation); + new_quad->set_line_orientation( + f, + orientation == + ReferenceCell:: + default_combined_face_orientation()); } #ifdef DEBUG AssertDimension(s.size(), 3); @@ -6972,7 +6976,11 @@ namespace internal make_array_view(vertices_0), make_array_view(vertices_1)); - new_quad->set_line_orientation(l, orientation); + new_quad->set_line_orientation( + l, + orientation == + ReferenceCell:: + default_combined_face_orientation()); // on a hex, inject the status of the current line // also to the line on the other quad along the @@ -6981,7 +6989,11 @@ namespace internal ReferenceCells::Hexahedron) new_quads[representative_lines[q % 4][1] + q - (q % 4)] - ->set_line_orientation(l, orientation); + ->set_line_orientation( + l, + orientation == + ReferenceCell:: + default_combined_face_orientation()); } } } -- 2.39.5