]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Redo output for physics/step-44* tests 3749/head
authorJean-Paul Pelteret <jppelteret@gmail.com>
Sun, 8 Jan 2017 21:36:37 +0000 (22:36 +0100)
committerJean-Paul Pelteret <jppelteret@gmail.com>
Sun, 8 Jan 2017 21:36:37 +0000 (22:36 +0100)
Instead of checking the convergence rates for each Newton iteration
against some benchmark result that is machine specific, we now monitor
the solution at a fixed point. This should remain unchanged if the
problem converges correctly for each timestep. The nonlinear solver
output is now send to the file "stdout".

Fixes #3747

tests/physics/step-44-standard_tensors-material_push_forward.cc
tests/physics/step-44-standard_tensors-material_push_forward.output
tests/physics/step-44-standard_tensors-spatial.cc
tests/physics/step-44-standard_tensors-spatial.output
tests/physics/step-44.cc
tests/physics/step-44.output

index 5fec894681e8e6ed2d57dfe2f23c0f6dc4139299..f9c8eaff2f6f6dc1ea6def7398ba53d4d8ac6ba4 100644 (file)
@@ -807,7 +807,7 @@ namespace Step44
     parameters(input_file),
     triangulation(Triangulation<dim>::maximum_smoothing),
     time(parameters.end_time, parameters.delta_t),
-    timer(deallog.get_file_stream(),
+    timer(std::cout,
           TimerOutput::never,
           TimerOutput::wall_times),
     degree(parameters.poly_degree),
@@ -824,7 +824,7 @@ namespace Step44
     qf_face(parameters.quad_order),
     n_q_points (qf_cell.size()),
     n_q_points_f (qf_face.size()),
-    pcout(deallog.get_file_stream())
+    pcout(std::cout)
   {
     Assert(dim==2 || dim==3, ExcMessage("This problem only works in 2 or 3 space dimensions."));
     determine_component_extractors();
@@ -850,7 +850,7 @@ namespace Step44
                             J_mask,
                             solution_n);
     }
-    output_results();
+    // output_results();
     time.increment();
     BlockVector<double> solution_delta(dofs_per_block);
     while (time.current() < time.end())
@@ -858,7 +858,24 @@ namespace Step44
         solution_delta = 0.0;
         solve_nonlinear_timestep(solution_delta);
         solution_n += solution_delta;
-        output_results();
+        // output_results();
+        // Output displacement at centre of traction surface
+        {
+          const Point<dim> soln_pt (dim==3 ? Point<dim>(0.0, 1.0, 0.0)*parameters.scale : Point<dim>(0.0, 1.0)*parameters.scale);
+          typename DoFHandler<dim>::active_cell_iterator cell =
+            dof_handler_ref.begin_active(), endc = dof_handler_ref.end();
+          for (; cell != endc; ++cell)
+            for (unsigned int v=0; v<GeometryInfo<dim>::vertices_per_cell; ++v)
+              if (cell->vertex(v).distance(soln_pt) < 1e-6*parameters.scale)
+                {
+                  Tensor<1,dim> soln;
+                  for (unsigned int d=0; d<dim; ++d)
+                    soln[d] = solution_n(cell->vertex_dof_index(v,u_dof+d));
+                  deallog
+                      << "Timestep " << time.get_timestep()
+                      << ": " << soln << std::endl;
+                }
+        }
         time.increment();
       }
   }
@@ -1895,7 +1912,7 @@ namespace Step44
                                      * parameters.max_iterations_lin;
               const double tol_sol = parameters.tol_lin
                                      * system_rhs.block(u_dof).l2_norm();
-              SolverControl solver_control(solver_its, tol_sol);
+              SolverControl solver_control(solver_its, tol_sol, false, false);
               GrowingVectorMemory<Vector<double> > GVM;
               SolverCG<Vector<double> > solver_CG(solver_control, GVM);
               PreconditionSelector<SparseMatrix<double>, Vector<double> >
@@ -2045,8 +2062,8 @@ namespace Step44
 int main (int argc,char **argv)
 {
   std::ofstream logfile("output");
+  deallog << std::setprecision(3);
   deallog.attach(logfile);
-  deallog.depth_file(0);
   deallog.threshold_double(1.e-10);
 
   Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv, testing_max_num_threads());
