From: David Wells <drwells@email.unc.edu>
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];