From: kanschat Date: Wed, 20 Oct 2010 22:20:46 +0000 (+0000) Subject: hanging node constraints work now with new FETools::compute_face_embedding_matrices X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9f7d6f96c803e55515d7f8b09c0fb4df7e231689;p=dealii-svn.git hanging node constraints work now with new FETools::compute_face_embedding_matrices git-svn-id: https://svn.dealii.org/trunk@22403 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_bdm.cc b/deal.II/deal.II/source/fe/fe_bdm.cc index 1822c2788a..276ee00f8e 100644 --- a/deal.II/deal.II/source/fe/fe_bdm.cc +++ b/deal.II/deal.II/source/fe/fe_bdm.cc @@ -74,20 +74,20 @@ FE_BDM::FE_BDM (const unsigned int deg) this->reinit_restriction_and_prolongation_matrices(true, true); FETools::compute_embedding_matrices (*this, this->prolongation, true); -// FullMatrix face_embeddings[GeometryInfo::subfaces_per_face]; -// for (unsigned int i=0; i::subfaces_per_face; ++i) -// face_embeddings[i].reinit (this->dofs_per_face, this->dofs_per_face); -// FETools::compute_face_embedding_matrices(*this, face_embeddings, 0, 0); -// this->interface_constraints.reinit((1<<(dim-1)) * this->dofs_per_face, -// this->dofs_per_face); -// unsigned int target_row=0; -// for (unsigned int d=0;d::subfaces_per_face;++d) -// for (unsigned int i=0;iinterface_constraints(target_row,j) = face_embeddings[d](i,j); -// ++target_row; -// } + FullMatrix face_embeddings[GeometryInfo::max_children_per_face]; + for (unsigned int i=0; i::max_children_per_face; ++i) + face_embeddings[i].reinit (this->dofs_per_face, this->dofs_per_face); + FETools::compute_face_embedding_matrices(*this, face_embeddings, 0, 0); + this->interface_constraints.reinit((1<<(dim-1)) * this->dofs_per_face, + this->dofs_per_face); + unsigned int target_row=0; + for (unsigned int d=0;d::max_children_per_face;++d) + for (unsigned int i=0;iinterface_constraints(target_row,j) = face_embeddings[d](i,j); + ++target_row; + } }