index d0fa5247d0edd3159fd08837006d054a32450d1e..7ed68e5e882e60634c428e199699675afc9fadc6 100644 (file)
 
-Grid:
-        Reference volume: 1.00000e-09
-Triangulation:
-        Number of active cells: 64
-        Number of degrees of freedom: 503
-    Setting up quadrature point data...
-
-Timestep 1 @ 0.100000s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      1.879e-06  1.000e+00  1.000e+00  0.000e+00  0.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       8.427e-04  8.085e-02  8.085e-02  7.076e-13  3.053e-08  9.395e-02  3.598e-02  9.395e-02  9.397e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 122      4.107e-06  8.861e-04  8.861e-04  2.254e-15  1.997e-10  8.585e-04  1.975e-04  8.585e-04  8.586e-04  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      4.885e-11  1.098e-06  1.098e-06  4.211e-20  2.153e-14  9.953e-07  1.498e-07  9.953e-07  9.953e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  1.498e-07
-Force:                 1.303e-11
-Dilatation:    3.484e-07
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 2 @ 2.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      3.594e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       8.305e-04  8.472e-02  8.472e-02  2.288e+13  6.462e+06  1.096e-01  4.023e-02  1.096e-01  1.096e-01  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 114      5.935e-06  9.844e-04  9.844e-04  8.970e+10  6.258e+04  1.154e-03  2.984e-04  1.154e-03  1.154e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 152      8.425e-11  1.584e-06  1.584e-06  2.131e+06  8.412e+00  2.765e-06  4.710e-07  2.765e-06  2.765e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  4.710e-07
-Force:                 2.247e-11
-Dilatation:    7.035e-07
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 3 @ 3.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      2.941e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 92       1.169e-03  9.356e-02  9.356e-02  4.455e+12  4.460e+06  1.131e-01  4.961e-02  1.131e-01  1.131e-01  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 111      3.956e-06  9.355e-04  9.355e-04  2.153e+10  4.528e+04  1.417e-03  4.916e-04  1.417e-03  1.417e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 147      2.074e-10  1.047e-06  1.047e-06  1.149e+06  7.409e+00  1.115e-06  2.439e-07  1.115e-06  1.115e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  2.439e-07
-Force:                 5.530e-11
-Dilatation:    1.059e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 4 @ 4.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      2.007e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 91       1.042e-03  1.069e-01  1.069e-01  8.138e+12  4.014e+06  9.761e-02  6.072e-02  9.761e-02  9.762e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 110      6.023e-06  1.449e-03  1.449e-03  4.399e+10  2.545e+04  1.897e-03  8.054e-04  1.897e-03  1.897e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 143      9.161e-10  1.660e-06  1.660e-06  6.590e+06  1.212e+01  1.847e-06  7.736e-07  1.847e-06  1.848e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  7.736e-07
-Force:                 2.443e-10
-Dilatation:    1.407e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 5 @ 5.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 129      2.191e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 94       6.965e-04  1.179e-01  1.179e-01  1.267e+12  3.437e+06  7.328e-02  6.923e-02  7.328e-02  7.329e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 114      4.389e-06  2.309e-03  2.309e-03  7.096e+09  8.901e+03  2.289e-03  1.160e-03  2.289e-03  2.289e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 138      1.580e-09  1.557e-06  1.557e-06  2.454e+06  1.600e+01  1.620e-06  5.738e-07  1.620e-06  1.620e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  5.738e-07
-Force:                 4.213e-10
-Dilatation:    1.742e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 6 @ 6.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      1.414e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 97       6.939e-04  1.230e-01  1.230e-01  2.558e+12  3.978e+06  5.486e-02  7.377e-02  5.486e-02  5.488e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 118      5.820e-06  2.713e-03  2.713e-03  1.385e+10  4.458e+03  2.263e-03  1.405e-03  2.263e-03  2.263e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 142      2.281e-09  2.218e-06  2.218e-06  7.015e+06  1.794e+01  2.396e-06  6.484e-07  2.396e-06  2.397e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.484e-07
-Force:                 6.081e-10
-Dilatation:    2.057e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 7 @ 7.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      2.699e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 93       7.507e-04  1.234e-01  1.234e-01  1.426e+12  3.501e+06  4.453e-02  7.478e-02  4.453e-02  4.455e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 116      4.617e-06  2.691e-03  2.691e-03  7.526e+09  2.451e+03  2.060e-03  1.491e-03  2.060e-03  2.060e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 147      1.854e-09  1.977e-06  1.977e-06  4.103e+06  1.274e+01  2.026e-06  7.866e-07  2.026e-06  2.026e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  7.866e-07
-Force:                 4.943e-10
-Dilatation:    2.353e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 8 @ 8.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 135      3.661e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 89       5.904e-04  1.207e-01  1.207e-01  1.432e+12  3.128e+06  3.820e-02  7.299e-02  3.820e-02  3.822e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 123      1.810e-06  2.483e-03  2.483e-03  7.548e+09  1.749e+03  1.678e-03  1.443e-03  1.678e-03  1.678e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      9.941e-10  1.233e-06  1.233e-06  3.771e+06  7.563e+00  1.039e-06  6.005e-07  1.039e-06  1.039e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.005e-07
-Force:                 2.651e-10
-Dilatation:    2.632e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 9 @ 9.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 131      3.539e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 92       4.407e-04  1.155e-01  1.155e-01  1.759e+12  3.580e+06  3.352e-02  6.924e-02  3.352e-02  3.353e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      1.937e-06  2.235e-03  2.235e-03  8.944e+09  1.778e+03  1.331e-03  1.252e-03  1.331e-03  1.331e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 154      9.421e-10  1.043e-06  1.043e-06  3.554e+06  5.187e+00  1.008e-06  4.678e-07  1.008e-06  1.008e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  4.678e-07
-Force:                 2.512e-10
-Dilatation:    2.899e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 10 @ 1.000e+00s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 137      1.947e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 92       3.365e-04  1.085e-01  1.085e-01  2.423e+12  3.045e+06  3.001e-02  6.447e-02  3.001e-02  3.003e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 129      1.208e-06  1.904e-03  1.904e-03  1.130e+10  1.284e+03  9.595e-04  1.011e-03  9.595e-04  9.595e-04  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 159      4.322e-10  6.904e-07  6.904e-07  3.233e+06  2.368e+00  5.895e-07  3.490e-07  5.895e-07  5.897e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  3.490e-07
-Force:                 1.153e-10
-Dilatation:    3.157e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 1: 0.00 -0.000110 0.00
+DEAL::Timestep 2: 0.00 -0.000220 0.00
+DEAL::Timestep 3: 0.00 -0.000325 0.00
+DEAL::Timestep 4: 0.00 -0.000422 0.00
+DEAL::Timestep 5: 0.00 -0.000509 0.00
+DEAL::Timestep 6: 0.00 -0.000585 0.00
+DEAL::Timestep 7: 0.00 -0.000650 0.00
+DEAL::Timestep 8: 0.00 -0.000705 0.00
+DEAL::Timestep 9: 0.00 -0.000752 0.00
+DEAL::Timestep 10: 0.00 -0.000790 0.00
index fd86f3cab8a8c4e9ac44f6117d82ce8844e43aa9..ac150b1eeb60edd7da553ee21384f1ce3616e602 100644 (file)
@@ -693,7 +693,7 @@ namespace Step44
     parameters(input_file),
     triangulation(Triangulation<dim>::maximum_smoothing),
     time(parameters.end_time, parameters.delta_t),
