From 5845144098eac854c8f7d20f0ea469ddcedc461b Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Thu, 3 Jul 2003 11:49:00 +0000 Subject: [PATCH] test dg numbering git-svn-id: https://svn.dealii.org/trunk@7835 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/deal.II/dof_renumbering.cc | 86 ++++++++++++++++++++++---------- 1 file changed, 59 insertions(+), 27 deletions(-) diff --git a/tests/deal.II/dof_renumbering.cc b/tests/deal.II/dof_renumbering.cc index 72c55b5dcf..4968214a0a 100644 --- a/tests/deal.II/dof_renumbering.cc +++ b/tests/deal.II/dof_renumbering.cc @@ -29,6 +29,7 @@ #include #include #include +#include #include #include @@ -70,44 +71,40 @@ print_dofs (const MGDoFHandler &dof, unsigned int level) } - template void -check () +check_renumbering(MGDoFHandler& mgdof, bool discontinuous) { - Functions::CosineFunction cosine; - - Triangulation tr; - if (dim==2) - GridGenerator::hyper_ball(tr, Point(), 1); - else - GridGenerator::hyper_cube(tr, -1,1); - tr.refine_global (1); - tr.begin_active()->set_refine_flag (); - tr.execute_coarsening_and_refinement (); - if (dim==1) - tr.refine_global(2); - - FESystem element (FE_Q(2), 2, FE_DGQ(1), 1); - MGDoFHandler mgdof(tr); + const FiniteElement& element = mgdof.get_fe(); DoFHandler& dof = mgdof; - dof.distribute_dofs(element); - + // Prepare a reordering of // components for later use std::vector order(element.n_components()); for (unsigned int i=0; i direction; + for (unsigned int i=0;i +void +check () +{ + Triangulation tr; + if (dim==2) + GridGenerator::hyper_ball(tr, Point(), 1); + else + GridGenerator::hyper_cube(tr, -1,1); + tr.refine_global (1); + tr.begin_active()->set_refine_flag (); + tr.execute_coarsening_and_refinement (); + if (dim==1) + tr.refine_global(2); + + MGDoFHandler mgdof(tr); + + FESystem e1 (FE_Q(2), 2, FE_DGQ(1), 1); + mgdof.distribute_dofs(e1); + check_renumbering(mgdof, false); + mgdof.clear(); + + FESystem e2 (FE_DGP(2), 2, FE_DGQ(1), 1); + mgdof.distribute_dofs(e2); + check_renumbering(mgdof, true); + mgdof.clear(); } -- 2.39.5