From cbcbdecc30a23537620c4369975c3ce724030191 Mon Sep 17 00:00:00 2001 From: Jean-Paul Pelteret Date: Sun, 8 Jan 2017 22:36:37 +0100 Subject: [PATCH] Redo output for physics/step-44* tests 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 --- ...-standard_tensors-material_push_forward.cc | 29 ++- ...ndard_tensors-material_push_forward.output | 176 +----------------- .../step-44-standard_tensors-spatial.cc | 29 ++- .../step-44-standard_tensors-spatial.output | 176 +----------------- tests/physics/step-44.cc | 29 ++- tests/physics/step-44.output | 176 +----------------- 6 files changed, 99 insertions(+), 516 deletions(-) diff --git a/tests/physics/step-44-standard_tensors-material_push_forward.cc b/tests/physics/step-44-standard_tensors-material_push_forward.cc index 5fec894681..f9c8eaff2f 100644 --- a/tests/physics/step-44-standard_tensors-material_push_forward.cc +++ b/tests/physics/step-44-standard_tensors-material_push_forward.cc @@ -807,7 +807,7 @@ namespace Step44 parameters(input_file), triangulation(Triangulation::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 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 soln_pt (dim==3 ? Point(0.0, 1.0, 0.0)*parameters.scale : Point(0.0, 1.0)*parameters.scale); + typename DoFHandler::active_cell_iterator cell = + dof_handler_ref.begin_active(), endc = dof_handler_ref.end(); + for (; cell != endc; ++cell) + for (unsigned int v=0; v::vertices_per_cell; ++v) + if (cell->vertex(v).distance(soln_pt) < 1e-6*parameters.scale) + { + Tensor<1,dim> soln; + for (unsigned int d=0; dvertex_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 > GVM; SolverCG > solver_CG(solver_control, GVM); PreconditionSelector, Vector > @@ -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()); diff --git a/tests/physics/step-44-standard_tensors-material_push_forward.output b/tests/physics/step-44-standard_tensors-material_push_forward.output index d0fa5247d0..7ed68e5e88 100644 --- a/tests/physics/step-44-standard_tensors-material_push_forward.output +++ b/tests/physics/step-44-standard_tensors-material_push_forward.output @@ -1,167 +1,11 @@ -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 diff --git a/tests/physics/step-44-standard_tensors-spatial.cc b/tests/physics/step-44-standard_tensors-spatial.cc index fd86f3cab8..ac150b1eeb 100644 --- a/tests/physics/step-44-standard_tensors-spatial.cc +++ b/tests/physics/step-44-standard_tensors-spatial.cc @@ -693,7 +693,7 @@ namespace Step44 parameters(input_file), triangulation(Triangulation::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 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 soln_pt (dim==3 ? Point(0.0, 1.0, 0.0)*parameters.scale : Point(0.0, 1.0)*parameters.scale); + typename DoFHandler::active_cell_iterator cell = + dof_handler_ref.begin_active(), endc = dof_handler_ref.end(); + for (; cell != endc; ++cell) + for (unsigned int v=0; v::vertices_per_cell; ++v) + if (cell->vertex(v).distance(soln_pt) < 1e-6*parameters.scale) + { + Tensor<1,dim> soln; + for (unsigned int d=0; dvertex_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 > GVM; SolverCG > solver_CG(solver_control, GVM); PreconditionSelector, Vector > @@ -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()); diff --git a/tests/physics/step-44-standard_tensors-spatial.output b/tests/physics/step-44-standard_tensors-spatial.output index 546f183a94..7ed68e5e88 100644 --- a/tests/physics/step-44-standard_tensors-spatial.output +++ b/tests/physics/step-44-standard_tensors-spatial.output @@ -1,167 +1,11 @@ -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 diff --git a/tests/physics/step-44.cc b/tests/physics/step-44.cc index ff475be415..cc2e0653ba 100644 --- a/tests/physics/step-44.cc +++ b/tests/physics/step-44.cc @@ -710,7 +710,7 @@ namespace Step44 parameters(input_file), triangulation(Triangulation::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 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 soln_pt (dim==3 ? Point(0.0, 1.0, 0.0)*parameters.scale : Point(0.0, 1.0)*parameters.scale); + typename DoFHandler::active_cell_iterator cell = + dof_handler_ref.begin_active(), endc = dof_handler_ref.end(); + for (; cell != endc; ++cell) + for (unsigned int v=0; v::vertices_per_cell; ++v) + if (cell->vertex(v).distance(soln_pt) < 1e-6*parameters.scale) + { + Tensor<1,dim> soln; + for (unsigned int d=0; dvertex_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 > GVM; SolverCG > solver_CG(solver_control, GVM); PreconditionSelector, Vector > @@ -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()); diff --git a/tests/physics/step-44.output b/tests/physics/step-44.output index a5bd04fc20..7ed68e5e88 100644 --- a/tests/physics/step-44.output +++ b/tests/physics/step-44.output @@ -1,167 +1,11 @@ -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 -- 2.39.5