-    timer(deallog.get_file_stream(),
+    timer(std::cout,
           TimerOutput::never,
           TimerOutput::wall_times),
     degree(parameters.poly_degree),
@@ -710,7 +710,7 @@ namespace Step44
     qf_face(parameters.quad_order),
     n_q_points (qf_cell.size()),
     n_q_points_f (qf_face.size()),
-    pcout(deallog.get_file_stream())
+    pcout(std::cout)
   {
     Assert(dim==2 || dim==3, ExcMessage("This problem only works in 2 or 3 space dimensions."));
     determine_component_extractors();
@@ -736,7 +736,7 @@ namespace Step44
                             J_mask,
                             solution_n);
     }
-    output_results();
+    // output_results();
     time.increment();
     BlockVector<double> solution_delta(dofs_per_block);
     while (time.current() < time.end())
@@ -744,7 +744,24 @@ namespace Step44
         solution_delta = 0.0;
         solve_nonlinear_timestep(solution_delta);
         solution_n += solution_delta;
-        output_results();
+        // output_results();
+        // Output displacement at centre of traction surface
+        {
+          const Point<dim> soln_pt (dim==3 ? Point<dim>(0.0, 1.0, 0.0)*parameters.scale : Point<dim>(0.0, 1.0)*parameters.scale);
+          typename DoFHandler<dim>::active_cell_iterator cell =
+            dof_handler_ref.begin_active(), endc = dof_handler_ref.end();
+          for (; cell != endc; ++cell)
+            for (unsigned int v=0; v<GeometryInfo<dim>::vertices_per_cell; ++v)
+              if (cell->vertex(v).distance(soln_pt) < 1e-6*parameters.scale)
+                {
+                  Tensor<1,dim> soln;
+                  for (unsigned int d=0; d<dim; ++d)
+                    soln[d] = solution_n(cell->vertex_dof_index(v,u_dof+d));
+                  deallog
+                      << "Timestep " << time.get_timestep()
+                      << ": " << soln << std::endl;
+                }
+        }
         time.increment();
       }
   }
