]> https://gitweb.dealii.org/ - dealii.git/commitdiff
multigrid test for uniform meshes and sufficiently many cells
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Sun, 26 May 2013 21:41:02 +0000 (21:41 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Sun, 26 May 2013 21:41:02 +0000 (21:41 +0000)
git-svn-id: https://svn.dealii.org/trunk@29622 0785d39b-7218-0410-832d-ea1e28bc413d

tests/mpi/multigrid_uniform.cc [new file with mode: 0644]
tests/mpi/multigrid_uniform/ncpu_1/cmp/generic [new file with mode: 0644]
tests/mpi/multigrid_uniform/ncpu_4/cmp/generic [new file with mode: 0644]

diff --git a/tests/mpi/multigrid_uniform.cc b/tests/mpi/multigrid_uniform.cc
new file mode 100644 (file)
index 0000000..f1a75ef
--- /dev/null
@@ -0,0 +1,536 @@
+//---------------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 2008, 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.
+//
+//---------------------------------------------------------------------------
+
+#include "../tests.h"
+
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/base/function.h>
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/utilities.h>
+
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/lac/full_matrix.h>
+#include <deal.II/lac/sparse_matrix.h>
+#include <deal.II/lac/solver_cg.h>
+#include <deal.II/lac/solver_gmres.h>
+#include <deal.II/lac/precondition.h>
+
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_refinement.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/dofs/dof_tools.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_values.h>
+
+#include <deal.II/numerics/vector_tools.h>
+#include <deal.II/numerics/data_out.h>
+#include <deal.II/numerics/error_estimator.h>
+
+#include <deal.II/base/index_set.h>
+#include <deal.II/distributed/tria.h>
+#include <deal.II/distributed/grid_refinement.h>
+
+#include <deal.II/lac/trilinos_vector.h>
+#include <deal.II/lac/trilinos_sparse_matrix.h>
+#include <deal.II/lac/trilinos_precondition.h>
+
+#include <deal.II/multigrid/mg_dof_handler.h>
+#include <deal.II/multigrid/mg_constrained_dofs.h>
+#include <deal.II/multigrid/multigrid.h>
+#include <deal.II/multigrid/mg_transfer.h>
+#include <deal.II/multigrid/mg_tools.h>
+#include <deal.II/multigrid/mg_coarse.h>
+#include <deal.II/multigrid/mg_smoother.h>
+#include <deal.II/multigrid/mg_matrix.h>
+
+// This is C++:
+#include <fstream>
+#include <sstream>
+
+namespace Step50
+{
+  using namespace dealii;
+
+  template <int dim>
+  class LaplaceProblem
+  {
+  public:
+    LaplaceProblem (const unsigned int deg);
+    void run ();
+
+  private:
+    void setup_system ();
+    void assemble_system ();
+    void assemble_multigrid ();
+    void solve ();
+    void refine_grid ();
+    void output_results (const unsigned int cycle) const;
+
+    parallel::distributed::Triangulation<dim>   triangulation;
+    FE_Q<dim>            fe;
+    DoFHandler<dim>    mg_dof_handler;
+
+    typedef TrilinosWrappers::SparseMatrix matrix_t;
+    typedef TrilinosWrappers::MPI::Vector vector_t;
+
+    matrix_t system_matrix;
+
+    ConstraintMatrix     hanging_node_constraints;
+    ConstraintMatrix     constraints;
+
+    vector_t       solution;
+    vector_t       system_rhs;
+
+    const unsigned int degree;
+
+    MGLevelObject<matrix_t> mg_matrices;
+    MGLevelObject<matrix_t> mg_interface_matrices;
+    MGConstrainedDoFs                    mg_constrained_dofs;
+  };
+
+  template <int dim>
+  class Coefficient : public Function<dim>
+  {
+  public:
+    Coefficient () : Function<dim>() {}
+
+    virtual double value (const Point<dim>   &p,
+                          const unsigned int  component = 0) const;
+
+    virtual void value_list (const std::vector<Point<dim> > &points,
+                             std::vector<double>            &values,
+                             const unsigned int              component = 0) const;
+  };
+
+
+
+  template <int dim>
+  double Coefficient<dim>::value (const Point<dim> &p,
+                                  const unsigned int) const
+  {
+    if (p.square() < 0.5*0.5)
+      return 20;
+    else
+      return 1;
+  }
+
+
+
+  template <int dim>
+  void Coefficient<dim>::value_list (const std::vector<Point<dim> > &points,
+                                     std::vector<double>            &values,
+                                     const unsigned int              component) const
+  {
+    const unsigned int n_points = points.size();
+
+    Assert (values.size() == n_points,
+            ExcDimensionMismatch (values.size(), n_points));
+
+    Assert (component == 0,
+            ExcIndexRange (component, 0, 1));
+
+    for (unsigned int i=0; i<n_points; ++i)
+      values[i] = Coefficient<dim>::value (points[i]);
+  }
+
+  template <int dim>
+  LaplaceProblem<dim>::LaplaceProblem (const unsigned int degree)
+    :
+    triangulation (MPI_COMM_WORLD,Triangulation<dim>::
+                   limit_level_difference_at_vertices,
+                   parallel::distributed::Triangulation<dim>::construct_multigrid_hierarchy),
+    fe (degree),
+    mg_dof_handler (triangulation),
+    degree(degree)
+  {
+    if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD)!=0)
+      deallog.depth_console(0);
+  }
+
+
+  template <int dim>
+  void LaplaceProblem<dim>::setup_system ()
+  {
+    mg_dof_handler.distribute_dofs (fe);
+    mg_dof_handler.distribute_mg_dofs (fe);
+    
+    deallog << "Number of degrees of freedom: "
+            << mg_dof_handler.n_dofs();
+
+    for (unsigned int l=0; l<triangulation.n_levels(); ++l)
+      deallog << "   " << 'L' << l << ": "
+              << mg_dof_handler.n_dofs(l);
+    deallog  << std::endl;
+
+
+    //solution.reinit (mg_dof_handler.n_dofs());
+    //system_rhs.reinit (mg_dof_handler.n_dofs());
+    solution.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
+    system_rhs.reinit(mg_dof_handler.locally_owned_dofs(), MPI_COMM_WORLD);
+
+    constraints.clear ();
+    hanging_node_constraints.clear ();
+    DoFTools::make_hanging_node_constraints (mg_dof_handler, hanging_node_constraints);
+    DoFTools::make_hanging_node_constraints (mg_dof_handler, constraints);
+
+    typename FunctionMap<dim>::type      dirichlet_boundary;
+    ZeroFunction<dim>                    homogeneous_dirichlet_bc (1);
+    dirichlet_boundary[0] = &homogeneous_dirichlet_bc;
+    VectorTools::interpolate_boundary_values (static_cast<const DoFHandler<dim>&>(mg_dof_handler),
+                                              dirichlet_boundary,
+                                              constraints);
+    constraints.close ();
+    hanging_node_constraints.close ();
+
+    CompressedSimpleSparsityPattern csp(mg_dof_handler.n_dofs(), mg_dof_handler.n_dofs());
+    DoFTools::make_sparsity_pattern (mg_dof_handler, csp, constraints);
+    system_matrix.reinit (mg_dof_handler.locally_owned_dofs(), csp, MPI_COMM_WORLD, true);
+
+    mg_constrained_dofs.clear();
+    mg_constrained_dofs.initialize(mg_dof_handler, dirichlet_boundary);
+
+    const unsigned int n_levels = triangulation.n_global_levels();
+
+    mg_interface_matrices.resize(0, n_levels-1);
+    mg_interface_matrices.clear ();
+    mg_matrices.resize(0, n_levels-1);
+    mg_matrices.clear ();
+
+    for (unsigned int level=0; level<n_levels; ++level)
+      {
+        CompressedSparsityPattern csp;
+        csp.reinit(mg_dof_handler.n_dofs(level),
+                   mg_dof_handler.n_dofs(level));
+        MGTools::make_sparsity_pattern(mg_dof_handler, csp, level);
+
+        mg_matrices[level].reinit(mg_dof_handler.locally_owned_mg_dofs(level),
+            mg_dof_handler.locally_owned_mg_dofs(level),
+            csp,
+            MPI_COMM_WORLD, true);
+
+        mg_interface_matrices[level].reinit(mg_dof_handler.locally_owned_mg_dofs(level),
+            mg_dof_handler.locally_owned_mg_dofs(level),
+            csp,
+            MPI_COMM_WORLD, true);
+      }
+  }
+
+  template <int dim>
+  void LaplaceProblem<dim>::assemble_system ()
+  {
+    const QGauss<dim>  quadrature_formula(degree+1);
+
+    FEValues<dim> fe_values (fe, quadrature_formula,
+                             update_values    |  update_gradients |
+                             update_quadrature_points  |  update_JxW_values);
+
+    const unsigned int   dofs_per_cell = fe.dofs_per_cell;
+    const unsigned int   n_q_points    = quadrature_formula.size();
+
+    FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
+    Vector<double>       cell_rhs (dofs_per_cell);
+
+    std::vector<unsigned int> local_dof_indices (dofs_per_cell);
+
+    const Coefficient<dim> coefficient;
+    std::vector<double>    coefficient_values (n_q_points);
+
+    typename DoFHandler<dim>::active_cell_iterator
+    cell = mg_dof_handler.begin_active(),
+    endc = mg_dof_handler.end();
+    for (; cell!=endc; ++cell)
+      if (cell->is_locally_owned())
+        {
+          cell_matrix = 0;
+          cell_rhs = 0;
+
+          fe_values.reinit (cell);
+
+          coefficient.value_list (fe_values.get_quadrature_points(),
+                                  coefficient_values);
+
+          for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
+            for (unsigned int i=0; i<dofs_per_cell; ++i)
+              {
+                for (unsigned int j=0; j<dofs_per_cell; ++j)
+                  cell_matrix(i,j) += (coefficient_values[q_point] *
+                                       fe_values.shape_grad(i,q_point) *
+                                       fe_values.shape_grad(j,q_point) *
+                                       fe_values.JxW(q_point));
+
+                cell_rhs(i) += (fe_values.shape_value(i,q_point) *
+                                1.0 *
+                                fe_values.JxW(q_point));
+              }
+
+          cell->get_dof_indices (local_dof_indices);
+          constraints.distribute_local_to_global (cell_matrix, cell_rhs,
+                                                  local_dof_indices,
+                                                  system_matrix, system_rhs);
+        }
+
+    system_matrix.compress(VectorOperation::add);
+    system_rhs.compress(VectorOperation::add);
+  }
+
+  template <int dim>
+  void LaplaceProblem<dim>::assemble_multigrid ()
+  {
+    QGauss<dim>  quadrature_formula(1+degree);
+
+    FEValues<dim> fe_values (fe, quadrature_formula,
+                             update_values   | update_gradients |
+                             update_quadrature_points | update_JxW_values);
+
+    const unsigned int   dofs_per_cell   = fe.dofs_per_cell;
+    const unsigned int   n_q_points      = quadrature_formula.size();
+
+    FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
+
+    std::vector<unsigned int> local_dof_indices (dofs_per_cell);
+
+    const Coefficient<dim> coefficient;
+    std::vector<double>    coefficient_values (n_q_points);
+
+    std::vector<std::vector<bool> > interface_dofs
+      = mg_constrained_dofs.get_refinement_edge_indices ();
+    std::vector<std::vector<bool> > boundary_interface_dofs
+      = mg_constrained_dofs.get_refinement_edge_boundary_indices ();
+
+    std::vector<ConstraintMatrix> boundary_constraints (triangulation.n_levels());
+    std::vector<ConstraintMatrix> boundary_interface_constraints (triangulation.n_levels());
+    for (unsigned int level=0; level<triangulation.n_levels(); ++level)
+      {
+        boundary_constraints[level].add_lines (interface_dofs[level]);
+        boundary_constraints[level].add_lines (mg_constrained_dofs.get_boundary_indices()[level]);
+        boundary_constraints[level].close ();
+
+        boundary_interface_constraints[level]
+        .add_lines (boundary_interface_dofs[level]);
+        boundary_interface_constraints[level].close ();
+      }
+
+    typename DoFHandler<dim>::cell_iterator cell = mg_dof_handler.begin(),
+                                            endc = mg_dof_handler.end();
+
+    for (; cell!=endc; ++cell)
+      if (cell->level_subdomain_id()==triangulation.locally_owned_subdomain())
+        {
+          cell_matrix = 0;
+          fe_values.reinit (cell);
+
+          coefficient.value_list (fe_values.get_quadrature_points(),
+                                  coefficient_values);
+
+          for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
+            for (unsigned int i=0; i<dofs_per_cell; ++i)
+              for (unsigned int j=0; j<dofs_per_cell; ++j)
+                cell_matrix(i,j) += (coefficient_values[q_point] *
+                                     fe_values.shape_grad(i,q_point) *
+                                     fe_values.shape_grad(j,q_point) *
+                                     fe_values.JxW(q_point));
+
+          cell->get_mg_dof_indices (local_dof_indices);
+         
+          boundary_constraints[cell->level()]
+          .distribute_local_to_global (cell_matrix,
+                                       local_dof_indices,
+                                       mg_matrices[cell->level()]);
+
+          for (unsigned int i=0; i<dofs_per_cell; ++i)
+            for (unsigned int j=0; j<dofs_per_cell; ++j)
+              if ( !(interface_dofs[cell->level()][local_dof_indices[i]]==true &&
+                     interface_dofs[cell->level()][local_dof_indices[j]]==false))
+                cell_matrix(i,j) = 0;
+
+          boundary_interface_constraints[cell->level()]
+          .distribute_local_to_global (cell_matrix,
+                                       local_dof_indices,
+                                       mg_interface_matrices[cell->level()]);
+        }
+
+    for (unsigned int i=0;i<triangulation.n_global_levels();++i)
+      {
+        mg_matrices[i].compress(VectorOperation::add);
+        mg_interface_matrices[i].compress(VectorOperation::add);
+      }
+  }
+
+  template <int dim>
+  void LaplaceProblem<dim>::solve ()
+  {
+    MGTransferPrebuilt<vector_t> mg_transfer(hanging_node_constraints, mg_constrained_dofs);
+    mg_transfer.build_matrices(mg_dof_handler);
+
+    matrix_t coarse_matrix;
+    coarse_matrix.copy_from (mg_matrices[0]);
+    SolverControl coarse_solver_control (1000, 1e-10, false, false);
+    SolverGMRES<vector_t> coarse_solver(coarse_solver_control);
+    PreconditionIdentity id;
+    MGCoarseGridLACIteration<SolverGMRES<vector_t>,vector_t> coarse_grid_solver(coarse_solver,
+        coarse_matrix,
+        id);
+
+    typedef TrilinosWrappers::PreconditionJacobi Smoother;
+    MGSmootherPrecondition<matrix_t, Smoother, vector_t> mg_smoother;
+    mg_smoother.initialize(mg_matrices);
+    mg_smoother.set_steps(2);
+    MGMatrix<matrix_t,vector_t> mg_matrix(&mg_matrices);
+    MGMatrix<matrix_t,vector_t> mg_interface_up(&mg_interface_matrices);
+    MGMatrix<matrix_t,vector_t> mg_interface_down(&mg_interface_matrices);
+
+    // Now, we are ready to set up the
+    // V-cycle operator and the
+    // multilevel preconditioner.
+    Multigrid<vector_t > mg(mg_dof_handler,
+                            mg_matrix,
+                            coarse_grid_solver,
+                            mg_transfer,
+                            mg_smoother,
+                            mg_smoother);
+    mg.set_debug(2);
+    mg.set_edge_matrices(mg_interface_down, mg_interface_up);
+
+    PreconditionMG<dim, vector_t, MGTransferPrebuilt<vector_t> >
+    preconditioner(mg_dof_handler, mg, mg_transfer);
+
+    // With all this together, we can finally
+    // get about solving the linear system in
+    // the usual way:
+    SolverControl solver_control (50, 1e-12, false);
+    SolverCG<vector_t>    cg (solver_control);
+
+    solution = 0;
+
+    cg.solve (system_matrix, solution, system_rhs,
+              preconditioner);
+    constraints.distribute (solution);
+  }
+
+  template <int dim>
+  void LaplaceProblem<dim>::refine_grid ()
+  {
+    Vector<float> estimated_error_per_cell (triangulation.n_active_cells());
+
+    KellyErrorEstimator<dim>::estimate (static_cast<DoFHandler<dim>&>(mg_dof_handler),
+                                        QGauss<dim-1>(3),
+                                        typename FunctionMap<dim>::type(),
+                                        solution,
+                                        estimated_error_per_cell);
+    GridRefinement::refine_and_coarsen_fixed_number (triangulation,
+                                                     estimated_error_per_cell,
+                                                     0.3, 0.03);
+    triangulation.execute_coarsening_and_refinement ();
+  }
+
+  template <int dim>
+  void LaplaceProblem<dim>::run ()
+  {
+    for (unsigned int cycle=0; cycle<3; ++cycle)
+      {
+        deallog << "Cycle " << cycle << ':' << std::endl;
+
+        if (cycle == 0)
+          {
+            GridGenerator::hyper_cube (triangulation);
+
+            // static const HyperBallBoundary<dim> boundary;
+            // triangulation.set_boundary (0, boundary);
+
+            triangulation.refine_global (3);
+          }
+        else
+          triangulation.refine_global (1);
+          //refine_grid ();
+
+
+        deallog << "   Number of active cells:       "
+                << triangulation.n_global_active_cells()
+                << std::endl;
+
+       MPI_Barrier(MPI_COMM_WORLD);
+        setup_system ();
+       MPI_Barrier(MPI_COMM_WORLD);
+
+        deallog << "   Number of degrees of freedom: "
+                << mg_dof_handler.n_dofs()
+                << " (by level: ";
+        for (unsigned int level=0; level<triangulation.n_levels(); ++level)
+          deallog << mg_dof_handler.n_dofs(level)
+                  << (level == triangulation.n_levels()-1
+                      ? ")" : ", ");
+        deallog << std::endl;
+
+       MPI_Barrier(MPI_COMM_WORLD);
+        assemble_system ();
+       MPI_Barrier(MPI_COMM_WORLD);
+        assemble_multigrid ();
+       MPI_Barrier(MPI_COMM_WORLD);
+        solve ();
+      }
+  }
+}
+
+
+// @sect3{The main() function}
+//
+// This is again the same function as
+// in step-6:
+int main (int argc, char *argv[])
+{
+  dealii::Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv);
+  mpi_initlog(__FILE__);
+
+  try
+    {
+      using namespace dealii;
+      using namespace Step50;
+
+      LaplaceProblem<2> laplace_problem(1);
+      laplace_problem.run ();
+    }
+  catch (std::exception &exc)
+    {
+      std::cerr << std::endl << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      std::cerr << "Exception on processing: " << std::endl
+                << exc.what() << std::endl
+                << "Aborting!" << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+
+      return 1;
+    }
+  catch (...)
+    {
+      std::cerr << std::endl << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      std::cerr << "Unknown exception!" << std::endl
+                << "Aborting!" << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      return 1;
+    }
+
+  return 0;
+}
diff --git a/tests/mpi/multigrid_uniform/ncpu_1/cmp/generic b/tests/mpi/multigrid_uniform/ncpu_1/cmp/generic
new file mode 100644 (file)
index 0000000..7761fc8
--- /dev/null
@@ -0,0 +1,557 @@
+
+DEAL::Cycle 0:
+DEAL::   Number of active cells:       64
+DEAL::Number of degrees of freedom: 81   L0: 4   L1: 9   L2: 25   L3: 81
+DEAL::   Number of degrees of freedom: 81 (by level: 4, 9, 25, 81)
+DEAL:cg::Starting value 0.109375
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Convergence step 7 value 2.00480e-14
+DEAL::Cycle 1:
+DEAL::   Number of active cells:       256
+DEAL::Number of degrees of freedom: 289   L0: 4   L1: 9   L2: 25   L3: 81   L4: 289
+DEAL::   Number of degrees of freedom: 289 (by level: 4, 9, 25, 81, 289)
+DEAL:cg::Starting value 0.0585938
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Convergence step 8 value 4.39876e-14
+DEAL::Cycle 2:
+DEAL::   Number of active cells:       1024
+DEAL::Number of degrees of freedom: 1089   L0: 4   L1: 9   L2: 25   L3: 81   L4: 289   L5: 1089
+DEAL::   Number of degrees of freedom: 1089 (by level: 4, 9, 25, 81, 289, 1089)
+DEAL:cg::Starting value 0.0302734
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::Convergence step 9 value 6.07087e-14
diff --git a/tests/mpi/multigrid_uniform/ncpu_4/cmp/generic b/tests/mpi/multigrid_uniform/ncpu_4/cmp/generic
new file mode 100644 (file)
index 0000000..7761fc8
--- /dev/null
@@ -0,0 +1,557 @@
+
+DEAL::Cycle 0:
+DEAL::   Number of active cells:       64
+DEAL::Number of degrees of freedom: 81   L0: 4   L1: 9   L2: 25   L3: 81
+DEAL::   Number of degrees of freedom: 81 (by level: 4, 9, 25, 81)
+DEAL:cg::Starting value 0.109375
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Convergence step 7 value 2.00480e-14
+DEAL::Cycle 1:
+DEAL::   Number of active cells:       256
+DEAL::Number of degrees of freedom: 289   L0: 4   L1: 9   L2: 25   L3: 81   L4: 289
+DEAL::   Number of degrees of freedom: 289 (by level: 4, 9, 25, 81, 289)
+DEAL:cg::Starting value 0.0585938
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Convergence step 8 value 4.39876e-14
+DEAL::Cycle 2:
+DEAL::   Number of active cells:       1024
+DEAL::Number of degrees of freedom: 1089   L0: 4   L1: 9   L2: 25   L3: 81   L4: 289   L5: 1089
+DEAL::   Number of degrees of freedom: 1089 (by level: 4, 9, 25, 81, 289, 1089)
+DEAL:cg::Starting value 0.0302734
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::V-cycle entering level 5
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::Residual on      level 5
+DEAL:cg::V-cycle entering level 4
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::Residual on      level 4
+DEAL:cg::V-cycle entering level 3
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::Residual on      level 3
+DEAL:cg::V-cycle entering level 2
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::Residual on      level 2
+DEAL:cg::V-cycle entering level 1
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::Residual on      level 1
+DEAL:cg::V-cycle entering level 0
+DEAL:cg::Coarse level           0
+DEAL:cg::Smoothing on     level 1
+DEAL:cg::V-cycle leaving  level 1
+DEAL:cg::Smoothing on     level 2
+DEAL:cg::V-cycle leaving  level 2
+DEAL:cg::Smoothing on     level 3
+DEAL:cg::V-cycle leaving  level 3
+DEAL:cg::Smoothing on     level 4
+DEAL:cg::V-cycle leaving  level 4
+DEAL:cg::Smoothing on     level 5
+DEAL:cg::V-cycle leaving  level 5
+DEAL:cg::Convergence step 9 value 6.07087e-14

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.