Fix FESystem::initialize for pyramids/wedges

diff --git a/source/fe/fe_system.cc b/source/fe/fe_system.cc
index 212a8636e8ad8d4727ca1c8c867ca6408f208e28..2975be1af8460e00ffd75b95ca6ab44518ae8b26 100644 (file)
--- 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(