@@ -1781,7 +1798,7 @@ namespace Step44
                                      * parameters.max_iterations_lin;
               const double tol_sol = parameters.tol_lin
                                      * system_rhs.block(u_dof).l2_norm();
-              SolverControl solver_control(solver_its, tol_sol);
+              SolverControl solver_control(solver_its, tol_sol, false, false);
               GrowingVectorMemory<Vector<double> > GVM;
               SolverCG<Vector<double> > solver_CG(solver_control, GVM);
               PreconditionSelector<SparseMatrix<double>, Vector<double> >
@@ -1931,8 +1948,8 @@ namespace Step44
 int main (int argc,char **argv)
 {
   std::ofstream logfile("output");
+  deallog << std::setprecision(3);
   deallog.attach(logfile);
-  deallog.depth_file(0);
   deallog.threshold_double(1.e-10);
 
   Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv, testing_max_num_threads());
index 546f183a9421e03c2078f481b80e4b9d001ddf37..7ed68e5e882e60634c428e199699675afc9fadc6 100644 (file)
 
-Grid:
-        Reference volume: 1.00000e-09
-Triangulation:
-        Number of active cells: 64
-        Number of degrees of freedom: 503
-    Setting up quadrature point data...
-
-Timestep 1 @ 0.100000s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      1.879e-06  1.000e+00  1.000e+00  0.000e+00  0.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       1.114e-03  8.085e-02  8.085e-02  7.076e-13  3.053e-08  9.396e-02  3.598e-02  9.396e-02  9.397e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 123      2.213e-06  9.284e-04  9.284e-04  2.254e-15  1.996e-10  8.786e-04  1.975e-04  8.786e-04  8.788e-04  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      3.555e-11  5.917e-07  5.917e-07  4.211e-20  2.312e-14  6.259e-07  1.495e-07  6.259e-07  6.259e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  1.495e-07
-Force:                 9.480e-12
-Dilatation:    3.484e-07
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 2 @ 2.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      1.864e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 92       1.055e-03  8.472e-02  8.472e-02  2.306e+13  6.107e+06  1.096e-01  4.023e-02  1.096e-01  1.097e-01  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 115      4.858e-06  9.620e-04  9.620e-04  9.042e+10  5.919e+04  1.142e-03  3.104e-04  1.142e-03  1.142e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      8.675e-11  1.296e-06  1.296e-06  2.300e+06  7.668e+00  2.098e-06  4.733e-07  2.098e-06  2.098e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  4.733e-07
-Force:                 2.313e-11
-Dilatation:    7.035e-07
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 3 @ 3.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      2.213e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 95       1.501e-03  9.356e-02  9.356e-02  4.481e+12  4.677e+06  1.130e-01  4.961e-02  1.130e-01  1.131e-01  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 110      7.323e-06  1.013e-03  1.013e-03  2.165e+10  4.746e+04  1.479e-03  4.834e-04  1.479e-03  1.480e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 143      2.610e-10  1.948e-06  1.948e-06  1.157e+06  7.005e+00  2.141e-06  2.425e-07  2.141e-06  2.142e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  2.425e-07
-Force:                 6.959e-11
-Dilatation:    1.059e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 4 @ 4.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      1.859e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 91       1.185e-03  1.069e-01  1.069e-01  8.200e+12  3.902e+06  9.765e-02  6.072e-02  9.765e-02  9.767e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 110      6.518e-06  1.440e-03  1.440e-03  4.432e+10  2.477e+04  1.897e-03  8.054e-04  1.897e-03  1.897e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 143      9.331e-10  1.764e-06  1.764e-06  6.641e+06  1.187e+01  2.206e-06  7.739e-07  2.206e-06  2.207e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  7.739e-07
-Force:                 2.488e-10
-Dilatation:    1.407e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 5 @ 5.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 129      3.277e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 94       7.466e-04  1.179e-01  1.179e-01  1.282e+12  3.596e+06  7.328e-02  6.923e-02  7.328e-02  7.329e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 113      6.856e-06  2.300e-03  2.300e-03  7.177e+09  9.311e+03  2.297e-03  1.160e-03  2.297e-03  2.297e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 140      1.920e-09  2.134e-06  2.134e-06  2.482e+06  1.673e+01  5.052e-06  5.744e-07  5.052e-06  5.053e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  5.744e-07
-Force:                 5.119e-10
-Dilatation:    1.742e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 6 @ 6.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      1.418e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 98       5.446e-04  1.230e-01  1.230e-01  2.598e+12  3.341e+06  5.486e-02  7.377e-02  5.486e-02  5.488e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 118      4.807e-06  2.689e-03  2.689e-03  1.406e+10  3.742e+03  2.236e-03  1.405e-03  2.236e-03  2.236e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 144      3.117e-09  1.929e-06  1.929e-06  7.125e+06  1.493e+01  2.063e-06  6.477e-07  2.063e-06  2.063e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.477e-07
-Force:                 8.310e-10
-Dilatation:    2.057e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 7 @ 7.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      2.877e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 93       7.438e-04  1.234e-01  1.234e-01  1.458e+12  3.481e+06  4.451e-02  7.478e-02  4.451e-02  4.453e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 124      1.896e-06  2.709e-03  2.709e-03  7.695e+09  2.434e+03  2.080e-03  1.491e-03  2.080e-03  2.080e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 147      1.532e-09  1.449e-06  1.449e-06  4.195e+06  1.295e+01  1.160e-06  6.105e-07  1.160e-06  1.160e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.105e-07
-Force:                 4.084e-10
-Dilatation:    2.353e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 8 @ 8.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      3.355e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       4.568e-04  1.207e-01  1.207e-01  1.347e+12  2.889e+06  3.818e-02  7.299e-02  3.818e-02  3.820e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 124      1.981e-06  2.507e-03  2.507e-03  7.100e+09  1.614e+03  1.731e-03  1.428e-03  1.731e-03  1.731e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      9.525e-10  1.147e-06  1.147e-06  3.504e+06  7.450e+00  1.053e-06  6.004e-07  1.053e-06  1.053e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.004e-07
-Force:                 2.540e-10
-Dilatation:    2.632e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 9 @ 9.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 138      2.296e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       3.492e-04  1.155e-01  1.155e-01  1.888e+12  3.080e+06  3.352e-02  6.924e-02  3.352e-02  3.353e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 129      1.002e-06  2.234e-03  2.234e-03  9.602e+09  1.530e+03  1.336e-03  1.247e-03  1.336e-03  1.336e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 153      6.417e-10  8.777e-07  8.777e-07  3.805e+06  4.375e+00  6.542e-07  4.685e-07  6.542e-07  6.543e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  4.685e-07
-Force:                 1.711e-10
-Dilatation:    2.899e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 10 @ 1.000e+00s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 137      2.024e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       2.728e-04  1.085e-01  1.085e-01  2.721e+12  3.072e+06  3.002e-02  6.447e-02  3.002e-02  3.004e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 129      1.029e-06  1.897e-03  1.897e-03  1.269e+10  1.296e+03  9.517e-04  1.008e-03  9.517e-04  9.517e-04  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 159      3.020e-10  6.524e-07  6.524e-07  3.633e+06  2.365e+00  5.342e-07  3.493e-07  5.342e-07  5.344e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  3.493e-07
-Force:                 8.052e-11
-Dilatation:    3.157e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 1: 0.00 -0.000110 0.00
+DEAL::Timestep 2: 0.00 -0.000220 0.00
+DEAL::Timestep 3: 0.00 -0.000325 0.00
+DEAL::Timestep 4: 0.00 -0.000422 0.00
+DEAL::Timestep 5: 0.00 -0.000509 0.00
+DEAL::Timestep 6: 0.00 -0.000585 0.00
+DEAL::Timestep 7: 0.00 -0.000650 0.00
+DEAL::Timestep 8: 0.00 -0.000705 0.00
+DEAL::Timestep 9: 0.00 -0.000752 0.00
+DEAL::Timestep 10: 0.00 -0.000790 0.00
index ff475be415d5fbb63c8de4411d9aeb84ecebbd59..cc2e0653ba4d788996d9d4a045d2bf4766efc1d3 100644 (file)
@@ -710,7 +710,7 @@ namespace Step44
     parameters(input_file),
     triangulation(Triangulation<dim>::maximum_smoothing),
     time(parameters.end_time, parameters.delta_t),
