]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix parallel mg test
authorTimo Heister <timo.heister@gmail.com>
Fri, 3 Jan 2014 13:31:55 +0000 (13:31 +0000)
committerTimo Heister <timo.heister@gmail.com>
Fri, 3 Jan 2014 13:31:55 +0000 (13:31 +0000)
git-svn-id: https://svn.dealii.org/trunk@32149 0785d39b-7218-0410-832d-ea1e28bc413d

tests/mpi/step-39.cc
tests/mpi/step-39.with_trilinos=true.mpirun=2.output
tests/mpi/step-39.with_trilinos=true.mpirun=6.output

index 9410d32187cb5911347630ed7c1bfeb78f780899..9736e2e4e801e60151a3f9355d9dcce155cf69ce 100644 (file)
@@ -209,6 +209,8 @@ namespace Step39
       dinfo.value(0) += penalty * (boundary_values[k] - uh[k]) * (boundary_values[k] - uh[k])
                         * fe.JxW(k);
     dinfo.value(0) = std::sqrt(dinfo.value(0));
+    if (!dinfo.cell->is_locally_owned())
+      dinfo.value(0) = 0.0;
   }
 
 
@@ -239,6 +241,10 @@ namespace Step39
       }
     dinfo1.value(0) = std::sqrt(dinfo1.value(0));
     dinfo2.value(0) = dinfo1.value(0);
+    if (!dinfo1.cell->is_locally_owned())
+      dinfo1.value(0) = 0.0;
+    if (!dinfo2.cell->is_locally_owned())
+      dinfo2.value(0) = 0.0;
   }
 
 
@@ -624,16 +630,16 @@ namespace Step39
     out_data.add(est, "cells");
     assembler.initialize(out_data, false);
 
-    FilteredIterator<typename DoFHandler<dim>::active_cell_iterator>
-    begin(IteratorFilters::LocallyOwnedCell(), dof_handler.begin_active());
-    FilteredIterator<typename DoFHandler<dim>::active_cell_iterator>
-    end(IteratorFilters::LocallyOwnedCell(), dof_handler.end());
-
     Estimator<dim> integrator;
+    MeshWorker::LoopControl lctrl;
+    // assemble all faces adjacent to ghost cells to get the full
+    // information for all own cells without communication
+    lctrl.faces_to_ghost = MeshWorker::LoopControl::both;
+    
     MeshWorker::integration_loop<dim, dim> (
-      begin, end,
+        dof_handler.begin_active(), dof_handler.end(),
       dof_info, info_box,
-      integrator, assembler); // TODO: do this on all faces touching ghost cells
+      integrator, assembler, lctrl);
 
     triangulation.load_user_indices(old_user_indices);
     // estimates is a BlockVector<double> (so serial) on each processor
@@ -733,12 +739,12 @@ namespace Step39
           }
 
         deallog << "Triangulation "
-                << triangulation.n_active_cells() << " cells, "
-                << triangulation.n_levels() << " levels" << std::endl;
+                << triangulation.n_global_active_cells() << " cells, "
+                << triangulation.n_global_levels() << " levels" << std::endl;
 
         setup_system();
         deallog << "DoFHandler " << dof_handler.n_dofs() << " dofs, level dofs";
-        for (unsigned int l=0; l<triangulation.n_levels(); ++l)
+        for (unsigned int l=0; l<triangulation.n_global_levels(); ++l)
           deallog << ' ' << dof_handler.n_dofs(l);
         deallog << std::endl;
 
index 97cb925384ee2cc574da70a16d295c80b62996e0..decf2f7cfdc38b3f188c69ff5faaab6994c8b998 100644 (file)
@@ -9,27 +9,27 @@ DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 20.6576
 DEAL:0:cg::Convergence step 15 value 8.11943e-13
-DEAL:0::Estimate 1.03616
+DEAL:0::Estimate 1.10045
 DEAL:0::Step 1
-DEAL:0::Triangulation 21 cells, 3 levels
-DEAL:0::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:0::Triangulation 18 cells, 3 levels
+DEAL:0::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 21.1045
-DEAL:0:cg::Convergence step 16 value 4.50442e-13
-DEAL:0::Estimate 0.847236
+DEAL:0:cg::Convergence step 16 value 8.12709e-13
+DEAL:0::Estimate 1.03046
 DEAL:0::Step 2
