From: Peter Munch <peterrmuench@gmail.com>
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<dim, spacedim>::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<dim, spacedim>::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(