From d764eddff04092da4ae1d894d3c47db6403ed81b Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Fri, 27 May 2022 23:13:49 +0200 Subject: [PATCH] Fix FESystem::initialize for pyramids/wedges --- source/fe/fe_system.cc | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) 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( -- 2.39.5