-DEAL:0::Triangulation 30 cells, 4 levels
-DEAL:0::DoFHandler 270 dofs, level dofs 27 108 108 108
+DEAL:0::Triangulation 27 cells, 4 levels
+DEAL:0::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 21.3244
-DEAL:0:cg::Convergence step 16 value 8.57299e-13
-DEAL:0::Estimate 0.707708
+DEAL:0:cg::Convergence step 17 value 6.05005e-13
+DEAL:0::Estimate 0.836764
 DEAL:0::Step 3
 DEAL:0::Triangulation 39 cells, 5 levels
 DEAL:0::DoFHandler 351 dofs, level dofs 27 108 144 108 72
@@ -39,7 +39,7 @@ DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 21.4584
 DEAL:0:cg::Convergence step 16 value 3.55455e-13
-DEAL:0::Estimate 0.608314
+DEAL:0::Estimate 0.644552
 DEAL:0::Step 4
 DEAL:0::Triangulation 54 cells, 6 levels
 DEAL:0::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
@@ -49,17 +49,17 @@ DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 21.5437
 DEAL:0:cg::Convergence step 16 value 8.06366e-13
-DEAL:0::Estimate 0.446476
+DEAL:0::Estimate 0.474987
 DEAL:0::Step 5
-DEAL:0::Triangulation 78 cells, 7 levels
-DEAL:0::DoFHandler 729 dofs, level dofs 27 108 288 216 144 108 72
+DEAL:0::Triangulation 84 cells, 7 levels
+DEAL:0::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
-DEAL:0:cg::Starting value 24.9158
-DEAL:0:cg::Convergence step 17 value 6.10713e-13
-DEAL:0::Estimate 0.336695
+DEAL:0:cg::Starting value 25.5435
+DEAL:0:cg::Convergence step 17 value 8.48262e-13
+DEAL:0::Estimate 0.356044
 
 DEAL:1::Element: FE_DGQ<2>(2)
 DEAL:1::Step 0
@@ -71,27 +71,27 @@ DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 20.6576
 DEAL:1:cg::Convergence step 15 value 8.11943e-13
-DEAL:1::Estimate 1.03616
+DEAL:1::Estimate 1.10045
 DEAL:1::Step 1
-DEAL:1::Triangulation 21 cells, 3 levels
-DEAL:1::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:1::Triangulation 18 cells, 3 levels
+DEAL:1::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 21.1045
-DEAL:1:cg::Convergence step 16 value 4.50442e-13
-DEAL:1::Estimate 0.847236
+DEAL:1:cg::Convergence step 16 value 8.12709e-13
+DEAL:1::Estimate 1.03046
 DEAL:1::Step 2
-DEAL:1::Triangulation 30 cells, 4 levels
-DEAL:1::DoFHandler 270 dofs, level dofs 27 108 108 108
+DEAL:1::Triangulation 27 cells, 4 levels
+DEAL:1::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 21.3244
-DEAL:1:cg::Convergence step 16 value 8.57299e-13
-DEAL:1::Estimate 0.707708
+DEAL:1:cg::Convergence step 17 value 6.05005e-13
+DEAL:1::Estimate 0.836764
 DEAL:1::Step 3
 DEAL:1::Triangulation 39 cells, 5 levels
 DEAL:1::DoFHandler 351 dofs, level dofs 27 108 144 108 72
@@ -101,9 +101,9 @@ DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 21.4584
 DEAL:1:cg::Convergence step 16 value 3.55455e-13
-DEAL:1::Estimate 0.608314
+DEAL:1::Estimate 0.644552
 DEAL:1::Step 4
-DEAL:1::Triangulation 51 cells, 6 levels
+DEAL:1::Triangulation 54 cells, 6 levels
 DEAL:1::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
@@ -111,15 +111,15 @@ DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 21.5437
 DEAL:1:cg::Convergence step 16 value 8.06366e-13
-DEAL:1::Estimate 0.446476
+DEAL:1::Estimate 0.474987
 DEAL:1::Step 5
