From: Guido Kanschat Date: Fri, 1 Aug 2003 21:55:41 +0000 (+0000) Subject: first multigrid test X-Git-Tag: v8.0.0~16316 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2dd52e090c3ded25943cc85281fa81f9b6b14062;p=dealii.git first multigrid test git-svn-id: https://svn.dealii.org/trunk@7894 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/multigrid/transfer.cc b/tests/multigrid/transfer.cc index 4b9f4c1f27..77353a49e2 100644 --- a/tests/multigrid/transfer.cc +++ b/tests/multigrid/transfer.cc @@ -67,7 +67,8 @@ void check_simple(const FiniteElement& fe) template void check_select(const FiniteElement& fe, - std::vector target_component) + std::vector target_component, + std::vector mg_target_component) { deallog << fe.get_name() << std::endl; @@ -79,10 +80,10 @@ void check_select(const FiniteElement& fe, mgdof.distribute_dofs(fe); DoFRenumbering::component_wise(mgdof, target_component); for (unsigned int l=0;l > ndofs(mgdof.get_tria().n_levels()); - MGTools::count_dofs_per_component(mgdof, ndofs, target_component); + MGTools::count_dofs_per_component(mgdof, ndofs, mg_target_component); for (unsigned int l=0;l& fe, MGTransferSelect transfer; transfer.build_matrices(mgdof, 0, - target_component); + mg_target_component); Vector u2(ndofs[2][0]); Vector u1(ndofs[1][0]); @@ -117,6 +118,11 @@ void check_select(const FiniteElement& fe, int main() { + std::ofstream logfile("transfer.output"); + logfile.precision(3); + deallog.attach(logfile); + deallog.depth_console(100); + // check_simple (FE_DGP<2>(0)); // check_simple (FE_DGP<2>(1)); check_simple (FE_DGQ<2>(1)); @@ -134,7 +140,7 @@ int main() v3[i] = i/2; } - check_select (FESystem<2>(FE_DGQ<2>(1), 4), v1); - check_select (FESystem<2>(FE_DGQ<2>(1), 4), v2); - check_select (FESystem<2>(FE_DGQ<2>(1), 4), v3); + check_select (FESystem<2>(FE_DGQ<2>(1), 4), v1, v1); + check_select (FESystem<2>(FE_DGQ<2>(1), 4), v2, v2); + check_select (FESystem<2>(FE_DGQ<2>(1), 4), v3, v3); }