From: Peter Munch <peterrmuench@gmail.com>
Date: Sat, 5 Dec 2020 11:37:42 +0000 (+0100)
Subject: Impelment ReferenceCell::internal::Info::Pyramid::standard_to_real_face_line
X-Git-Tag: v9.3.0-rc1~773^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ef20aaa6e8ed892063aeb2904bc21d4d78c6102f;p=dealii.git

Impelment ReferenceCell::internal::Info::Pyramid::standard_to_real_face_line
---

diff --git a/include/deal.II/grid/reference_cell.h b/include/deal.II/grid/reference_cell.h
index 5402c73267..87a4cc6c84 100644
--- a/include/deal.II/grid/reference_cell.h
+++ b/include/deal.II/grid/reference_cell.h
@@ -834,19 +834,26 @@ namespace ReferenceCell
           const unsigned int  face,
           const unsigned char face_orientation) const override
         {
-          Assert(false, ExcNotImplemented());
-
-          (void)face;
-
-          static const std::array<std::array<unsigned int, 3>, 6> table = {
-            {{{2, 1, 0}},
-             {{0, 1, 2}},
-             {{1, 2, 0}},
-             {{0, 2, 1}},
-             {{1, 0, 2}},
-             {{2, 0, 1}}}};
+          if (face == 0) // QUAD
+            {
+              return GeometryInfo<3>::standard_to_real_face_line(
+                line,
+                get_bit(face_orientation, 0),
+                get_bit(face_orientation, 2),
+                get_bit(face_orientation, 1));
+            }
+          else // TRI
+            {
+              static const std::array<std::array<unsigned int, 3>, 6> table = {
+                {{{2, 1, 0}},
+                 {{0, 1, 2}},
+                 {{1, 2, 0}},
+                 {{0, 2, 1}},
+                 {{1, 0, 2}},
+                 {{2, 0, 1}}}};
 
-          return table[face_orientation][line];
+              return table[face_orientation][line];
+            }
         }
 
         bool