-DEAL:1::Triangulation 72 cells, 7 levels
-DEAL:1::DoFHandler 729 dofs, level dofs 27 108 288 216 144 108 72
+DEAL:1::Triangulation 84 cells, 7 levels
+DEAL:1::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
-DEAL:1:cg::Starting value 24.9158
-DEAL:1:cg::Convergence step 17 value 6.10713e-13
-DEAL:1::Estimate 0.336695
+DEAL:1:cg::Starting value 25.5435
+DEAL:1:cg::Convergence step 17 value 8.48262e-13
+DEAL:1::Estimate 0.356044
 
index 7addc5299b5384fb9965d6e1c930719e58703e42..b25cdce1c68692f877076598d07e49a8fb459c3d 100644 (file)
@@ -9,119 +9,119 @@ DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 20.6576
 DEAL:0:cg::Convergence step 15 value 8.11943e-13
-DEAL:0::Estimate 0.961852
+DEAL:0::Estimate 1.10045
 DEAL:0::Step 1
 DEAL:0::Triangulation 18 cells, 3 levels
-DEAL:0::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:0::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 21.1045
-DEAL:0:cg::Convergence step 16 value 4.50442e-13
-DEAL:0::Estimate 0.773951
+DEAL:0:cg::Convergence step 16 value 8.12709e-13
+DEAL:0::Estimate 1.03046
 DEAL:0::Step 2
-DEAL:0::Triangulation 21 cells, 3 levels
-DEAL:0::DoFHandler 297 dofs, level dofs 27 108 144
+DEAL:0::Triangulation 27 cells, 4 levels
+DEAL:0::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
-DEAL:0:cg::Starting value 21.3494
-DEAL:0:cg::Convergence step 17 value 4.76110e-13
-DEAL:0::Estimate 0.599181
+DEAL:0:cg::Starting value 21.3244
+DEAL:0:cg::Convergence step 17 value 6.05005e-13
+DEAL:0::Estimate 0.836764
 DEAL:0::Step 3
-DEAL:0::Triangulation 30 cells, 4 levels
-DEAL:0::DoFHandler 378 dofs, level dofs 27 108 144 108
+DEAL:0::Triangulation 39 cells, 5 levels
+DEAL:0::DoFHandler 351 dofs, level dofs 27 108 144 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
 DEAL:0:cg::Starting value 21.4584
-DEAL:0:cg::Convergence step 17 value 3.20961e-13
-DEAL:0::Estimate 0.498225
+DEAL:0:cg::Convergence step 16 value 3.55455e-13
+DEAL:0::Estimate 0.644552
 DEAL:0::Step 4
-DEAL:0::Triangulation 36 cells, 4 levels
-DEAL:0::DoFHandler 540 dofs, level dofs 27 108 216 144
+DEAL:0::Triangulation 54 cells, 6 levels
+DEAL:0::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
-DEAL:0:cg::Starting value 21.9802
-DEAL:0:cg::Convergence step 16 value 6.52149e-13
-DEAL:0::Estimate 0.385295
+DEAL:0:cg::Starting value 21.5437
+DEAL:0:cg::Convergence step 16 value 8.06366e-13
+DEAL:0::Estimate 0.474987
 DEAL:0::Step 5
-DEAL:0::Triangulation 51 cells, 5 levels
-DEAL:0::DoFHandler 837 dofs, level dofs 27 108 360 252 144
+DEAL:0::Triangulation 84 cells, 7 levels
+DEAL:0::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:0::Assemble matrix
 DEAL:0::Assemble multilevel matrix
 DEAL:0::Assemble right hand side
 DEAL:0::Solve
-DEAL:0:cg::Starting value 26.2654
-DEAL:0:cg::Convergence step 15 value 5.00033e-13
-DEAL:0::Estimate 0.286622
+DEAL:0:cg::Starting value 25.5435
+DEAL:0:cg::Convergence step 17 value 8.48262e-13
+DEAL:0::Estimate 0.356044
 
 DEAL:1::Element: FE_DGQ<2>(2)
 DEAL:1::Step 0
