]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Don't build hanging node constraints if base element doesn't do so.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 17 Feb 2003 00:05:14 +0000 (00:05 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 17 Feb 2003 00:05:14 +0000 (00:05 +0000)
git-svn-id: https://svn.dealii.org/trunk@7104 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_system.cc

index e146ac1fbc26c4a97d5f7c6b8d2897b929b26772..53aca37dc641b29cbe3ce969fa865e81e040dd9b 100644 (file)
@@ -1224,6 +1224,18 @@ FESystem<dim>::build_face_tables()
 template <int dim>
 void FESystem<dim>::build_interface_constraints () 
 {
+                                   // check whether all base elements
+                                   // implement their interface
+                                   // constraint matrices. if this is
+                                   // not the case, then leave the
+                                   // interface costraints of this
+                                   // composed element empty as well;
+                                   // however, the rest of the element
+                                   // is usable
+  for (unsigned int base=0; base<n_base_elements(); ++base)
+    if (base_element(base).constraints_are_implemented() == false)
+      return;
+  
   this->interface_constraints.
     TableBase<2,double>::reinit (this->interface_constraints_size());
   

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.