Skip to content

Commit

Permalink
update test
Browse files Browse the repository at this point in the history
  • Loading branch information
@tjhei
tjhei committed Jun 1, 2024
1 parent ebfb354 commit 8c58da9
Show file tree
Hide file tree
Showing 4 changed files with 72 additions and 189 deletions.
2 changes: 1 addition & 1 deletion tests/nonlinear_failure_strategy_cut.cc
View file
Original file line number Diff line number Diff line change
@@ -1 +1 @@
#include "../benchmarks/nonlinear_channel_flow/simple_nonlinear.cc"
#include "../benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/simple_nonlinear.cc"
2 changes: 1 addition & 1 deletion tests/nonlinear_failure_strategy_cut.prm
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@
include $ASPECT_SOURCE_DIR/tests/iterated_advection_and_stokes_residual.prm

set Nonlinear solver failure strategy = cut timestep size
set Max nonlinear iterations = 7
set Max nonlinear iterations = 4

set Dimension = 2
set End time = 20000
249 changes: 66 additions & 183 deletions tests/nonlinear_failure_strategy_cut/screen-output
View file
Original file line number Diff line number Diff line change
@@ -1,254 +1,137 @@
-----------------------------------------------------------------------------
-----------------------------------------------------------------------------
-----------------------------------------------------------------------------

Loading shared library <./libnonlinear_failure_strategy_cut.so>
Loading shared library <./libnonlinear_failure_strategy_cut.debug.so>

-----------------------------------------------------------------------------
-----------------------------------------------------------------------------
Number of active cells: 512 (on 5 levels)
Number of degrees of freedom: 9,141 (4,290+561+2,145+2,145)

*** Timestep 0: t=0 years
*** Timestep 0: t=0 years, dt=0 years
Solving temperature system... 0 iterations.
Skipping porosity composition solve because RHS is zero.
Rebuilding Stokes preconditioner...
Solving Stokes system... 58+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.89643e-16, 0, 1
Solving Stokes system... 57+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.02441e-16, 0, 1
Relative nonlinear residual (total system) after nonlinear iteration 1: 1

Solving temperature system... 0 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 1+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.89643e-16, 0, 5.20807e-13
Relative nonlinear residual (total system) after nonlinear iteration 2: 5.20807e-13
Solving Stokes system... 3+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.02441e-16, 0, 1.00035e-12
Relative nonlinear residual (total system) after nonlinear iteration 2: 1.00035e-12


Postprocessing:

*** Timestep 1: t=13251 years
Solving temperature system... 10 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 55+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.821377, 0, 0.75783
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.821377

Solving temperature system... 10 iterations.
*** Timestep 1: t=13251 years, dt=13251 years
Solving temperature system... 14 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 53+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0977331, 0, 0.113623
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.113623

Solving temperature system... 9 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 49+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0168861, 0, 0.0205037
Relative nonlinear residual (total system) after nonlinear iteration 3: 0.0205037

Solving temperature system... 9 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 46+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00320707, 0, 0.00392882
Relative nonlinear residual (total system) after nonlinear iteration 4: 0.00392882
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.108621, 0, 0.131943
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.131943

Solving temperature system... 8 iterations.
Solving temperature system... 12 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 43+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000633249, 0, 0.000777886
Relative nonlinear residual (total system) after nonlinear iteration 5: 0.000777886
Solving Stokes system... 47+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00443791, 0, 0.0058187
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0058187

Solving temperature system... 8 iterations.
Solving temperature system... 11 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 40+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000127909, 0, 0.000157509
Relative nonlinear residual (total system) after nonlinear iteration 6: 0.000157509
Solving Stokes system... 42+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000368735, 0, 0.00049859
Relative nonlinear residual (total system) after nonlinear iteration 3: 0.00049859

Solving temperature system... 7 iterations.
Solving temperature system... 10 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 37+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.62433e-05, 0, 3.24128e-05
Relative nonlinear residual (total system) after nonlinear iteration 7: 3.24128e-05