-DEAL:1::Triangulation 3 cells, 1 levels
-DEAL:1::DoFHandler 108 dofs, level dofs 27
+DEAL:1::Triangulation 12 cells, 2 levels
+DEAL:1::DoFHandler 108 dofs, level dofs 27 108
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 20.6576
 DEAL:1:cg::Convergence step 15 value 8.11943e-13
-DEAL:1::Estimate 0.961852
+DEAL:1::Estimate 1.10045
 DEAL:1::Step 1
-DEAL:1::Triangulation 21 cells, 3 levels
-DEAL:1::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:1::Triangulation 18 cells, 3 levels
+DEAL:1::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 21.1045
-DEAL:1:cg::Convergence step 16 value 4.50442e-13
-DEAL:1::Estimate 0.773951
+DEAL:1:cg::Convergence step 16 value 8.12709e-13
+DEAL:1::Estimate 1.03046
 DEAL:1::Step 2
-DEAL:1::Triangulation 33 cells, 4 levels
-DEAL:1::DoFHandler 297 dofs, level dofs 27 108 144 108
+DEAL:1::Triangulation 27 cells, 4 levels
+DEAL:1::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
-DEAL:1:cg::Starting value 21.3494
-DEAL:1:cg::Convergence step 17 value 4.76110e-13
-DEAL:1::Estimate 0.599181
+DEAL:1:cg::Starting value 21.3244
+DEAL:1:cg::Convergence step 17 value 6.05005e-13
+DEAL:1::Estimate 0.836764
 DEAL:1::Step 3
 DEAL:1::Triangulation 39 cells, 5 levels
-DEAL:1::DoFHandler 378 dofs, level dofs 27 108 144 108 108
+DEAL:1::DoFHandler 351 dofs, level dofs 27 108 144 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
 DEAL:1:cg::Starting value 21.4584
-DEAL:1:cg::Convergence step 17 value 3.20961e-13
-DEAL:1::Estimate 0.498225
+DEAL:1:cg::Convergence step 16 value 3.55455e-13
+DEAL:1::Estimate 0.644552
 DEAL:1::Step 4
-DEAL:1::Triangulation 48 cells, 6 levels
-DEAL:1::DoFHandler 540 dofs, level dofs 27 108 216 144 108 108
+DEAL:1::Triangulation 54 cells, 6 levels
+DEAL:1::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
-DEAL:1:cg::Starting value 21.9802
-DEAL:1:cg::Convergence step 16 value 6.52149e-13
-DEAL:1::Estimate 0.385295
+DEAL:1:cg::Starting value 21.5437
+DEAL:1:cg::Convergence step 16 value 8.06366e-13
+DEAL:1::Estimate 0.474987
 DEAL:1::Step 5
-DEAL:1::Triangulation 69 cells, 7 levels
-DEAL:1::DoFHandler 837 dofs, level dofs 27 108 360 252 144 108 108
+DEAL:1::Triangulation 84 cells, 7 levels
+DEAL:1::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:1::Assemble matrix
 DEAL:1::Assemble multilevel matrix
 DEAL:1::Assemble right hand side
 DEAL:1::Solve
-DEAL:1:cg::Starting value 26.2654
-DEAL:1:cg::Convergence step 15 value 5.00033e-13
-DEAL:1::Estimate 0.286622
+DEAL:1:cg::Starting value 25.5435
+DEAL:1:cg::Convergence step 17 value 8.48262e-13
+DEAL:1::Estimate 0.356044
 
 
 DEAL:2::Element: FE_DGQ<2>(2)
@@ -134,120 +134,120 @@ DEAL:2::Assemble right hand side
 DEAL:2::Solve
 DEAL:2:cg::Starting value 20.6576
 DEAL:2:cg::Convergence step 15 value 8.11943e-13
-DEAL:2::Estimate 0.961852
+DEAL:2::Estimate 1.10045
 DEAL:2::Step 1
-DEAL:2::Triangulation 21 cells, 3 levels
-DEAL:2::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:2::Triangulation 18 cells, 3 levels
+DEAL:2::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:2::Assemble matrix
 DEAL:2::Assemble multilevel matrix
 DEAL:2::Assemble right hand side
 DEAL:2::Solve
 DEAL:2:cg::Starting value 21.1045
