From: Peter Munch Date: Fri, 27 May 2022 21:13:49 +0000 (+0200) Subject: Fix FESystem::initialize for pyramids/wedges X-Git-Tag: v9.4.0-rc1~125^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F13843%2Fhead;p=dealii.git Fix FESystem::initialize for pyramids/wedges --- diff --git a/source/fe/fe_system.cc b/source/fe/fe_system.cc index 212a8636e8..2975be1af8 100644 --- a/source/fe/fe_system.cc +++ b/source/fe/fe_system.cc @@ -1536,7 +1536,9 @@ FESystem::build_interface_constraints() { // TODO: the implementation makes the assumption that all faces have the // same number of dofs - AssertDimension(this->n_unique_faces(), 1); + if (this->n_unique_faces() != 1) + return; + const unsigned int face_no = 0; // check whether all base elements implement their interface constraint @@ -1809,7 +1811,7 @@ FESystem::initialize( for (unsigned int face_no = 0; face_no < this->n_unique_faces(); ++face_no) { - this->face_system_to_component_table[0].resize( + this->face_system_to_component_table[face_no].resize( this->n_dofs_per_face(face_no)); FETools::Compositing::build_face_tables(