From 1e69380a19152972704fa8d3d07505e8c7f2fdab Mon Sep 17 00:00:00 2001
From: David Wells <drwells@email.unc.edu>
Date: Mon, 6 Feb 2023 21:04:05 -0500
Subject: [PATCH] QProjector: minor orientation fixes.

---
 source/base/qprojector.cc | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/source/base/qprojector.cc b/source/base/qprojector.cc
index cafd48cf3d..60cf0a7442 100644
--- a/source/base/qprojector.cc
+++ b/source/base/qprojector.cc
@@ -644,7 +644,8 @@ QProjector<3>::project_to_all_faces(const ReferenceCell &     reference_cell,
                                     const hp::QCollection<2> &quadrature)
 {
   const auto support_points_tri =
-    [](const auto &face, const auto &orientation) -> std::vector<Point<3>> {
+    [](const auto &        face,
+       const unsigned char orientation) -> std::vector<Point<3>> {
     std::array<Point<3>, 3> vertices;
     std::copy_n(face.first.begin(), face.first.size(), vertices.begin());
     const auto temp =
@@ -655,7 +656,8 @@ QProjector<3>::project_to_all_faces(const ReferenceCell &     reference_cell,
   };
 
   const auto support_points_quad =
-    [](const auto &face, const auto &orientation) -> std::vector<Point<3>> {
+    [](const auto &        face,
+       const unsigned char orientation) -> std::vector<Point<3>> {
     std::array<Point<3>, 4> vertices;
     std::copy_n(face.first.begin(), face.first.size(), vertices.begin());
     const auto temp =
@@ -688,8 +690,8 @@ QProjector<3>::project_to_all_faces(const ReferenceCell &     reference_cell,
             static_cast<const ScalarPolynomialsBase<2> &>(poly_quad);
 
         // ... and over all possible orientations
-        for (unsigned int orientation = 0;
-             orientation < (n_shape_functions * 2);
+        for (unsigned char orientation = 0;
+             orientation < reference_cell.n_face_orientations(face_no);
              ++orientation)
           {
             const auto &face = faces[face_no];
-- 
2.39.5