-    timer(deallog.get_file_stream(),
+    timer(std::cout,
           TimerOutput::never,
           TimerOutput::wall_times),
     degree(parameters.poly_degree),
@@ -727,7 +727,7 @@ namespace Step44
     qf_face(parameters.quad_order),
     n_q_points (qf_cell.size()),
     n_q_points_f (qf_face.size()),
-    pcout(deallog.get_file_stream())
+    pcout(std::cout)
   {
     Assert(dim==2 || dim==3, ExcMessage("This problem only works in 2 or 3 space dimensions."));
     determine_component_extractors();
@@ -753,7 +753,7 @@ namespace Step44
                             J_mask,
                             solution_n);
     }
-    output_results();
+    // output_results();
     time.increment();
     BlockVector<double> solution_delta(dofs_per_block);
     while (time.current() < time.end())
@@ -761,7 +761,24 @@ namespace Step44
         solution_delta = 0.0;
         solve_nonlinear_timestep(solution_delta);
         solution_n += solution_delta;
-        output_results();
+        // output_results();
+        // Output displacement at centre of traction surface
+        {
+          const Point<dim> soln_pt (dim==3 ? Point<dim>(0.0, 1.0, 0.0)*parameters.scale : Point<dim>(0.0, 1.0)*parameters.scale);
+          typename DoFHandler<dim>::active_cell_iterator cell =
+            dof_handler_ref.begin_active(), endc = dof_handler_ref.end();
+          for (; cell != endc; ++cell)
+            for (unsigned int v=0; v<GeometryInfo<dim>::vertices_per_cell; ++v)
+              if (cell->vertex(v).distance(soln_pt) < 1e-6*parameters.scale)
+                {
+                  Tensor<1,dim> soln;
+                  for (unsigned int d=0; d<dim; ++d)
+                    soln[d] = solution_n(cell->vertex_dof_index(v,u_dof+d));
+                  deallog
+                      << "Timestep " << time.get_timestep()
+                      << ": " << soln << std::endl;
+                }
+        }
         time.increment();
       }
   }
