]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
hanging node constraints work now with new FETools::compute_face_embedding_matrices
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 20 Oct 2010 22:20:46 +0000 (22:20 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 20 Oct 2010 22:20:46 +0000 (22:20 +0000)
git-svn-id: https://svn.dealii.org/trunk@22403 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_bdm.cc

index 1822c2788a8e5638ce71ab1ed857b3088213b84b..276ee00f8e161c127073bd236481f4b091301fc5 100644 (file)
@@ -74,20 +74,20 @@ FE_BDM<dim>::FE_BDM (const unsigned int deg)
   this->reinit_restriction_and_prolongation_matrices(true, true);
   FETools::compute_embedding_matrices (*this, this->prolongation, true);
   
-//   FullMatrix<double> face_embeddings[GeometryInfo<dim>::subfaces_per_face];
-//   for (unsigned int i=0; i<GeometryInfo<dim>::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<GeometryInfo<dim>::subfaces_per_face;++d)
-//     for (unsigned int i=0;i<face_embeddings[d].m();++i)
-//       {
-//     for (unsigned int j=0;j<face_embeddings[d].n();++j)
-//       this->interface_constraints(target_row,j) = face_embeddings[d](i,j);
-//     ++target_row;
-//       }
+  FullMatrix<double> face_embeddings[GeometryInfo<dim>::max_children_per_face];
+  for (unsigned int i=0; i<GeometryInfo<dim>::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<GeometryInfo<dim>::max_children_per_face;++d)
+    for (unsigned int i=0;i<face_embeddings[d].m();++i)
+      {
+       for (unsigned int j=0;j<face_embeddings[d].n();++j)
+         this->interface_constraints(target_row,j) = face_embeddings[d](i,j);
+       ++target_row;
+      }
 }
 
 

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.