-DEAL:2:cg::Convergence step 16 value 4.50442e-13
-DEAL:2::Estimate 0.773951
+DEAL:2:cg::Convergence step 16 value 8.12709e-13
+DEAL:2::Estimate 1.03046
 DEAL:2::Step 2
 DEAL:2::Triangulation 27 cells, 4 levels
-DEAL:2::DoFHandler 297 dofs, level dofs 27 108 144 108
+DEAL:2::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:2::Assemble matrix
 DEAL:2::Assemble multilevel matrix
 DEAL:2::Assemble right hand side
 DEAL:2::Solve
-DEAL:2:cg::Starting value 21.3494
-DEAL:2:cg::Convergence step 17 value 4.76110e-13
-DEAL:2::Estimate 0.599181
+DEAL:2:cg::Starting value 21.3244
+DEAL:2:cg::Convergence step 17 value 6.05005e-13
+DEAL:2::Estimate 0.836764
 DEAL:2::Step 3
-DEAL:2::Triangulation 30 cells, 4 levels
-DEAL:2::DoFHandler 378 dofs, level dofs 27 108 144 108
+DEAL:2::Triangulation 39 cells, 5 levels
+DEAL:2::DoFHandler 351 dofs, level dofs 27 108 144 108 72
 DEAL:2::Assemble matrix
 DEAL:2::Assemble multilevel matrix
 DEAL:2::Assemble right hand side
 DEAL:2::Solve
 DEAL:2:cg::Starting value 21.4584
-DEAL:2:cg::Convergence step 17 value 3.20961e-13
-DEAL:2::Estimate 0.498225
+DEAL:2:cg::Convergence step 16 value 3.55455e-13
+DEAL:2::Estimate 0.644552
 DEAL:2::Step 4
-DEAL:2::Triangulation 36 cells, 4 levels
-DEAL:2::DoFHandler 540 dofs, level dofs 27 108 216 144
+DEAL:2::Triangulation 54 cells, 6 levels
+DEAL:2::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:2::Assemble matrix
 DEAL:2::Assemble multilevel matrix
 DEAL:2::Assemble right hand side
 DEAL:2::Solve
-DEAL:2:cg::Starting value 21.9802
-DEAL:2:cg::Convergence step 16 value 6.52149e-13
-DEAL:2::Estimate 0.385295
+DEAL:2:cg::Starting value 21.5437
+DEAL:2:cg::Convergence step 16 value 8.06366e-13
+DEAL:2::Estimate 0.474987
 DEAL:2::Step 5
-DEAL:2::Triangulation 45 cells, 5 levels
-DEAL:2::DoFHandler 837 dofs, level dofs 27 108 360 252 144
+DEAL:2::Triangulation 84 cells, 7 levels
+DEAL:2::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:2::Assemble matrix
 DEAL:2::Assemble multilevel matrix
 DEAL:2::Assemble right hand side
 DEAL:2::Solve
-DEAL:2:cg::Starting value 26.2654
-DEAL:2:cg::Convergence step 15 value 5.00033e-13
-DEAL:2::Estimate 0.286622
+DEAL:2:cg::Starting value 25.5435
+DEAL:2:cg::Convergence step 17 value 8.48262e-13
+DEAL:2::Estimate 0.356044
 
 
 DEAL:3::Element: FE_DGQ<2>(2)
 DEAL:3::Step 0
-DEAL:3::Triangulation 3 cells, 1 levels
-DEAL:3::DoFHandler 108 dofs, level dofs 27
+DEAL:3::Triangulation 12 cells, 2 levels
+DEAL:3::DoFHandler 108 dofs, level dofs 27 108
 DEAL:3::Assemble matrix
 DEAL:3::Assemble multilevel matrix
 DEAL:3::Assemble right hand side
 DEAL:3::Solve
 DEAL:3:cg::Starting value 20.6576
 DEAL:3:cg::Convergence step 15 value 8.11943e-13
-DEAL:3::Estimate 0.961852
+DEAL:3::Estimate 1.10045
 DEAL:3::Step 1
