----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.6.0-pre (master, 69faf2563) -- . using deal.II 9.4.2 -- . with 32 bit indices and vectorization level 3 (512 bits) -- . using Trilinos 13.2.0 -- . using p4est 2.3.2 -- . running in DEBUG mode -- . running with 96 MPI processes ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- The output directory provided in the input file appears not to exist. ASPECT will create it for you. ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&NewtonSolver=1&sha=69faf2563&src=code ----------------------------------------------------------------------------- Number of active cells: 3,200 (on 4 levels) Number of degrees of freedom: 107,649 (26,082+3,321+13,041+13,041+13,041+13,041+13,041+13,041) Number of mesh deformation degrees of freedom: 6642 Solving mesh displacement system... 0 iterations. *** Timestep 0: t=0 years, dt=0 years Solving mesh displacement system... 0 iterations. Solving temperature system... 0 iterations. Skipping noninitial_plastic_strain composition solve because RHS is zero. Solving plastic_strain system ... 0 iterations. Solving crust_upper system ... 0 iterations. Solving crust_lower system ... 0 iterations. Solving mantle_lithosphere system ... 0 iterations. Initial Newton Stokes residual = 2.13503e+15, v = 2.13385e+15, p = 7.11214e+13 Rebuilding Stokes preconditioner... Solving Stokes system... 0+21 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 1, norm of the rhs: 2.13503e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+24 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.00480521, norm of the rhs: 1.02593e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.00158891, norm of the rhs: 3.39236e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.00149651, norm of the rhs: 3.1951e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 0.00115166, norm of the rhs: 2.45883e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 0.00108, norm of the rhs: 2.30583e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 0.000886965, norm of the rhs: 1.8937e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 0.000760206, norm of the rhs: 1.62306e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 9: 0.000627395, norm of the rhs: 1.33951e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 10: 0.000530564, norm of the rhs: 1.13277e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 11: 0.000473717, norm of the rhs: 1.0114e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 12: 0.000408903, norm of the rhs: 8.7302e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 13: 0.000321457, norm of the rhs: 6.8632e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 14: 0.000242215, norm of the rhs: 5.17136e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 15: 0.000183669, norm of the rhs: 3.92139e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 16: 0.000152634, norm of the rhs: 3.25878e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 17: 0.000142276, norm of the rhs: 3.03764e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 18: 0.000118645, norm of the rhs: 2.53311e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 19: 0.000109462, norm of the rhs: 2.33705e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 20: 0.000106817, norm of the rhs: 2.28057e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 21: 9.93864e-05, norm of the rhs: 2.12193e+11, newton_derivative_scaling_factor: 0 Number of active cells: 6,080 (on 5 levels) Number of degrees of freedom: 204,549 (49,570+6,269+24,785+24,785+24,785+24,785+24,785+24,785) Number of mesh deformation degrees of freedom: 12538 Solving mesh displacement system... 0 iterations. *** Timestep 0: t=0 years, dt=0 years Solving mesh displacement system... 0 iterations. Solving temperature system... 0 iterations. Skipping noninitial_plastic_strain composition solve because RHS is zero. Solving plastic_strain system ... 0 iterations. Solving crust_upper system ... 0 iterations. Solving crust_lower system ... 0 iterations. Solving mantle_lithosphere system ... 0 iterations. Initial Newton Stokes residual = 2.13502e+15, v = 2.13385e+15, p = 7.0922e+13 Rebuilding Stokes preconditioner... Solving Stokes system... 0+21 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 1, norm of the rhs: 2.13502e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.00434243, norm of the rhs: 9.27118e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.00131223, norm of the rhs: 2.80163e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.00132611, norm of the rhs: 2.83127e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 0.001047, norm of the rhs: 2.23537e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 0.000938719, norm of the rhs: 2.00419e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 0.000722251, norm of the rhs: 1.54202e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 0.000606582, norm of the rhs: 1.29507e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 9: 0.000521357, norm of the rhs: 1.11311e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 10: 0.000446079, norm of the rhs: 9.52388e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 11: 0.000396383, norm of the rhs: 8.46287e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 12: 0.000327231, norm of the rhs: 6.98646e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 13: 0.000262042, norm of the rhs: 5.59465e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 14: 0.000210182, norm of the rhs: 4.48744e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 15: 0.000178976, norm of the rhs: 3.82119e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 16: 0.000155825, norm of the rhs: 3.32691e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 17: 0.000128441, norm of the rhs: 2.74225e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 18: 0.000111806, norm of the rhs: 2.38708e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 19: 9.89962e-05, norm of the rhs: 2.11359e+11, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: 0/0/0 // 0/1.497/7.424e+09 // -0.125/1/4.058e+09 // -0.125/1.125/4e+09 // 0/1.125/1.194e+10 Heat flux densities for boundaries: 0.1898 W/m, 0.1898 W/m, -0.347 W/m, 0.05506 W/m Heat fluxes through boundary parts: 1.898e+04 W, 1.898e+04 W, -6.94e+04 W, 1.101e+04 W Mass fluxes through boundary parts: 7.626e+05 kg/yr, 7.626e+05 kg/yr, -1.606e+06 kg/yr, 1.531 kg/yr Total system matrix memory consumption: 46.35 MB. Total system matrix nnz: 2,930,311 system matrix nnz by block: 1,543,642 302,693 96 96 96 96 96 96 302,693 183 96 96 96 96 96 96 96 96 381,553 95 95 95 95 95 96 96 95 394,393 95 95 95 95 96 96 95 95 0 95 95 95 96 96 95 95 95 0 95 95 96 96 95 95 95 95 0 95 96 96 95 95 95 95 95 0 Total system preconditioner matrix memory consumption: 22.73 MB. Total system preconditioner matrix nnz: 833,149 system preconditioner matrix nnz by block: 772,870 96 96 96 96 96 96 96 96 54,933 96 96 96 96 96 96 96 96 0 95 95 95 95 95 96 96 95 0 95 95 95 95 96 96 95 95 0 95 95 95 96 96 95 95 95 0 95 95 96 96 95 95 95 95 0 95 96 96 95 95 95 95 95 0 Pressure min/avg/max: -2.748e+07 Pa, 1.393e+09 Pa, 2.996e+09 Pa Temperature min/avg/max: 273 K, 998.3 K, 1613 K Topography min/max: 0 m, 0 m RMS, max velocity: 0.00236 m/year, 0.00348 m/year Writing graphical output: output-continental_extension/solution/solution-00000 *** Timestep 1: t=20000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 10 iterations. ============ Job Information =================================================== Submitted: 2023-08-02T16:58:20 Started: 2023-08-02T16:58:20 Ended: 2023-08-02T16:59:21 Elapsed: 2 min, Limit: 60 min, Difference: 58 min (HT-)CPUs: 192, Nodes: 1 Estimated core-h: 1.6226 ================================================================================