(++cell)->set_active_fe_index(1);
dof_handler.distribute_dofs (fe);
+
+ deallog << "n_dofs=" << dof_handler.n_dofs() << std::endl;
+ for (auto cell : dof_handler.active_cell_iterators())
+ {
+ deallog << cell << ": ";
+ std::vector<types::global_dof_index> dof_indices (cell->get_fe().dofs_per_cell);
+ cell->get_dof_indices (dof_indices);
+ for (auto i : dof_indices)
+ deallog << i << ' ';
+ deallog << std::endl;
+ }
+
ConstraintMatrix hanging_node_constraints;
DoFTools::make_hanging_node_constraints (dof_handler,
hanging_node_constraints);
hanging_node_constraints.close ();
// print constraints. there shouldn't be any
- deallog << hanging_node_constraints.n_constraints() << std::endl;
+ deallog << "n_constraints: "
+ << hanging_node_constraints.n_constraints() << std::endl;
hanging_node_constraints.print (deallog.get_file_stream());
}
-DEAL::0
-DEAL::0
-DEAL::0
+DEAL::n_dofs=5
+DEAL::1.0: 0 1
+DEAL::1.1: 1 2 3 4
+DEAL::n_constraints: 0
+DEAL::n_dofs=13
+DEAL::1.0: 0 1 2 3
+DEAL::1.1: 1 4 5 6 3 7 11 8
+DEAL::1.2: 2 3 9 10
+DEAL::1.3: 3 11 10 12
+DEAL::n_constraints: 0
+DEAL::n_dofs=35
+DEAL::1.0: 0 1 2 3 4 5 6 7
+DEAL::1.1: 1 8 9 10 3 11 21 12 5 13 29 14 7 15 23 16
+DEAL::1.2: 2 3 17 18 6 7 19 20
+DEAL::1.3: 3 21 18 22 7 23 20 24
+DEAL::1.4: 4 5 6 7 25 26 27 28
+DEAL::1.5: 5 29 7 23 26 30 28 31
+DEAL::1.6: 6 7 19 20 27 28 32 33
+DEAL::1.7: 7 23 20 24 28 31 33 34
+DEAL::n_constraints: 0