-DEAL:3::Triangulation 21 cells, 3 levels
-DEAL:3::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:3::Triangulation 18 cells, 3 levels
+DEAL:3::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:3::Assemble matrix
 DEAL:3::Assemble multilevel matrix
 DEAL:3::Assemble right hand side
 DEAL:3::Solve
 DEAL:3:cg::Starting value 21.1045
-DEAL:3:cg::Convergence step 16 value 4.50442e-13
-DEAL:3::Estimate 0.773951
+DEAL:3:cg::Convergence step 16 value 8.12709e-13
+DEAL:3::Estimate 1.03046
 DEAL:3::Step 2
-DEAL:3::Triangulation 33 cells, 4 levels
-DEAL:3::DoFHandler 297 dofs, level dofs 27 108 144 108
+DEAL:3::Triangulation 27 cells, 4 levels
+DEAL:3::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:3::Assemble matrix
 DEAL:3::Assemble multilevel matrix
 DEAL:3::Assemble right hand side
 DEAL:3::Solve
-DEAL:3:cg::Starting value 21.3494
-DEAL:3:cg::Convergence step 17 value 4.76110e-13
-DEAL:3::Estimate 0.599181
+DEAL:3:cg::Starting value 21.3244
+DEAL:3:cg::Convergence step 17 value 6.05005e-13
+DEAL:3::Estimate 0.836764
 DEAL:3::Step 3
 DEAL:3::Triangulation 39 cells, 5 levels
-DEAL:3::DoFHandler 378 dofs, level dofs 27 108 144 108 108
+DEAL:3::DoFHandler 351 dofs, level dofs 27 108 144 108 72
 DEAL:3::Assemble matrix
 DEAL:3::Assemble multilevel matrix
 DEAL:3::Assemble right hand side
 DEAL:3::Solve
 DEAL:3:cg::Starting value 21.4584
-DEAL:3:cg::Convergence step 17 value 3.20961e-13
-DEAL:3::Estimate 0.498225
+DEAL:3:cg::Convergence step 16 value 3.55455e-13
+DEAL:3::Estimate 0.644552
 DEAL:3::Step 4
-DEAL:3::Triangulation 57 cells, 6 levels
-DEAL:3::DoFHandler 540 dofs, level dofs 27 108 216 144 108 108
+DEAL:3::Triangulation 54 cells, 6 levels
+DEAL:3::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:3::Assemble matrix
 DEAL:3::Assemble multilevel matrix
 DEAL:3::Assemble right hand side
 DEAL:3::Solve
-DEAL:3:cg::Starting value 21.9802
-DEAL:3:cg::Convergence step 16 value 6.52149e-13
-DEAL:3::Estimate 0.385295
+DEAL:3:cg::Starting value 21.5437
+DEAL:3:cg::Convergence step 16 value 8.06366e-13
+DEAL:3::Estimate 0.474987
 DEAL:3::Step 5
-DEAL:3::Triangulation 75 cells, 7 levels
-DEAL:3::DoFHandler 837 dofs, level dofs 27 108 360 252 144 108 108
+DEAL:3::Triangulation 84 cells, 7 levels
+DEAL:3::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:3::Assemble matrix
 DEAL:3::Assemble multilevel matrix
 DEAL:3::Assemble right hand side
 DEAL:3::Solve
-DEAL:3:cg::Starting value 26.2654
-DEAL:3:cg::Convergence step 15 value 5.00033e-13
-DEAL:3::Estimate 0.286622
+DEAL:3:cg::Starting value 25.5435
+DEAL:3:cg::Convergence step 17 value 8.48262e-13
+DEAL:3::Estimate 0.356044
 
 
 DEAL:4::Element: FE_DGQ<2>(2)
@@ -260,118 +260,118 @@ DEAL:4::Assemble right hand side
 DEAL:4::Solve
 DEAL:4:cg::Starting value 20.6576
 DEAL:4:cg::Convergence step 15 value 8.11943e-13
-DEAL:4::Estimate 0.961852
+DEAL:4::Estimate 1.10045
 DEAL:4::Step 1
-DEAL:4::Triangulation 21 cells, 3 levels
-DEAL:4::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:4::Triangulation 18 cells, 3 levels
+DEAL:4::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:4::Assemble matrix
 DEAL:4::Assemble multilevel matrix
 DEAL:4::Assemble right hand side
 DEAL:4::Solve
 DEAL:4:cg::Starting value 21.1045
