]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix FESystem::initialize for pyramids/wedges 13843/head
authorPeter Munch <peterrmuench@gmail.com>
Fri, 27 May 2022 21:13:49 +0000 (23:13 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Fri, 27 May 2022 21:13:49 +0000 (23:13 +0200)
source/fe/fe_system.cc

index 212a8636e8ad8d4727ca1c8c867ca6408f208e28..2975be1af8460e00ffd75b95ca6ab44518ae8b26 100644 (file)
@@ -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(

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.