@@ -1798,7 +1815,7 @@ namespace Step44
                                      * parameters.max_iterations_lin;
               const double tol_sol = parameters.tol_lin
                                      * system_rhs.block(u_dof).l2_norm();
-              SolverControl solver_control(solver_its, tol_sol);
+              SolverControl solver_control(solver_its, tol_sol, false, false);
               GrowingVectorMemory<Vector<double> > GVM;
               SolverCG<Vector<double> > solver_CG(solver_control, GVM);
               PreconditionSelector<SparseMatrix<double>, Vector<double> >
@@ -1948,8 +1965,8 @@ namespace Step44
 int main (int argc,char **argv)
 {
   std::ofstream logfile("output");
+  deallog << std::setprecision(3);
   deallog.attach(logfile);
-  deallog.depth_file(0);
   deallog.threshold_double(1.e-10);
 
   Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv, testing_max_num_threads());
index a5bd04fc20037d5b098647796e6f132bf65ffb98..7ed68e5e882e60634c428e199699675afc9fadc6 100644 (file)
 
-Grid:
-        Reference volume: 1.00000e-09
-Triangulation:
-        Number of active cells: 64
-        Number of degrees of freedom: 503
-    Setting up quadrature point data...
-
-Timestep 1 @ 0.100000s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      1.879e-06  1.000e+00  1.000e+00  0.000e+00  0.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       9.785e-04  8.085e-02  8.085e-02  7.076e-13  3.053e-08  9.397e-02  3.598e-02  9.397e-02  9.399e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 119      4.017e-06  8.826e-04  8.826e-04  2.254e-15  1.998e-10  8.493e-04  1.975e-04  8.493e-04  8.494e-04  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 149      7.738e-11  1.069e-06  1.069e-06  4.212e-20  2.171e-14  9.915e-07  2.729e-07  9.915e-07  9.915e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  2.729e-07
-Force:                 2.063e-11
-Dilatation:    3.484e-07
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 2 @ 2.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 131      3.051e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 92       1.172e-03  8.472e-02  8.472e-02  8.757e+12  5.853e+06  1.096e-01  4.023e-02  1.096e-01  1.097e-01  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 111      5.758e-06  9.375e-04  9.375e-04  3.433e+10  5.674e+04  1.112e-03  3.103e-04  1.112e-03  1.112e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      1.003e-10  1.537e-06  1.537e-06  8.734e+05  6.606e+00  1.426e-06  7.544e-07  1.426e-06  1.426e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  7.544e-07
-Force:                 2.675e-11
-Dilatation:    7.035e-07
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 3 @ 3.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 127      2.879e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 91       1.498e-03  9.356e-02  9.356e-02  1.684e+12  4.976e+06  1.131e-01  4.961e-02  1.131e-01  1.131e-01  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 107      6.289e-06  9.500e-04  9.500e-04  8.139e+09  5.059e+04  1.448e-03  4.916e-04  1.448e-03  1.448e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 147      2.200e-10  1.687e-06  1.687e-06  4.346e+05  8.460e+00  1.502e-06  5.180e-07  1.502e-06  1.502e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  5.180e-07
-Force:                 5.866e-11
-Dilatation:    1.059e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 4 @ 4.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 130      1.876e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 91       1.124e-03  1.069e-01  1.069e-01  2.260e+12  4.315e+06  9.763e-02  6.072e-02  9.763e-02  9.764e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 110      6.459e-06  1.409e-03  1.409e-03  1.222e+10  2.736e+04  1.896e-03  8.054e-04  1.896e-03  1.896e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 144      7.557e-10  1.771e-06  1.771e-06  1.831e+06  1.301e+01  1.957e-06  7.733e-07  1.957e-06  1.957e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  7.733e-07
-Force:                 2.015e-10
-Dilatation:    1.407e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 5 @ 5.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 133      1.921e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 94       7.758e-04  1.179e-01  1.179e-01  1.270e+12  3.762e+06  7.329e-02  6.923e-02  7.329e-02  7.330e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 114      6.837e-06  2.265e-03  2.265e-03  7.113e+09  9.735e+03  2.289e-03  1.160e-03  2.289e-03  2.289e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      2.491e-09  2.033e-06  2.033e-06  2.460e+06  1.721e+01  2.811e-06  5.735e-07  2.811e-06  2.811e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  5.735e-07
-Force:                 6.642e-10
-Dilatation:    1.742e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 6 @ 6.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      1.570e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 93       8.254e-04  1.230e-01  1.230e-01  2.574e+12  3.411e+06  5.485e-02  7.377e-02  5.485e-02  5.486e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 118      6.328e-06  2.737e-03  2.737e-03  1.393e+10  3.815e+03  2.295e-03  1.409e-03  2.295e-03  2.295e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 141      1.786e-09  2.421e-06  2.421e-06  7.040e+06  1.550e+01  2.731e-06  6.528e-07  2.731e-06  2.731e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.528e-07
-Force:                 4.763e-10
-Dilatation:    2.057e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 7 @ 7.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      2.942e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 93       7.359e-04  1.234e-01  1.234e-01  1.412e+12  3.061e+06  4.450e-02  7.478e-02  4.450e-02  4.452e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 120      3.420e-06  2.699e-03  2.699e-03  7.451e+09  2.142e+03  2.079e-03  1.491e-03  2.079e-03  2.079e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 148      1.688e-09  1.625e-06  1.625e-06  4.062e+06  1.135e+01  2.104e-06  6.589e-07  2.104e-06  2.105e-06  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.589e-07
-Force:                 4.502e-10
-Dilatation:    2.353e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 8 @ 8.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 137      2.965e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       5.552e-04  1.207e-01  1.207e-01  1.543e+12  3.613e+06  3.817e-02  7.299e-02  3.817e-02  3.818e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 124      1.818e-06  2.512e-03  2.512e-03  8.134e+09  2.018e+03  1.745e-03  1.428e-03  1.745e-03  1.746e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 151      1.101e-09  1.149e-06  1.149e-06  4.014e+06  9.162e+00  9.495e-07  6.011e-07  9.495e-07  9.496e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  6.011e-07
-Force:                 2.936e-10
-Dilatation:    2.632e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 9 @ 9.000e-01s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      3.726e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 96       3.272e-04  1.155e-01  1.155e-01  1.793e+12  3.058e+06  3.352e-02  6.924e-02  3.352e-02  3.354e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 129      1.004e-06  2.232e-03  2.232e-03  9.116e+09  1.519e+03  1.335e-03  1.247e-03  1.335e-03  1.335e-03  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 153      9.268e-10  8.594e-07  8.594e-07  3.613e+06  4.382e+00  6.229e-07  4.688e-07  6.229e-07  6.230e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  4.688e-07
-Force:                 2.472e-10
-Dilatation:    2.899e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
-
-Timestep 10 @ 1.000e+00s
-___________________________________________________________________________________________________________________________________________________________
-                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
-___________________________________________________________________________________________________________________________________________________________
- 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 134      3.218e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
- 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 95       3.148e-04  1.085e-01  1.085e-01  2.377e+12  3.129e+06  3.002e-02  6.447e-02  3.002e-02  3.004e-02  
- 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 132      9.070e-07  1.905e-03  1.905e-03  1.109e+10  1.319e+03  9.672e-04  1.008e-03  9.672e-04  9.672e-04  
- 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 156      5.143e-10  6.087e-07  6.087e-07  3.174e+06  2.465e+00  4.277e-07  3.161e-07  4.277e-07  4.278e-07  
- 4   ASM_R  CONVERGED! 
-___________________________________________________________________________________________________________________________________________________________
-Relative errors:
-Displacement:  3.161e-07
-Force:                 1.372e-10
-Dilatation:    3.157e-06
-v / V_0:       1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 1: 0.00 -0.000110 0.00
+DEAL::Timestep 2: 0.00 -0.000220 0.00
+DEAL::Timestep 3: 0.00 -0.000325 0.00
+DEAL::Timestep 4: 0.00 -0.000422 0.00
+DEAL::Timestep 5: 0.00 -0.000509 0.00
+DEAL::Timestep 6: 0.00 -0.000585 0.00
+DEAL::Timestep 7: 0.00 -0.000650 0.00
+DEAL::Timestep 8: 0.00 -0.000705 0.00
+DEAL::Timestep 9: 0.00 -0.000752 0.00
+DEAL::Timestep 10: 0.00 -0.000790 0.00

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.