From: wolf Date: Mon, 17 Feb 2003 00:05:14 +0000 (+0000) Subject: Don't build hanging node constraints if base element doesn't do so. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7eaad12e88748a9375295dc494adff320b2d3508;p=dealii-svn.git Don't build hanging node constraints if base element doesn't do so. git-svn-id: https://svn.dealii.org/trunk@7104 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_system.cc b/deal.II/deal.II/source/fe/fe_system.cc index e146ac1fbc..53aca37dc6 100644 --- a/deal.II/deal.II/source/fe/fe_system.cc +++ b/deal.II/deal.II/source/fe/fe_system.cc @@ -1224,6 +1224,18 @@ FESystem::build_face_tables() template void FESystem::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; baseinterface_constraints. TableBase<2,double>::reinit (this->interface_constraints_size());