From: Peter Munch Date: Sun, 8 May 2022 14:32:12 +0000 (+0200) Subject: Implement ReferenceCell::face_to_cell_lines() for pyramids and wedges X-Git-Tag: v9.4.0-rc1~262^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F13692%2Fhead;p=dealii.git Implement ReferenceCell::face_to_cell_lines() for pyramids and wedges --- diff --git a/include/deal.II/grid/reference_cell.h b/include/deal.II/grid/reference_cell.h index 7e38aa9bc1..5431a59092 100644 --- a/include/deal.II/grid/reference_cell.h +++ b/include/deal.II/grid/reference_cell.h @@ -1404,6 +1404,8 @@ ReferenceCell::face_to_cell_lines(const unsigned int face, AssertIndexRange(face, n_faces()); AssertIndexRange(line, face_reference_cell(face).n_lines()); + static constexpr unsigned int X = numbers::invalid_unsigned_int; + if (*this == ReferenceCells::Vertex) { Assert(false, ExcNotImplemented()); @@ -1440,11 +1442,25 @@ ReferenceCell::face_to_cell_lines(const unsigned int face, } else if (*this == ReferenceCells::Pyramid) { - Assert(false, ExcNotImplemented()); + static const ndarray table = {{{{0, 1, 2, 3}}, + {{0, 6, 4, X}}, + {{1, 5, 7, X}}, + {{2, 4, 5, X}}, + {{3, 7, 6, 2}}}}; + + return table[face] + [standard_to_real_face_line(line, face, face_orientation)]; } else if (*this == ReferenceCells::Wedge) { - Assert(false, ExcNotImplemented()); + static const ndarray table = {{{{0, 2, 1, X}}, + {{3, 4, 5, X}}, + {{6, 7, 0, 3}}, + {{7, 8, 1, 4}}, + {{8, 6, 5, 2}}}}; + + return table[face] + [standard_to_real_face_line(line, face, face_orientation)]; } else if (*this == ReferenceCells::Hexahedron) {