From 1e69380a19152972704fa8d3d07505e8c7f2fdab Mon Sep 17 00:00:00 2001 From: David Wells 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> { + [](const auto & face, + const unsigned char orientation) -> std::vector> { std::array, 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> { + [](const auto & face, + const unsigned char orientation) -> std::vector> { std::array, 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 &>(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