In an assertion, it was checked whether vertex dof indices for

different children of a face matched each other, for consistency. This
failed if the FE in use does not have dofs on vertices at all. Check
for the vertex index instead of the vertex _dof_ index now.

git-svn-id: https://svn.dealii.org/trunk@5534 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/dofs/dof_tools.cc b/deal.II/deal.II/source/dofs/dof_tools.cc
index 9ce0174e7e97fe688aff54d84f205acf83578e34..85809786b0abc20593db879c9f60b8ec6cea0152 100644 (file)
--- a/deal.II/deal.II/source/dofs/dof_tools.cc
+++ b/deal.II/deal.II/source/dofs/dof_tools.cc
@@ -824,12 +824,15 @@ void DoFTools::make_hanging_node_constraints (const DoFHandler<3> &dof_handler,
                  ExcInternalError());
          
          next_index = 0;
-         Assert ((face->child(0)->vertex_dof_index(2,0) ==
-                  face->child(1)->vertex_dof_index(3,0)) &&
-                 (face->child(0)->vertex_dof_index(2,0) ==
-                  face->child(2)->vertex_dof_index(0,0)) &&
-                 (face->child(0)->vertex_dof_index(2,0) ==
-                  face->child(3)->vertex_dof_index(1,0)),
+
+                                          // assert some consistency
+                                          // assumptions
+         Assert ((face->child(0)->vertex_index(2) ==
+                  face->child(1)->vertex_index(3)) &&
+                 (face->child(0)->vertex_index(2) ==
+                  face->child(2)->vertex_index(0)) &&
+                 (face->child(0)->vertex_index(2) ==
+                  face->child(3)->vertex_index(1)),
                  ExcInternalError());
          for (unsigned int dof=0; dof!=fe.dofs_per_vertex; ++dof)
            dofs_on_children[next_index++]