-DEAL:4:cg::Convergence step 16 value 4.50442e-13
-DEAL:4::Estimate 0.773951
+DEAL:4:cg::Convergence step 16 value 8.12709e-13
+DEAL:4::Estimate 1.03046
 DEAL:4::Step 2
-DEAL:4::Triangulation 30 cells, 4 levels
-DEAL:4::DoFHandler 297 dofs, level dofs 27 108 144 108
+DEAL:4::Triangulation 27 cells, 4 levels
+DEAL:4::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:4::Assemble matrix
 DEAL:4::Assemble multilevel matrix
 DEAL:4::Assemble right hand side
 DEAL:4::Solve
-DEAL:4:cg::Starting value 21.3494
-DEAL:4:cg::Convergence step 17 value 4.76110e-13
-DEAL:4::Estimate 0.599181
+DEAL:4:cg::Starting value 21.3244
+DEAL:4:cg::Convergence step 17 value 6.05005e-13
+DEAL:4::Estimate 0.836764
 DEAL:4::Step 3
 DEAL:4::Triangulation 39 cells, 5 levels
-DEAL:4::DoFHandler 378 dofs, level dofs 27 108 144 108 108
+DEAL:4::DoFHandler 351 dofs, level dofs 27 108 144 108 72
 DEAL:4::Assemble matrix
 DEAL:4::Assemble multilevel matrix
 DEAL:4::Assemble right hand side
 DEAL:4::Solve
 DEAL:4:cg::Starting value 21.4584
-DEAL:4:cg::Convergence step 17 value 3.20961e-13
-DEAL:4::Estimate 0.498225
+DEAL:4:cg::Convergence step 16 value 3.55455e-13
+DEAL:4::Estimate 0.644552
 DEAL:4::Step 4
-DEAL:4::Triangulation 48 cells, 6 levels
-DEAL:4::DoFHandler 540 dofs, level dofs 27 108 216 144 108 108
+DEAL:4::Triangulation 54 cells, 6 levels
+DEAL:4::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:4::Assemble matrix
 DEAL:4::Assemble multilevel matrix
 DEAL:4::Assemble right hand side
 DEAL:4::Solve
-DEAL:4:cg::Starting value 21.9802
-DEAL:4:cg::Convergence step 16 value 6.52149e-13
-DEAL:4::Estimate 0.385295
+DEAL:4:cg::Starting value 21.5437
+DEAL:4:cg::Convergence step 16 value 8.06366e-13
+DEAL:4::Estimate 0.474987
 DEAL:4::Step 5
-DEAL:4::Triangulation 66 cells, 7 levels
-DEAL:4::DoFHandler 837 dofs, level dofs 27 108 360 252 144 108 108
+DEAL:4::Triangulation 84 cells, 7 levels
+DEAL:4::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:4::Assemble matrix
 DEAL:4::Assemble multilevel matrix
 DEAL:4::Assemble right hand side
 DEAL:4::Solve
-DEAL:4:cg::Starting value 26.2654
-DEAL:4:cg::Convergence step 15 value 5.00033e-13
-DEAL:4::Estimate 0.286622
+DEAL:4:cg::Starting value 25.5435
+DEAL:4:cg::Convergence step 17 value 8.48262e-13
+DEAL:4::Estimate 0.356044
 
 
 DEAL:5::Element: FE_DGQ<2>(2)
 DEAL:5::Step 0
-DEAL:5::Triangulation 3 cells, 1 levels
-DEAL:5::DoFHandler 108 dofs, level dofs 27
+DEAL:5::Triangulation 12 cells, 2 levels
+DEAL:5::DoFHandler 108 dofs, level dofs 27 108
 DEAL:5::Assemble matrix
 DEAL:5::Assemble multilevel matrix
 DEAL:5::Assemble right hand side
 DEAL:5::Solve
 DEAL:5:cg::Starting value 20.6576
 DEAL:5:cg::Convergence step 15 value 8.11943e-13
-DEAL:5::Estimate 0.961852
+DEAL:5::Estimate 1.10045
 DEAL:5::Step 1
