From: David Wells Date: Mon, 25 Mar 2024 16:25:35 +0000 (-0400) Subject: ReferenceCell: permit face index calculations for lines. X-Git-Tag: v9.6.0-rc1~431^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=17ddc95f8406ea151b01f6ab8606004f5a89e27f;p=dealii.git ReferenceCell: permit face index calculations for lines. This way we can call standard_to_real_face_vertex for all supported topological dimensions. --- diff --git a/include/deal.II/grid/reference_cell.h b/include/deal.II/grid/reference_cell.h index 19b760fa1b..803a5a00fb 100644 --- a/include/deal.II/grid/reference_cell.h +++ b/include/deal.II/grid/reference_cell.h @@ -1877,7 +1877,7 @@ ReferenceCell::face_to_cell_vertices( if (get_dimension() == 1) Assert(combined_face_orientation == ReferenceCell::default_combined_face_orientation(), - ExcMessage("In 1D, all cells must have the default orientation.")); + ExcMessage("In 1D, all faces must have the default orientation.")); else AssertIndexRange(combined_face_orientation, n_face_orientations(face)); @@ -1988,9 +1988,13 @@ ReferenceCell::standard_to_real_face_vertex( switch (this->kind) { case ReferenceCells::Vertex: - case ReferenceCells::Line: DEAL_II_NOT_IMPLEMENTED(); break; + case ReferenceCells::Line: + Assert(face_orientation == default_combined_face_orientation(), + ExcMessage( + "In 1D, all faces must have the default orientation.")); + return vertex; case ReferenceCells::Triangle: case ReferenceCells::Quadrilateral: return line_vertex_permutations[face_orientation][vertex];