< const unsigned int n_cycles = 8;
---
> const unsigned int n_cycles = 3;
-659c661
-< output_results(cycle);
----
-> // output_results(cycle);
661c663
< computing_timer.print_summary();
---
< for (unsigned int iteration = 0; iteration <= solution.size(); ++iteration)
---
> for (unsigned int iteration = 0; iteration <= 5; ++iteration)
-602c606
-< output_results(iteration);
----
-> // output_results(iteration);
Updating active set...
Size of active set: 13164
Residual of the non-contact part of the system: 9.42452e-05
+ Writing graphical output...
Newton iteration 1
Assembling system...
Updating active set...
Size of active set: 12363
Residual of the non-contact part of the system: 3.9373
+ Writing graphical output...
Newton iteration 2
Assembling system...
Updating active set...
Size of active set: 11584
Residual of the non-contact part of the system: 2.11652
+ Writing graphical output...
Newton iteration 3
Assembling system...
Updating active set...
Size of active set: 10827
Residual of the non-contact part of the system: 1.46846
+ Writing graphical output...
Newton iteration 4
Assembling system...
Updating active set...
Size of active set: 10092
Residual of the non-contact part of the system: 1.11513
+ Writing graphical output...
Newton iteration 5
Assembling system...
Updating active set...
Size of active set: 9381
Residual of the non-contact part of the system: 0.885955
-
+ Writing graphical output...
+
< for (; current_refinement_cycle < n_refinement_cycles;
---
> for (; current_refinement_cycle < 1;
-2149c2155
-< output_results(current_refinement_cycle);
----
-> // output_results(current_refinement_cycle);
2151,2152c2157,2158
< computing_timer.print_summary();
< computing_timer.reset();
Residual of the non-contact part of the system: 8.20097e-14
with a damping parameter alpha = 1
Active set did not change!
+ Writing graphical output... solution_00.pvtu
Contact force = 37.3058
---
> // std::cout << " ..." << solver_control.last_step()
> // << " CG iterations." << std::endl;
-2161,2162c2170,2171
-< if (timestep_number % 200 == 0)
-< output_results();
----
-> // if (timestep_number % 200 == 0)
-> // output_results();
2164c2173,2174
< if (timestep_number % 25 == 0)
---
< , timer(std::cout, TimerOutput::summary, TimerOutput::wall_times)
---
> , timer(std::cout, TimerOutput::never, TimerOutput::wall_times)
-1105c1105
-< output_results();
----
-> // output_results();
1114c1114
< while (time.current() < time.end())
---
> while (time.current() < 0.3)
-1126c1126
-< output_results();
----
-> // output_results();
-1543c1543
-< std::cout << " Setting up quadrature point data..." << std::endl;
----
-> // std::cout << " Setting up quadrature point data..." << std::endl;
1660c1660
< print_conv_header();
---
Triangulation:
Number of active cells: 64
Number of degrees of freedom: 2699
+ Setting up quadrature point data...
Timestep 1 @ 0.1s
CST ASM_SYS ------ SLV PP UQPH CST ASM_SYS ------ SLV PP UQPH --- ASM_SYS ------ SLV PP UQPH --- ASM_SYS ------ SLV PP UQPH --- ASM_SYS CONVERGED!
-768c768
-< output_results(refinement_cycle);
----
-> // output_results(refinement_cycle);
-1011,1012c1011,1012
-< std::cout << " Writing output..." << std::endl;
-< output_results(refinement_cycle);
----
-> std::cout << " Writing output..." << std::endl;
-> // output_results(refinement_cycle);
-905c905
-< output_results(cycle);
----
-> // output_results(cycle);
Cycle 0 of 4
Number of degrees of freedom: 81
Solving system...
+ Writing graphical output...
Error in the L2 norm : 0.00878095
Error in the H1 seminorm : 0.0709546
Error in the broken H2 seminorm: 1.64542
Cycle 1 of 4
Number of degrees of freedom: 289
Solving system...
+ Writing graphical output...
Error in the L2 norm : 0.003515
Error in the H1 seminorm : 0.0217418
Error in the broken H2 seminorm: 0.81215
Cycle 2 of 4
Number of degrees of freedom: 1089
Solving system...
+ Writing graphical output...
Error in the L2 norm : 0.00110394
Error in the H1 seminorm : 0.00610667
Error in the broken H2 seminorm: 0.401517
Cycle 3 of 4
Number of degrees of freedom: 4225
Solving system...
+ Writing graphical output...
Error in the L2 norm : 0.000308433
Error in the H1 seminorm : 0.00162277
Error in the broken H2 seminorm: 0.199356
< , final_time(10.)
---
> , final_time(-5.)
-495,497c495,497
-< data_out.write_vtu_with_pvtu_record(
-< "./", "solution", timestep_number, MPI_COMM_WORLD, 3);
+497c497
<
---
-> // data_out.write_vtu_with_pvtu_record(
-> // "./", "solution", timestep_number, MPI_COMM_WORLD, 3);
> (void)timestep_number;
521,533d520
< {
-91,91c91,92
-< std::cout << " written to " << filename << std::endl << std::endl;
----
-> // std::cout << " written to " << filename << std::endl << std::endl;
-> (void)filename;
-131c132
+131c131
< std::ifstream f("example.msh");
---
> std::ifstream f("../../../source/step-49/example.msh");
dimension: 2
no. of cells: 11904
boundary indicators: 0(624 times) 1(48 times) 2(32 times) 3(32 times)
+ written to grid-1.vtu
+
Mesh info:
dimension: 2
no. of cells: 14
boundary indicators: 0(22 times)
+ written to grid-2.vtu
+
Mesh info:
dimension: 2
no. of cells: 128
boundary indicators: 0(32 times) 1(32 times)
+ written to grid-3.vtu
+
Mesh info:
dimension: 3
no. of cells: 16
boundary indicators: 0(16 times) 1(16 times) 2(8 times) 3(8 times)
+ written to grid-4.vtu
+
Mesh info:
dimension: 2
no. of cells: 28
boundary indicators: 0(32 times)
+ written to grid-5.vtu
+
Mesh info:
dimension: 2
no. of cells: 1600
boundary indicators: 0(160 times)
+ written to grid-6.vtu
+
Mesh info:
dimension: 2
no. of cells: 256
- boundary indicators: 0(64 times)
+ boundary indicators: 0(64 times)
+ written to grid-7.vtu
+