ERROR: the solver in the current timestep failed to converge.
Reducing timestep to 2.0908e+11

*** Timestep 1: t=6625.49 years
Solving temperature system... 9 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 54+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.682328, 0, 0.623091
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.682328
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.52488e-05, 0, 3.48705e-05
Relative nonlinear residual (total system) after nonlinear iteration 4: 3.48705e-05

Solving temperature system... 8 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 51+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.067174, 0, 0.0714663
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0714663

Solving temperature system... 7 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 46+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00881208, 0, 0.00950715
Relative nonlinear residual (total system) after nonlinear iteration 3: 0.00950715

Solving temperature system... 6 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 43+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0012324, 0, 0.00132848
Relative nonlinear residual (total system) after nonlinear iteration 4: 0.00132848

Solving temperature system... 6 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 39+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000176949, 0, 0.000190452
Relative nonlinear residual (total system) after nonlinear iteration 5: 0.000190452

Solving temperature system... 5 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 35+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.5804e-05, 0, 2.77522e-05
Relative nonlinear residual (total system) after nonlinear iteration 6: 2.77522e-05

Solving temperature system... 4 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 31+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.80311e-06, 0, 4.08965e-06
Relative nonlinear residual (total system) after nonlinear iteration 7: 4.08965e-06


Postprocessing:

*** Timestep 2: t=19439.3 years
WARNING: The nonlinear solver in the current timestep failed to converge.
Acting according to the parameter 'Nonlinear solver failure strategy'...
Repeating the current time step based on the time stepping manager ...
*** Timestep 1: t=6625.49 years, dt=6625.49 years
Solving temperature system... 11 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 54+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.349648, 0, 0.479292
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.479292

Solving temperature system... 10 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 52+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0186076, 0, 0.0495985
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0495985
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.055211, 0, 0.0682093
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0682093

Solving temperature system... 9 iterations.
Solving temperature system... 10 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 47+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0022294, 0, 0.00863281
Relative nonlinear residual (total system) after nonlinear iteration 3: 0.00863281
Solving Stokes system... 45+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00119069, 0, 0.00153548
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00153548

Solving temperature system... 8 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 44+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000434729, 0, 0.0017423
Relative nonlinear residual (total system) after nonlinear iteration 4: 0.0017423
Solving Stokes system... 38+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 5.20727e-05, 0, 6.84988e-05
Relative nonlinear residual (total system) after nonlinear iteration 3: 6.84988e-05

Solving temperature system... 7 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 41+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 9.08033e-05, 0, 0.000363433
Relative nonlinear residual (total system) after nonlinear iteration 5: 0.000363433
Solving Stokes system... 32+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.82298e-06, 0, 2.43116e-06
Relative nonlinear residual (total system) after nonlinear iteration 4: 2.43116e-06

Solving temperature system... 6 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 38+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.92857e-05, 0, 7.7081e-05
Relative nonlinear residual (total system) after nonlinear iteration 6: 7.7081e-05

Solving temperature system... 6 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 35+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.13709e-06, 0, 1.65133e-05
Relative nonlinear residual (total system) after nonlinear iteration 7: 1.65133e-05
Postprocessing:

ERROR: the solver in the current timestep failed to converge.
Reducing timestep to 2.02182e+11
*** Timestep 2: t=19280.2 years, dt=12654.7 years
Solving temperature system... 12 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 50+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0196072, 0, 0.0239797
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0239797

*** Timestep 2: t=13032.4 years
Solving temperature system... 10 iterations.
Solving temperature system... 11 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 53+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.140486, 0, 0.273969
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.273969
Solving Stokes system... 43+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000910078, 0, 0.00126542
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00126542

Solving temperature system... 9 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 49+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00973167, 0, 0.0248216
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0248216
Solving Stokes system... 37+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.6147e-05, 0, 3.79486e-05
Relative nonlinear residual (total system) after nonlinear iteration 3: 3.79486e-05

