[cig-commits] [commit] master: add a test for the iterated IMPES solver (86d8c86)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Mon Dec 1 12:55:04 PST 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/822ab1757dd06d1edde2754abfa69ba0a0774b9e...66c0a9b8dcd358b4709abe5984f05de5e376e87b
>---------------------------------------------------------------
commit 86d8c86ed0dbc285a7ff2040a7f2a4286caa0364
Author: Juliane Dannberg <dannberg at gfz-potsdam.de>
Date: Mon Nov 17 20:43:06 2014 +0100
add a test for the iterated IMPES solver
>---------------------------------------------------------------
86d8c86ed0dbc285a7ff2040a7f2a4286caa0364
tests/iterated_IMPES.prm | 138 +++++++++
tests/iterated_IMPES/screen-output | 574 +++++++++++++++++++++++++++++++++++++
tests/iterated_IMPES/statistics | 96 +++++++
3 files changed, 808 insertions(+)
diff --git a/tests/iterated_IMPES.prm b/tests/iterated_IMPES.prm
new file mode 100644
index 0000000..5c70a8e
--- /dev/null
+++ b/tests/iterated_IMPES.prm
@@ -0,0 +1,138 @@
+#########################################################
+# This is a model setup to test the iterated IMPES solver.
+# Both the calculation of the residuum and the convergence
+# with prescribed boundary conditions are tested.
+
+set Adiabatic surface temperature = 0 # default: 0
+set CFL number = 0.5
+set Composition solver tolerance = 1e-12
+set Nonlinear solver scheme = iterated IMPES
+set Max nonlinear iterations = 10
+set Linear solver tolerance = 1e-12
+
+# The number of space dimensions you want to run this program in.
+set Dimension = 2
+
+# The end time of the simulation. Units: years if the 'Use years in output
+# instead of seconds' parameter is set; seconds otherwise.
+# This end time is chosen in such a way that the solitary wave travels
+# approximately 5 times its wavelength during the model time.
+set End time = 1e5
+
+set Pressure normalization = surface
+set Surface pressure = 0
+set Resume computation = false
+set Start time = 0
+
+set Use years in output instead of seconds = true
+set Use direct solver for Stokes system = false
+
+subsection Discretization
+ set Stokes velocity polynomial degree = 2
+
+ subsection Stabilization parameters
+ set beta = 0.5
+ end
+end
+
+subsection Compositional fields
+ set Number of fields = 1
+ set Names of fields = porosity
+end
+
+
+subsection Boundary temperature model
+ set Model name = box
+
+ subsection Box
+ set Top temperature = 0 # default: 6000
+ set Bottom temperature = 1 # default: 0
+ end
+end
+
+subsection Boundary composition model
+ set Model name = initial composition
+end
+
+subsection Geometry model
+ set Model name = box # default:
+
+ subsection Box
+ set X extent = 150000
+ set Y extent = 65000
+ set X repetitions = 2
+ end
+
+end
+
+
+subsection Gravity model
+ set Model name = vertical
+ subsection Vertical
+ set Magnitude = 10.0
+ end
+
+end
+
+subsection Initial conditions
+ set Model name = function
+ subsection Function
+ set Function constants = pi=3.1415926,x0=150000,a = 0.2, b = 75000, c = 10000, d=0.0
+ set Function expression = d + a * exp(-((x-b)*(x-b)+y*y)/(2*c*c))
+ set Variable names = x,y
+ end
+end
+
+subsection Compositional initial conditions
+ set Model name = function
+ subsection Function
+ set Function constants = pi=3.1415926,x0=150000,a = 0.2, b = 75000, c = 10000, d=0.0
+ set Function expression = d + a * exp(-((x-b)*(x-b)+y*y)/(2*c*c))
+ set Variable names = x,y
+ end
+end
+
+
+subsection Material model
+ set Model name = simple
+ subsection Simple model
+ set Thermal conductivity = 0
+ set Thermal expansion coefficient = 1
+ set Thermal viscosity exponent = 0.0
+ set Reference specific heat = 1.0
+ set Reference density = 1.0
+ set Reference temperature = 0
+ set Viscosity = 1e16
+ end
+end
+
+
+subsection Mesh refinement
+ set Coarsening fraction = 0.0
+ set Refinement fraction = 0.0
+
+ set Initial adaptive refinement = 0 # default: 2
+ set Initial global refinement = 4 # default: 2
+ set Strategy = composition
+ set Time steps between mesh refinement = 0
+end
+
+
+subsection Model settings
+ set Fixed temperature boundary indicators = 2,3
+ set Fixed composition boundary indicators = 2,3
+ set Prescribed velocity boundary indicators =
+
+ set Tangential velocity boundary indicators = 0,1,2,3
+ set Zero velocity boundary indicators =
+
+ set Include shear heating = false
+end
+
+
+subsection Postprocess
+ set List of postprocessors =
+end
+
+
+
diff --git a/tests/iterated_IMPES/screen-output b/tests/iterated_IMPES/screen-output
new file mode 100644
index 0000000..febcfb6
--- /dev/null
+++ b/tests/iterated_IMPES/screen-output
@@ -0,0 +1,574 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
+-- . version 1.2.pre
+-- . running in OPTIMIZED mode
+-- . running with 1 MPI process
+-- . using Trilinos
+-----------------------------------------------------------------------------
+
+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
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 9.23139e-10, 8.98822e+07, 1.13069e-09
+
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 2 iterations.
+ Nonlinear residuals: 9.23139e-10, 7.27522e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Solving Stokes system... 0 iterations.
+ Nonlinear residuals: 9.23139e-10, 2.73618e-05, 1.13069e-09
+
+ residual: 1
+
+ Postprocessing:
+
+*** Timestep 1: t=13251 years
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 1.30257e+07, 7.76212e+07, 362042
+
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 1.54989e+06, 1.16379e+07, 23807.8
+
+ residual: 0.149931
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30+3 iterations.
+ Nonlinear residuals: 267787, 2.1001e+06, 3251.24
+
+ residual: 0.0270558
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 50859, 402412, 654.641
+
+ residual: 0.0051843
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 10042.3, 79675.5, 148.515
+
+ residual: 0.00102647
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 30+1 iterations.
+ Nonlinear residuals: 2028.44, 16133, 34.3545
+
+ residual: 0.000207843
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 29 iterations.
+ Nonlinear residuals: 416.177, 3319.9, 7.88745
+
+ residual: 4.27706e-05
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 26 iterations.
+ Nonlinear residuals: 86.3598, 691.174, 1.78676
+
+ residual: 8.90445e-06
+
+ Postprocessing:
+
+*** Timestep 2: t=25357.2 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 8.74328e+06, 2.77541e+07, 68345.8
+
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving Stokes system... 30+3 iterations.
+ Nonlinear residuals: 534360, 3.36554e+06, 3864.78
+
+ residual: 0.121263
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 73992, 562190, 679.828
+
+ residual: 0.0202561
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 13458.8, 102450, 185.172
+
+ residual: 0.00369136
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 2554.88, 19327.3, 48.2406
+
+ residual: 0.000705831
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 29 iterations.
+ Nonlinear residuals: 494.904, 3733.66, 11.6979
+
+ residual: 0.000171158
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 27 iterations.
+ Nonlinear residuals: 97.4349, 734.428, 2.7041
+
+ residual: 3.9565e-05
+ Solving temperature system... 5 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 24 iterations.
+ Nonlinear residuals: 19.4318, 146.478, 0.605517
+
+ residual: 8.85961e-06
+
+ Postprocessing:
+
+*** Timestep 3: t=36192.1 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 8.78636e+06, 3.96054e+07, 52145.3
+
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 298794, 1.66045e+06, 3943.3
+
+ residual: 0.0756213
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 13880.1, 90098, 184.183
+
+ residual: 0.0035321
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 30+1 iterations.
+ Nonlinear residuals: 1194.47, 8793.38, 24.9373
+
+ residual: 0.000478228
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 27 iterations.
+ Nonlinear residuals: 150.623, 1120.11, 2.9169
+
+ residual: 5.5938e-05
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 24 iterations.
+ Nonlinear residuals: 20.2838, 150.152, 0.314073
+
+ residual: 6.02303e-06
+
+ Postprocessing:
+
+*** Timestep 4: t=46125 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 3.47276e+06, 3.10297e+07, 55955.2
+
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 121429, 787205, 7098.88
+
+ residual: 0.126867
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 7066.23, 54337, 552.304
+
+ residual: 0.00987046
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30+1 iterations.
+ Nonlinear residuals: 662.283, 5212.76, 56.8192
+
+ residual: 0.00101544
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 27 iterations.
+ Nonlinear residuals: 68.9605, 551.256, 6.20146
+
+ residual: 0.000110829
+ Solving temperature system... 5 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 23 iterations.
+ Nonlinear residuals: 7.34032, 58.7242, 0.636452
+
+ residual: 1.13743e-05
+ Solving temperature system... 4 iterations.
+ Solving composition system 1... 4 iterations.
+ Solving Stokes system... 20 iterations.
+ Nonlinear residuals: 0.788769, 6.30559, 0.0642126
+
+ residual: 1.14757e-06
+
+ Postprocessing:
+
+*** Timestep 5: t=55855.4 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 2.20444e+06, 1.81351e+07, 56884.3
+
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 71532.3, 568086, 8058.51
+
+ residual: 0.141665
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 5106.63, 43981.5, 1031.27
+
+ residual: 0.0181293
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30 iterations.
+ Nonlinear residuals: 482.133, 4266.08, 107.856
+
+ residual: 0.00189606
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 27 iterations.
+ Nonlinear residuals: 47.9479, 427.409, 10.9268
+
+ residual: 0.000192088
+ Solving temperature system... 5 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 23 iterations.
+ Nonlinear residuals: 4.82082, 43.0829, 1.08872
+
+ residual: 1.91392e-05
+ Solving temperature system... 4 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 19 iterations.
+ Nonlinear residuals: 0.488538, 4.36167, 0.107761
+
+ residual: 1.89438e-06
+
+ Postprocessing:
+
+*** Timestep 6: t=65167.2 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 1.47905e+06, 1.27436e+07, 60360.2
+
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 45392.5, 397633, 9396.55
+
+ residual: 0.155675
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 3292.3, 31095.3, 987.531
+
+ residual: 0.0163606
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30 iterations.
+ Nonlinear residuals: 305.891, 2975.95, 98.0598
+
+ residual: 0.00162458
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 26 iterations.
+ Nonlinear residuals: 29.6601, 293.924, 9.63021
+
+ residual: 0.000159546
+ Solving temperature system... 5 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 22 iterations.
+ Nonlinear residuals: 2.88929, 29.0718, 0.944393
+
+ residual: 1.5646e-05
+ Solving temperature system... 4 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 19 iterations.
+ Nonlinear residuals: 0.281158, 2.86638, 0.0926398
+
+ residual: 1.53478e-06
+
+ Postprocessing:
+
+*** Timestep 7: t=73771.7 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 940943, 8.43616e+06, 64026.9
+
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 29828.1, 290742, 8702.22
+
+ residual: 0.135915
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 2227.04, 23199.9, 826.527
+
+ residual: 0.0129091
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 29 iterations.
+ Nonlinear residuals: 202.161, 2177.94, 77.3223
+
+ residual: 0.00120765
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 25 iterations.
+ Nonlinear residuals: 19.0296, 209.277, 7.2626
+
+ residual: 0.00011343
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 22 iterations.
+ Nonlinear residuals: 1.7984, 20.0627, 0.685094
+
+ residual: 1.07001e-05
+ Solving temperature system... 4 iterations.
+ Solving composition system 1... 4 iterations.
+ Solving Stokes system... 18 iterations.
+ Nonlinear residuals: 0.169731, 1.91361, 0.064782
+
+ residual: 1.01179e-06
+
+ Postprocessing:
+
+*** Timestep 8: t=81630.4 years
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 595027, 5.46878e+06, 66651.8
+
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 21593.1, 227819, 7503.19
+
+ residual: 0.112573
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 1686, 18497.5, 666.472
+
+ residual: 0.00999931
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 29 iterations.
+ Nonlinear residuals: 148.219, 1669.58, 58.9207
+
+ residual: 0.000884007
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 25 iterations.
+ Nonlinear residuals: 13.333, 152.215, 5.25042
+
+ residual: 7.87738e-05
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 21 iterations.
+ Nonlinear residuals: 1.20102, 13.8173, 0.469926
+
+ residual: 7.05047e-06
+
+ Postprocessing:
+
+*** Timestep 9: t=88740.3 years
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 391928, 3.62238e+06, 67395.2
+
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 16341.3, 180511, 6246.83
+
+ residual: 0.0926896
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 1276.77, 14484.7, 509.505
+
+ residual: 0.00755997
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 28 iterations.
+ Nonlinear residuals: 106.112, 1230.68, 41.9683
+
+ residual: 0.00062272
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 24 iterations.
+ Nonlinear residuals: 8.96319, 104.714, 3.50566
+
+ residual: 5.20165e-05
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 21 iterations.
+ Nonlinear residuals: 0.756761, 8.87878, 0.294157
+
+ residual: 4.36466e-06
+
+ Postprocessing:
+
+*** Timestep 10: t=95202.3 years
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving Stokes system... 30+4 iterations.
+ Nonlinear residuals: 292799, 2.74823e+06, 67420
+
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 13530.5, 148178, 5354.25
+
+ residual: 0.0794163
+ Solving temperature system... 9 iterations.
+ Solving composition system 1... 9 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 984.759, 11239.9, 390.392
+
+ residual: 0.00579044
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 28 iterations.
+ Nonlinear residuals: 75.5843, 878.741, 29.5862
+
+ residual: 0.000438834
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 6 iterations.
+ Solving Stokes system... 23 iterations.
+ Nonlinear residuals: 5.89602, 68.8871, 2.28763
+
+ residual: 3.3931e-05
+ Solving temperature system... 5 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 20 iterations.
+ Nonlinear residuals: 0.460968, 5.3964, 0.178174
+
+ residual: 2.64274e-06
+
+ Postprocessing:
+
+*** Timestep 11: t=100000 years
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving Stokes system... 30+3 iterations.
+ Nonlinear residuals: 180229, 1.78002e+06, 47625.8
+
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving Stokes system... 30+2 iterations.
+ Nonlinear residuals: 6885.89, 78122.5, 2763.21
+
+ residual: 0.0580192
+ Solving temperature system... 8 iterations.
+ Solving composition system 1... 8 iterations.
+ Solving Stokes system... 30 iterations.
+ Nonlinear residuals: 382.889, 4559.44, 150.945
+
+ residual: 0.0031694
+ Solving temperature system... 7 iterations.
+ Solving composition system 1... 7 iterations.
+ Solving Stokes system... 26 iterations.
+ Nonlinear residuals: 22.6769, 274.308, 8.80057
+
+ residual: 0.000184786
+ Solving temperature system... 6 iterations.
+ Solving composition system 1... 5 iterations.
+ Solving Stokes system... 21 iterations.
+ Nonlinear residuals: 1.37197, 16.6415, 0.528167
+
+ residual: 1.10899e-05
+ Solving temperature system... 4 iterations.
+ Solving composition system 1... 4 iterations.
+ Solving Stokes system... 17 iterations.
+ Nonlinear residuals: 0.0835044, 1.01291, 0.032067
+
+ residual: 6.73312e-07
+
+ Postprocessing:
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 10.1s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 84 | 3.23s | 32% |
+| Assemble composition system | 84 | 1.54s | 15% |
+| Assemble temperature system | 84 | 1.83s | 18% |
+| Build Stokes preconditioner | 1 | 0.0315s | 0.31% |
+| Build composition preconditioner| 84 | 0.254s | 2.5% |
+| Build temperature preconditioner| 12 | 0.0366s | 0.36% |
+| Solve Stokes system | 84 | 2.86s | 28% |
+| Solve composition system | 84 | 0.125s | 1.2% |
+| Solve temperature system | 84 | 0.127s | 1.3% |
+| Initialization | 2 | 0.0226s | 0.22% |
+| Postprocessing | 12 | 0.000754s | 0.0074% |
+| Setup dof systems | 1 | 0.0289s | 0.28% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/iterated_IMPES/statistics b/tests/iterated_IMPES/statistics
new file mode 100644
index 0000000..0cb16cf
--- /dev/null
+++ b/tests/iterated_IMPES/statistics
@@ -0,0 +1,96 @@
+# 1: Time step number
+# 2: Time (years)
+# 3: Number of mesh cells
+# 4: Number of Stokes degrees of freedom
+# 5: Number of temperature degrees of freedom
+# 6: Number of degrees of freedom for all compositions
+# 7: Iterations for temperature solver
+# 8: Iterations for composition solver 1
+# 9: Iterations for Stokes solver
+# 10: Velocity iterations in Stokes preconditioner
+# 11: Schur complement iterations in Stokes preconditioner
+# 12: Time step size (years)
+ 0 0.0000e+00 512 4851 2145 2145 0 0 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 2 3 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 0 0 0 0 0 0.0000e+00
+ 1 1.3251e+04 512 4851 2145 2145 0 0 0 0 0 1.3251e+04
+ 0 0.0000e+00 0 0 0 0 10 9 34 28 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 11 9 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 8 33 23 4 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 8 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 7 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 7 31 12 2 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 6 29 30 30 0.0000e+00
+ 2 2.5357e+04 512 4851 2145 2145 7 6 26 27 27 1.2106e+04
+ 0 0.0000e+00 0 0 0 0 11 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 10 33 23 4 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 9 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 29 30 30 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 6 27 28 28 0.0000e+00
+ 3 3.6192e+04 512 4851 2145 2145 5 5 24 25 25 1.0835e+04
+ 0 0.0000e+00 0 0 0 0 11 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 10 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 9 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 31 12 2 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 6 27 28 28 0.0000e+00
+ 4 4.6125e+04 512 4851 2145 2145 6 5 24 25 25 9.9330e+03
+ 0 0.0000e+00 0 0 0 0 11 12 34 28 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 11 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 9 32 17 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 8 31 12 2 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 7 27 28 28 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 5 5 23 24 24 0.0000e+00
+ 5 5.5855e+04 512 4851 2145 2145 4 4 20 21 21 9.7304e+03
+ 0 0.0000e+00 0 0 0 0 11 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 11 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 9 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 8 30 32 32 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 7 27 28 28 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 5 6 23 24 24 0.0000e+00
+ 6 6.5167e+04 512 4851 2145 2145 4 5 19 20 20 9.3117e+03
+ 0 0.0000e+00 0 0 0 0 11 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 11 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 10 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 30 32 32 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 7 26 27 27 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 5 6 22 23 23 0.0000e+00
+ 7 7.3772e+04 512 4851 2145 2145 4 5 19 20 20 8.6045e+03
+ 0 0.0000e+00 0 0 0 0 11 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 11 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 10 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 29 30 30 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 25 26 26 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 6 22 23 23 0.0000e+00
+ 8 8.1630e+04 512 4851 2145 2145 4 4 18 19 19 7.8586e+03
+ 0 0.0000e+00 0 0 0 0 12 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 11 11 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 9 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 29 30 30 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 25 26 26 0.0000e+00
+ 9 8.8740e+04 512 4851 2145 2145 6 6 21 22 22 7.1099e+03
+ 0 0.0000e+00 0 0 0 0 12 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 11 11 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 9 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 28 29 29 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 24 25 25 0.0000e+00
+10 9.5202e+04 512 4851 2145 2145 6 5 21 22 22 6.4620e+03
+ 0 0.0000e+00 0 0 0 0 12 12 34 29 5 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 11 10 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 9 9 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 28 29 29 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 6 23 24 24 0.0000e+00
+11 1.0000e+05 512 4851 2145 2145 5 5 20 21 21 4.7977e+03
+ 0 0.0000e+00 0 0 0 0 11 11 33 23 4 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 10 10 32 18 3 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 8 8 30 32 32 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 7 7 26 27 27 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 6 5 21 22 22 0.0000e+00
+ 0 0.0000e+00 0 0 0 0 4 4 17 18 18 5.4709e+03
More information about the CIG-COMMITS
mailing list