From 08758f986728d7e4afe9d4e7175c10c545999a8d Mon Sep 17 00:00:00 2001 From: janssen Date: Tue, 21 Sep 2010 11:20:35 +0000 Subject: [PATCH] use the set_precision in the assembling for the output git-svn-id: https://svn.dealii.org/trunk@22103 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/multigrid/mg_renumbered_02.cc | 59 ++--------------------------- tests/multigrid/mg_renumbered_03.cc | 1 + 2 files changed, 4 insertions(+), 56 deletions(-) diff --git a/tests/multigrid/mg_renumbered_02.cc b/tests/multigrid/mg_renumbered_02.cc index 5b98abc1d1..82d074127e 100644 --- a/tests/multigrid/mg_renumbered_02.cc +++ b/tests/multigrid/mg_renumbered_02.cc @@ -126,25 +126,6 @@ void diff (Vector &diff, const MGDoFHandler &dof, } } -template -void print(const MGDoFHandler &dof, std::vector > &interface_dofs) -{ - const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; - std::vector dof_indices(dofs_per_cell); - for(unsigned int l=0; l::cell_iterator - cell = dof.begin(l); - cell != dof.end(l); ++cell) - { - cell->get_mg_dof_indices(dof_indices); - for(unsigned int i=0; i class OutputCreator : public Subscriptor @@ -188,15 +169,14 @@ class LaplaceProblem private: void setup_system (); - void test (); void test_boundary (); void output_gpl(const MGDoFHandler &dof, MGLevelObject > &v); void refine_local (); Triangulation triangulation; - const MappingQ1 mapping; - FESystem fe; + const MappingQ1 mapping; + FESystem fe; MGDoFHandler mg_dof_handler_renumbered; const unsigned int degree; @@ -263,6 +243,7 @@ LaplaceProblem::output_gpl(const MGDoFHandler &dof, MeshWorker::Assembler::GnuplotPatch assembler; assembler.initialize(dim, quadrature.size(), dim+dof.get_fe().n_components()); + assembler.set_precision(4); for(unsigned int l=0; l::output_gpl(const MGDoFHandler &dof, } } - - template -void LaplaceProblem::test () -{ - typename FunctionMap::type dirichlet_boundary; - ZeroFunction dirichlet_bc(fe.n_components()); - dirichlet_boundary[0] = &dirichlet_bc; - MGTools::make_boundary_list (mg_dof_handler_renumbered, dirichlet_boundary, - boundary_indices_renumbered); - MGTransferPrebuilt > mg_transfer_renumbered; - mg_transfer_renumbered.build_matrices(mg_dof_handler_renumbered, boundary_indices_renumbered); - - Vector test; - test.reinit(mg_dof_handler_renumbered.n_dofs()); - - MGLevelObject > u(0, triangulation.n_levels()-1); - MGLevelObject > d(0, triangulation.n_levels()-1); - - initialize(mg_dof_handler_renumbered, test); - mg_transfer_renumbered.copy_to_mg(mg_dof_handler_renumbered, u, test); - - for(unsigned int l=0; l void LaplaceProblem::test_boundary () { diff --git a/tests/multigrid/mg_renumbered_03.cc b/tests/multigrid/mg_renumbered_03.cc index 709599ae67..669c5ffb94 100644 --- a/tests/multigrid/mg_renumbered_03.cc +++ b/tests/multigrid/mg_renumbered_03.cc @@ -290,6 +290,7 @@ LaplaceProblem::output_gpl(const MGDoFHandler &dof, MeshWorker::Assembler::GnuplotPatch assembler; assembler.initialize(dim, quadrature.size(), dim+dof.get_fe().n_components()); + assembler.set_precision(4); for(unsigned int l=0; l