Solving temperature system... 8 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 44+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00110007, 0, 0.00289087
Relative nonlinear residual (total system) after nonlinear iteration 3: 0.00289087

Solving temperature system... 7 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 40+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000141653, 0, 0.000369765
Relative nonlinear residual (total system) after nonlinear iteration 4: 0.000369765

Solving temperature system... 6 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 35+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.88893e-05, 0, 4.92024e-05
Relative nonlinear residual (total system) after nonlinear iteration 5: 4.92024e-05

Solving temperature system... 5 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 32+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.56717e-06, 0, 6.69796e-06
Relative nonlinear residual (total system) after nonlinear iteration 6: 6.69796e-06
Solving Stokes system... 30+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.02452e-06, 0, 1.49967e-06
Relative nonlinear residual (total system) after nonlinear iteration 4: 1.49967e-06


Postprocessing:

*** Timestep 3: t=20000 years
Solving temperature system... 10 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 53+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.153077, 0, 0.374788
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.374788

Solving temperature system... 8 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 47+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00697232, 0, 0.020623
Relative nonlinear residual (total system) after nonlinear iteration 2: 0.020623

Solving temperature system... 7 iterations.
*** Timestep 3: t=20000 years, dt=719.818 years
Solving temperature system... 9 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 43+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000414911, 0, 0.00143987
Relative nonlinear residual (total system) after nonlinear iteration 3: 0.00143987
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000709478, 0, 0.000896481
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.000896481

Solving temperature system... 6 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 38+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.212e-05, 0, 0.000119257
Relative nonlinear residual (total system) after nonlinear iteration 4: 0.000119257

Solving temperature system... 5 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 34+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.86401e-06, 0, 1.0762e-05
Relative nonlinear residual (total system) after nonlinear iteration 5: 1.0762e-05

Solving temperature system... 4 iterations.
Skipping porosity composition solve because RHS is zero.
Solving Stokes system... 29+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.72973e-07, 0, 1.02263e-06
Relative nonlinear residual (total system) after nonlinear iteration 6: 1.02263e-06
Solving Stokes system... 32+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.11247e-06, 0, 4.34375e-06
Relative nonlinear residual (total system) after nonlinear iteration 2: 4.34375e-06


Postprocessing:

Termination requested by criterion: end time


+---------------------------------------------+------------+------------+
+---------------------------------+-----------+------------+------------+
+---------------------------------+-----------+------------+------------+
+----------------------------------------------+------------+------------+
+----------------------------------+-----------+------------+------------+
+----------------------------------+-----------+------------+------------+

-----------------------------------------------------------------------------
-----------------------------------------------------------------------------
8 changes: 4 additions & 4 deletions tests/nonlinear_failure_strategy_cut/statistics
View file
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@
# 11: Iterations for Stokes solver
# 12: Velocity iterations in Stokes preconditioner
# 13: Schur complement iterations in Stokes preconditioner
0 0.000000000000e+00 0.000000000000e+00 512 4851 2145 2145 2 0 0 59 62 61
1 6.625491997651e+03 6.625491997651e+03 512 4851 2145 2145 14 106 0 622 640 640
2 1.303237182611e+04 6.406879828463e+03 512 4851 2145 2145 13 102 0 564 580 580
3 2.000000000000e+04 6.967628173886e+03 512 4851 2145 2145 6 40 0 244 251 251
0 0.000000000000e+00 0.000000000000e+00 512 4851 2145 2145 2 0 0 58 60 59
1 6.625491997652e+03 6.625491997652e+03 512 4851 2145 2145 8 83 0 338 346 346
2 1.928018171317e+04 1.265468971551e+04 512 4851 2145 2145 4 40 0 156 160 160
3 2.000000000000e+04 7.198182868337e+02 512 4851 2145 2145 2 15 0 73 75 75

0 comments on commit 8c58da9

Please sign in to comment.