template <int dim>
void
-test_simple(MGDoFHandler<dim>& mgdofs, bool faces)
+test_simple(DoFHandler<dim>& dofs, bool faces)
{
SparsityPattern pattern;
SparseMatrix<double> matrix;
- const DoFHandler<dim>& dofs = mgdofs;
const FiniteElement<dim>& fe = dofs.get_fe();
pattern.reinit (dofs.n_dofs(), dofs.n_dofs(),
(GeometryInfo<dim>::faces_per_cell
cell != tr.end(); ++cell, ++cn)
cell->set_user_index(cn);
- MGDoFHandler<dim> dofs(tr);
+ DoFHandler<dim> dofs(tr);
dofs.distribute_dofs(fe);
dofs.initialize_local_block_info();
- deallog << "DoFHandler " << dofs.n_dofs() << " levels";
- for (unsigned int l=0;l<tr.n_levels();++l)
- deallog << ' ' << l << ':' << dofs.n_dofs(l);
- deallog << std::endl;
+ deallog << "DoFHandler " << dofs.n_dofs() << std::endl;
test_simple(dofs, false);
+ deallog << "now with jump terms" << std::endl;
test_simple(dofs, true);
}
DEAL::Triangulation levels 0:1 1:4 2:16
-DEAL::DoFHandler 16 levels 0:1 1:4 2:16
+DEAL::DoFHandler 16
(0,0) 10.0000
(0,1) 0.00000
(0,2) 0.00000
(15,13) 0.00000
(15,14) 0.00000
(15,15) 10.0000
+DEAL::now with jump terms
(0,0) 14.0000
(0,1) -1.00000
(0,2) -1.00000
(15,14) -1.00000
(15,15) 14.0000
DEAL::Triangulation levels 0:1 1:4 2:16
-DEAL::DoFHandler 25 levels 0:4 1:9 2:25
+DEAL::DoFHandler 25
(0,0) 10.0000
(0,1) 10.0000
(0,2) 10.0000
(24,22) 10.0000
(24,23) 10.0000
(24,24) 10.0000
+DEAL::now with jump terms
(0,0) 14.0000
(0,1) 13.0000
(0,2) 13.0000
(24,23) 13.0000
(24,24) 14.0000
DEAL::Triangulation levels 0:1 1:4 2:16
-DEAL::DoFHandler 41 levels 0:5 1:13 2:41
+DEAL::DoFHandler 41
(0,0) 10.0000
(0,1) 10.0000
(0,2) 10.0000
(40,38) 0.00000
(40,39) 0.00000
(40,40) 10.0000
+DEAL::now with jump terms
(0,0) 14.0000
(0,1) 13.0000
(0,2) 13.0000