From 9ebae5a8e16063b67613292fad6ccd3026536886 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Tue, 2 Mar 2010 21:39:35 +0000 Subject: [PATCH] Add test. git-svn-id: https://svn.dealii.org/trunk@20723 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/multigrid/transfer_system_01.cc | 132 ++++++++++++++++++ .../multigrid/transfer_system_01/cmp/generic | 73 ++++++++++ 2 files changed, 205 insertions(+) create mode 100644 tests/multigrid/transfer_system_01.cc create mode 100644 tests/multigrid/transfer_system_01/cmp/generic diff --git a/tests/multigrid/transfer_system_01.cc b/tests/multigrid/transfer_system_01.cc new file mode 100644 index 0000000000..2808b7ac21 --- /dev/null +++ b/tests/multigrid/transfer_system_01.cc @@ -0,0 +1,132 @@ +//---------------------------------------------------------------------------- +// transfer.cc,v 1.13 2005/12/30 16:07:03 guido Exp +// Version: +// +// Copyright (C) 2000 - 2007, 2009, 2010 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------------- + +// test TransferPrebuilt with system elements + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include +#include + +using namespace std; + +template +void +reinit_vector (const dealii::MGDoFHandler &mg_dof, + MGLevelObject > &v) +{ + for (unsigned int level=v.get_minlevel(); + level<=v.get_maxlevel();++level) + { + unsigned int n = mg_dof.n_dofs (level); + v[level].reinit(n); + } + +} + + +template +void +make_matrix (const Transfer &transfer, + const unsigned int high_level, + FullMatrix &matrix) +{ + Vector src (matrix.n()); + Vector dst (matrix.m()); + for (unsigned int i=0; i &m) +{ + for (unsigned int i=0; i +void check (const FiniteElement& fe) +{ + deallog << fe.get_name() << std::endl; + + Triangulation tr; + GridGenerator::hyper_cube(tr); + tr.refine_global(2); + + MGDoFHandler mg_dof_handler(tr); + mg_dof_handler.distribute_dofs(fe); + + DoFRenumbering::component_wise (mg_dof_handler); + for (unsigned int level=0; level > transfer; + transfer.build_matrices(mg_dof_handler); + + FullMatrix prolong_0_1 (mg_dof_handler.n_dofs(1), + mg_dof_handler.n_dofs(0)); + FullMatrix prolong_1_2 (mg_dof_handler.n_dofs(2), + mg_dof_handler.n_dofs(1)); + + deallog << "Level 0->1" << std::endl; + make_matrix (transfer, 1, prolong_0_1); + print_matrix (prolong_0_1); + + deallog << std::endl; + + deallog << "Level 1->2" << std::endl; + make_matrix (transfer, 2, prolong_1_2); + print_matrix (prolong_1_2); +} + + +int main() +{ + std::ofstream logfile("transfer_system_01/output"); + deallog << std::setprecision(4); + deallog.attach(logfile); + deallog.depth_console(0); + deallog.threshold_double(1.e-10); + +//TODO: do in 1d + check (FESystem<2>(FE_Q<2>(1),2)); +} diff --git a/tests/multigrid/transfer_system_01/cmp/generic b/tests/multigrid/transfer_system_01/cmp/generic new file mode 100644 index 0000000000..40097b80da --- /dev/null +++ b/tests/multigrid/transfer_system_01/cmp/generic @@ -0,0 +1,73 @@ + +DEAL::FESystem<2>[FE_Q<2>(1)^2] +DEAL::Level 0->1 +DEAL::1.000 0 0 0 0 0 0 0 +DEAL::0.5000 0.5000 0 0 0 0 0 0 +DEAL::0.5000 0 0.5000 0 0 0 0 0 +DEAL::0.2500 0.2500 0.2500 0.2500 0 0 0 0 +DEAL::0 1.000 0 0 0 0 0 0 +DEAL::0 0.5000 0 0.5000 0 0 0 0 +DEAL::0 0 1.000 0 0 0 0 0 +DEAL::0 0 0.5000 0.5000 0 0 0 0 +DEAL::0 0 0 1.000 0 0 0 0 +DEAL::0 0 0 0 1.000 0 0 0 +DEAL::0 0 0 0 0.5000 0.5000 0 0 +DEAL::0 0 0 0 0.5000 0 0.5000 0 +DEAL::0 0 0 0 0.2500 0.2500 0.2500 0.2500 +DEAL::0 0 0 0 0 1.000 0 0 +DEAL::0 0 0 0 0 0.5000 0 0.5000 +DEAL::0 0 0 0 0 0 1.000 0 +DEAL::0 0 0 0 0 0 0.5000 0.5000 +DEAL::0 0 0 0 0 0 0 1.000 +DEAL:: +DEAL::Level 1->2 +DEAL::1.000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0.5000 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0.5000 0 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0.2500 0.2500 0.2500 0.2500 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 1.000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0.5000 0 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 1.000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0.5000 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 1.000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0.5000 0 0 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0.2500 0 0.2500 0.2500 0.2500 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 1.000 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0.5000 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0.5000 0 0.5000 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 1.000 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0.5000 0 0 0 0.5000 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0.2500 0.2500 0 0 0.2500 0.2500 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0.5000 0 0 0 0.5000 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 1.000 0 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0.5000 0.5000 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 1.000 0 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0.2500 0 0.2500 0 0.2500 0.2500 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0.5000 0 0 0.5000 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0.5000 0.5000 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 1.000 0 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 1.000 0 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0.5000 0.5000 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0.5000 0 0.5000 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0.2500 0.2500 0.2500 0.2500 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 1.000 0 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0.5000 0 0.5000 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 1.000 0 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0.5000 0.5000 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 1.000 0 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0.5000 0 0 0.5000 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0.2500 0 0.2500 0.2500 0.2500 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 1.000 0 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0.5000 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0 0.5000 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.000 0 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0.5000 0 0 0 0.5000 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0.2500 0.2500 0 0 0.2500 0.2500 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0 0 0 0.5000 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.000 0 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0.5000 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.000 0 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0.2500 0 0.2500 0 0.2500 0.2500 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0 0 0.5000 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0.5000 +DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.000 -- 2.39.5