From: David Wells Date: Tue, 7 Feb 2023 02:04:05 +0000 (-0500) Subject: QProjector: minor orientation fixes. X-Git-Tag: v9.5.0-rc1~577^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=1e69380a19152972704fa8d3d07505e8c7f2fdab;p=dealii.git QProjector: minor orientation fixes. --- 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];