From: guido Date: Fri, 1 Aug 2003 21:55:41 +0000 (+0000) Subject: first multigrid test X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3bf4d291cc5fa07a9e1f5632764ea40fe584d1c9;p=dealii-svn.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); } diff --git a/tests/results/i686-pc-linux-gnu+gcc2.95/multigrid/transfer.output b/tests/results/i686-pc-linux-gnu+gcc2.95/multigrid/transfer.output new file mode 100644 index 0000000000..c9b7473583 --- /dev/null +++ b/tests/results/i686-pc-linux-gnu+gcc2.95/multigrid/transfer.output @@ -0,0 +1,18 @@ + +DEAL::FE_DGQ<2>(1) +DEAL::u0 2.00 u1 4.00 u2 8.00 u1 16.0 u0 32.0 +DEAL::FESystem<2>[FE_DGQ<2>(1)^4] +DEAL::Level 0 dofs: 4 4 4 4 +DEAL::Level 1 dofs: 16 16 16 16 +DEAL::Level 2 dofs: 64 64 64 64 +DEAL::u0 2.00 u1 4.00 u2 8.00 u1 16.0 u0 32.0 +DEAL::FESystem<2>[FE_DGQ<2>(1)^4] +DEAL::Level 0 dofs: 16 0 0 0 +DEAL::Level 1 dofs: 64 0 0 0 +DEAL::Level 2 dofs: 256 0 0 0 +DEAL::u0 4.00 u1 8.00 u2 16.0 u1 32.0 u0 64.0 +DEAL::FESystem<2>[FE_DGQ<2>(1)^4] +DEAL::Level 0 dofs: 8 8 0 0 +DEAL::Level 1 dofs: 32 32 0 0 +DEAL::Level 2 dofs: 128 128 0 0 +DEAL::u0 2.83 u1 5.66 u2 11.3 u1 22.6 u0 45.3