-DEAL:5::Triangulation 15 cells, 3 levels
-DEAL:5::DoFHandler 189 dofs, level dofs 27 108 108
+DEAL:5::Triangulation 18 cells, 3 levels
+DEAL:5::DoFHandler 162 dofs, level dofs 27 108 72
 DEAL:5::Assemble matrix
 DEAL:5::Assemble multilevel matrix
 DEAL:5::Assemble right hand side
 DEAL:5::Solve
 DEAL:5:cg::Starting value 21.1045
-DEAL:5:cg::Convergence step 16 value 4.50442e-13
-DEAL:5::Estimate 0.773951
+DEAL:5:cg::Convergence step 16 value 8.12709e-13
+DEAL:5::Estimate 1.03046
 DEAL:5::Step 2
-DEAL:5::Triangulation 24 cells, 4 levels
-DEAL:5::DoFHandler 297 dofs, level dofs 27 108 144 108
+DEAL:5::Triangulation 27 cells, 4 levels
+DEAL:5::DoFHandler 243 dofs, level dofs 27 108 108 72
 DEAL:5::Assemble matrix
 DEAL:5::Assemble multilevel matrix
 DEAL:5::Assemble right hand side
 DEAL:5::Solve
-DEAL:5:cg::Starting value 21.3494
-DEAL:5:cg::Convergence step 17 value 4.76110e-13
-DEAL:5::Estimate 0.599181
+DEAL:5:cg::Starting value 21.3244
+DEAL:5:cg::Convergence step 17 value 6.05005e-13
+DEAL:5::Estimate 0.836764
 DEAL:5::Step 3
-DEAL:5::Triangulation 33 cells, 5 levels
-DEAL:5::DoFHandler 378 dofs, level dofs 27 108 144 108 108
+DEAL:5::Triangulation 39 cells, 5 levels
+DEAL:5::DoFHandler 351 dofs, level dofs 27 108 144 108 72
 DEAL:5::Assemble matrix
 DEAL:5::Assemble multilevel matrix
 DEAL:5::Assemble right hand side
 DEAL:5::Solve
 DEAL:5:cg::Starting value 21.4584
-DEAL:5:cg::Convergence step 17 value 3.20961e-13
-DEAL:5::Estimate 0.498225
+DEAL:5:cg::Convergence step 16 value 3.55455e-13
+DEAL:5::Estimate 0.644552
 DEAL:5::Step 4
-DEAL:5::Triangulation 42 cells, 6 levels
-DEAL:5::DoFHandler 540 dofs, level dofs 27 108 216 144 108 108
+DEAL:5::Triangulation 54 cells, 6 levels
+DEAL:5::DoFHandler 486 dofs, level dofs 27 108 180 144 108 72
 DEAL:5::Assemble matrix
 DEAL:5::Assemble multilevel matrix
 DEAL:5::Assemble right hand side
 DEAL:5::Solve
-DEAL:5:cg::Starting value 21.9802
-DEAL:5:cg::Convergence step 16 value 6.52149e-13
-DEAL:5::Estimate 0.385295
+DEAL:5:cg::Starting value 21.5437
+DEAL:5:cg::Convergence step 16 value 8.06366e-13
+DEAL:5::Estimate 0.474987
 DEAL:5::Step 5
-DEAL:5::Triangulation 60 cells, 7 levels
-DEAL:5::DoFHandler 837 dofs, level dofs 27 108 360 252 144 108 108
+DEAL:5::Triangulation 84 cells, 7 levels
+DEAL:5::DoFHandler 756 dofs, level dofs 27 108 324 216 144 108 72
 DEAL:5::Assemble matrix
 DEAL:5::Assemble multilevel matrix
 DEAL:5::Assemble right hand side
 DEAL:5::Solve
-DEAL:5:cg::Starting value 26.2654
-DEAL:5:cg::Convergence step 15 value 5.00033e-13
-DEAL:5::Estimate 0.286622
+DEAL:5:cg::Starting value 25.5435
+DEAL:5:cg::Convergence step 17 value 8.48262e-13
+DEAL:5::Estimate 0.356044
 

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.