----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.6.0-pre (main, 0ca502825) -- . using deal.II 9.4.2 -- . with 32 bit indices and vectorization level 1 (128 bits) -- . using Trilinos 13.2.0 -- . using p4est 2.3.2 -- . running in DEBUG mode -- . running with 4 MPI processes ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&NewtonSolver=1&sha=0ca502825&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+22 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.00480522, 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.00158539, norm of the rhs: 3.38485e+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.00149458, norm of the rhs: 3.19097e+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.00115145, norm of the rhs: 2.45838e+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.00108018, norm of the rhs: 2.30621e+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.000888132, norm of the rhs: 1.89619e+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.000760367, norm of the rhs: 1.62341e+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.000627739, norm of the rhs: 1.34024e+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.000530886, norm of the rhs: 1.13346e+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 11: 0.000474957, norm of the rhs: 1.01405e+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.000409088, norm of the rhs: 8.73415e+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.000321352, norm of the rhs: 6.86095e+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.000242415, norm of the rhs: 5.17564e+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.000183743, norm of the rhs: 3.92296e+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.000152685, norm of the rhs: 3.25987e+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.00014238, norm of the rhs: 3.03986e+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.000118598, norm of the rhs: 2.53211e+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.000109499, norm of the rhs: 2.33784e+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.000106795, norm of the rhs: 2.28011e+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.92354e-05, norm of the rhs: 2.11871e+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+19 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.00434246, norm of the rhs: 9.27125e+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.00132294, norm of the rhs: 2.8245e+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.00132772, norm of the rhs: 2.8347e+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.00105014, norm of the rhs: 2.24206e+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.000944385, norm of the rhs: 2.01628e+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.000722221, norm of the rhs: 1.54196e+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.000603706, norm of the rhs: 1.28893e+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.000521305, norm of the rhs: 1.113e+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.000445941, norm of the rhs: 9.52095e+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.000395981, norm of the rhs: 8.45429e+11, 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 12: 0.000326272, norm of the rhs: 6.96599e+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 13: 0.000262243, norm of the rhs: 5.59894e+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.000210468, norm of the rhs: 4.49355e+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.000178781, norm of the rhs: 3.81703e+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.000155844, norm of the rhs: 3.32731e+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.000128508, norm of the rhs: 2.74368e+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.00011179, norm of the rhs: 2.38675e+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.89635e-05, norm of the rhs: 2.11289e+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, 2.304 kg/yr Total system matrix memory consumption: 39.63 MB. Total system matrix nnz: 2,925,282 system matrix nnz by block: 1,543,642 302,693 4 4 4 4 4 4 302,693 92 4 4 4 4 4 4 4 4 381,553 4 4 4 4 4 4 4 4 394,393 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 Total system preconditioner matrix memory consumption: 16.01 MB. Total system preconditioner matrix nnz: 828,027 system preconditioner matrix nnz by block: 772,870 4 4 4 4 4 4 4 4 54,933 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 Pressure min/avg/max: -2.75e+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... 8 iterations. -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.6.0-pre (main, 0ca502825) -- . using deal.II 9.4.2 -- . with 32 bit indices and vectorization level 1 (128 bits) -- . using Trilinos 13.2.0 -- . using p4est 2.3.2 -- . running in DEBUG mode -- . running with 4 MPI processes ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&NewtonSolver=1&sha=0ca502825&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+22 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.00480522, 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.00158539, norm of the rhs: 3.38485e+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.00149458, norm of the rhs: 3.19097e+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.00115145, norm of the rhs: 2.45838e+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.00108018, norm of the rhs: 2.30621e+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.000888132, norm of the rhs: 1.89619e+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.000760367, norm of the rhs: 1.62341e+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.000627739, norm of the rhs: 1.34024e+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.000530886, norm of the rhs: 1.13346e+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 11: 0.000474957, norm of the rhs: 1.01405e+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.000409088, norm of the rhs: 8.73415e+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.000321352, norm of the rhs: 6.86095e+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.000242415, norm of the rhs: 5.17564e+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.000183743, norm of the rhs: 3.92296e+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.000152685, norm of the rhs: 3.25987e+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.00014238, norm of the rhs: 3.03986e+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.000118598, norm of the rhs: 2.53211e+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.000109499, norm of the rhs: 2.33784e+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.000106795, norm of the rhs: 2.28011e+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.92354e-05, norm of the rhs: 2.11871e+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+19 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.00434246, norm of the rhs: 9.27125e+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.00132294, norm of the rhs: 2.8245e+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.00132772, norm of the rhs: 2.8347e+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.00105014, norm of the rhs: 2.24206e+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.000944385, norm of the rhs: 2.01628e+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.000722221, norm of the rhs: 1.54196e+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.000603706, norm of the rhs: 1.28893e+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.000521305, norm of the rhs: 1.113e+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.000445941, norm of the rhs: 9.52095e+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.000395981, norm of the rhs: 8.45429e+11, 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 12: 0.000326272, norm of the rhs: 6.96599e+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 13: 0.000262243, norm of the rhs: 5.59894e+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.000210468, norm of the rhs: 4.49355e+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.000178781, norm of the rhs: 3.81703e+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.000155844, norm of the rhs: 3.32731e+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.000128508, norm of the rhs: 2.74368e+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.00011179, norm of the rhs: 2.38675e+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.89635e-05, norm of the rhs: 2.11289e+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, 2.304 kg/yr Total system matrix memory consumption: 39.63 MB. Total system matrix nnz: 2,925,282 system matrix nnz by block: 1,543,642 302,693 4 4 4 4 4 4 302,693 92 4 4 4 4 4 4 4 4 381,553 4 4 4 4 4 4 4 4 394,393 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 Total system preconditioner matrix memory consumption: 16.01 MB. Total system preconditioner matrix nnz: 828,027 system preconditioner matrix nnz by block: 772,870 4 4 4 4 4 4 4 4 54,933 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 4 4 4 4 4 4 4 4 0 Pressure min/avg/max: -2.75e+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... 8 iterations. -------------------------------------------------------- An error occurred in line <2721> of file in function void dealii::DoFCellAccessor::get_dof_values(const InputVector&, ForwardIterator, ForwardIterator) const [with InputVector = dealii::TrilinosWrappers::MPI::Vector; ForwardIterator = double*; int dimension_ = 2; int space_dimension_ = 2; bool level_dof_access = false] The violated condition was: this->is_artificial() == false Additional information: Can't ask for DoF indices on artificial cells. Stacktrace: ----------- #0 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, double*, double*) const #1 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, dealii::Vector&) const #2 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::get_vertices(dealii::TriaIterator > const&) const #3 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::compute_mapping_support_points(dealii::TriaIterator > const&) const #4 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ<2, 2>::get_bounding_box(dealii::TriaIterator > const&) const #5 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::GridTools::Cache<2, 2>::get_cell_bounding_boxes_rtree() const #6 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: std::tuple::active_cell_iterator, std::allocator::active_cell_iterator> >, std::vector, std::allocator > >, std::allocator, std::allocator > > > >, std::vector >, std::allocator > > >, std::vector > > dealii::GridTools::compute_point_locations_try_all<2, 2>(dealii::GridTools::Cache<2, 2> const&, std::vector, std::allocator > > const&, dealii::Triangulation<2, 2>::active_cell_iterator const&) #7 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::compute_point_locations(std::vector, std::allocator > > const&, std::vector >, std::allocator > > >&, std::vector, std::allocator > >, std::allocator, std::allocator > > > >&, std::vector >, std::allocator > > >&) const #8 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::vector_value_list(std::vector, std::allocator > > const&, std::vector, std::allocator > >&) const #9 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::value_list(std::vector, std::allocator > > const&, std::vector >&, unsigned int) -------------------------------------------------------- An error occurred in line <2721> of file in function void dealii::DoFCellAccessor::get_dof_values(const InputVector&, ForwardIterator, ForwardIterator) const [with InputVector = dealii::TrilinosWrappers::MPI::Vector; ForwardIterator = double*; int dimension_ = 2; int space_dimension_ = 2; bool level_dof_access = false] The violated condition was: this->is_artificial() == false Additional information: Can't ask for DoF indices on artificial cells. Stacktrace: ----------- #0 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, double*, double*) const #1 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, dealii::Vector&) const #2 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::get_vertices(dealii::TriaIterator > const&) const #3 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::compute_mapping_support_points(dealii::TriaIterator > const&) const #4 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ<2, 2>::get_bounding_box(dealii::TriaIterator > const&) const #5 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::GridTools::Cache<2, 2>::get_cell_bounding_boxes_rtree() const #6 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: std::tuple::active_cell_iterator, std::allocator::active_cell_iterator> >, std::vector, std::allocator > >, std::allocator, std::allocator > > > >, std::vector >, std::allocator > > >, std::vector > > dealii::GridTools::compute_point_locations_try_all<2, 2>(dealii::GridTools::Cache<2, 2> const&, std::vector, std::allocator > > const&, dealii::Triangulation<2, 2>::active_cell_iterator const&) #7 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::compute_point_locations(std::vector, std::allocator > > const&, std::vector >, std::allocator > > >&, std::vector, std::allocator > >, std::allocator, std::allocator > > > >&, std::vector >, std::allocator > > >&) const #8 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::vector_value_list(std::vector, std::allocator > > const&, std::vector, std::allocator > >&) const #9 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::value_list(std::vector, std::allocator > > const&, std::vector >&, unsigned int) -------------------------------------------------------- An error occurred in line <2721> of file in function void dealii::DoFCellAccessor::get_dof_values(const InputVector&, ForwardIterator, ForwardIterator) const [with InputVector = dealii::TrilinosWrappers::MPI::Vector; ForwardIterator = double*; int dimension_ = 2; int space_dimension_ = 2; bool level_dof_access = false] The violated condition was: this->is_artificial() == false Additional information: Can't ask for DoF indices on artificial cells. Stacktrace: ----------- #0 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, double*, double*) const #1 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, dealii::Vector&) const #2 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::get_vertices(dealii::TriaIterator > const&) const #3 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::compute_mapping_support_points(dealii::TriaIterator > const&) const #4 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ<2, 2>::get_bounding_box(dealii::TriaIterator > const&) const #5 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::GridTools::Cache<2, 2>::get_cell_bounding_boxes_rtree() const #6 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: std::tuple::active_cell_iterator, std::allocator::active_cell_iterator> >, std::vector, std::allocator > >, std::allocator, std::allocator > > > >, std::vector >, std::allocator > > >, std::vector > > dealii::GridTools::compute_point_locations_try_all<2, 2>(dealii::GridTools::Cache<2, 2> const&, std::vector, std::allocator > > const&, dealii::Triangulation<2, 2>::active_cell_iterator const&) #7 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::compute_point_locations(std::vector, std::allocator > > const&, std::vector >, std::allocator > > >&, std::vector, std::allocator > >, std::allocator, std::allocator > > > >&, std::vector >, std::allocator > > >&) const #8 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::vector_value_list(std::vector, std::allocator > > const&, std::vector, std::allocator > >&) const #9 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::value_list(std::vector, std::allocator > > const&, std::vector >&, unsigned int) -------------------------------------------------------- An error occurred in line <2721> of file in function void dealii::DoFCellAccessor::get_dof_values(const InputVector&, ForwardIterator, ForwardIterator) const [with InputVector = dealii::TrilinosWrappers::MPI::Vector; ForwardIterator = double*; int dimension_ = 2; int space_dimension_ = 2; bool level_dof_access = false] The violated condition was: this->is_artificial() == false Additional information: Can't ask for DoF indices on artificial cells. Stacktrace: ----------- #0 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, double*, double*) const #1 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: void dealii::DoFCellAccessor<2, 2, false>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, dealii::Vector&) const #2 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::get_vertices(dealii::TriaIterator > const&) const #3 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ1Eulerian<2, dealii::TrilinosWrappers::MPI::Vector, 2>::compute_mapping_support_points(dealii::TriaIterator > const&) const #4 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::MappingQ<2, 2>::get_bounding_box(dealii::TriaIterator > const&) const #5 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::GridTools::Cache<2, 2>::get_cell_bounding_boxes_rtree() const #6 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: std::tuple::active_cell_iterator, std::allocator::active_cell_iterator> >, std::vector, std::allocator > >, std::allocator, std::allocator > > > >, std::vector >, std::allocator > > >, std::vector > > dealii::GridTools::compute_point_locations_try_all<2, 2>(dealii::GridTools::Cache<2, 2> const&, std::vector, std::allocator > > const&, dealii::Triangulation<2, 2>::active_cell_iterator const&) #7 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::compute_point_locations(std::vector, std::allocator > > const&, std::vector >, std::allocator > > >&, std::vector, std::allocator > >, std::allocator, std::allocator > > > >&, std::vector >, std::allocator > > >&) const #8 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::vector_value_list(std::vector, std::allocator > > const&, std::vector, std::allocator > >&) const #9 /root/bin/deal.II-v9.4.2/lib/libdeal_II.g.so.9.4.2: dealii::Functions::FEFieldFunction<2, dealii::TrilinosWrappers::MPI::BlockVector, 2>::value_list(std::vector, std::allocator > > const&, std::vector >&, unsigned int) const #10 /root/aspect/build/aspect: aspect::MaterialModel::Rheology::StrainDependent<2>::fill_reaction_outputs(aspect::MaterialModel::MaterialModelInputs<2> const&, int, double, bool, aspect::MaterialModel::MaterialModelOutputs<2>&) const #11 /root/aspect/build/aspect: aspect::MaterialModel::ViscoPlastic<2>::evaluate(aspect::MaterialModel::MaterialModelInputs<2> const&, aspect::MaterialModel::MaterialModelOutputs<2>&) const #12 /root/aspect/build/aspect: void aspect::Simulator<2>::get_artificial_viscosity(dealii::Vector&, aspect::Simulator<2>::AdvectionField const&, bool) const #13 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_advection_system(aspect::Simulator<2>::AdvectionField const&) #14 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_and_solve_temperature(double const&, double*) #15 /root/aspect/build/aspect: aspect::Simulator<2>::solve_single_advection_iterated_defect_correction_stokes() #16 /root/aspect/build/aspect: aspect::Simulator<2>::solve_timestep() #17 /root/aspect/build/aspect: aspect::Simulator<2>::run() #18 /root/aspect/build/aspect: void run_simulator<2>(std::__cxx11::basic_string, std::allocator > const&, std::__cxx11::basic_string, std::allocator > const&, bool, bool, bool) #19 /root/aspect/build/aspect: main -------------------------------------------------------- Calling MPI_Abort now. To break execution in a GDB session, execute 'break MPI_Abort' before running. You can also put the following into your ~/.gdbinit: set breakpoint pending on break MPI_Abort set breakpoint pending auto const #10 /root/aspect/build/aspect: aspect::MaterialModel::Rheology::StrainDependent<2>::fill_reaction_outputs(aspect::MaterialModel::MaterialModelInputs<2> const&, int, double, bool, aspect::MaterialModel::MaterialModelOutputs<2>&) const #11 /root/aspect/build/aspect: aspect::MaterialModel::ViscoPlastic<2>::evaluate(aspect::MaterialModel::MaterialModelInputs<2> const&, aspect::MaterialModel::MaterialModelOutputs<2>&) const #12 /root/aspect/build/aspect: void aspect::Simulator<2>::get_artificial_viscosity(dealii::Vector&, aspect::Simulator<2>::AdvectionField const&, bool) const #13 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_advection_system(aspect::Simulator<2>::AdvectionField const&) #14 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_and_solve_temperature(double const&, double*) #15 /root/aspect/build/aspect: aspect::Simulator<2>::solve_single_advection_iterated_defect_correction_stokes() #16 /root/aspect/build/aspect: aspect::Simulator<2>::solve_timestep() #17 /root/aspect/build/aspect: aspect::Simulator<2>::run() #18 /root/aspect/build/aspect: void run_simulator<2>(std::__cxx11::basic_string, std::allocator > const&, std::__cxx11::basic_string, std::allocator > const&, bool, bool, bool) #19 /root/aspect/build/aspect: main -------------------------------------------------------- Calling MPI_Abort now. To break execution in a GDB session, execute 'break MPI_Abort' before running. You can also put the following into your ~/.gdbinit: set breakpoint pending on break MPI_Abort set breakpoint pending auto const #10 /root/aspect/build/aspect: aspect::MaterialModel::Rheology::StrainDependent<2>::fill_reaction_outputs(aspect::MaterialModel::MaterialModelInputs<2> const&, int, double, bool, aspect::MaterialModel::MaterialModelOutputs<2>&) const #11 /root/aspect/build/aspect: aspect::MaterialModel::ViscoPlastic<2>::evaluate(aspect::MaterialModel::MaterialModelInputs<2> const&, aspect::MaterialModel::MaterialModelOutputs<2>&) const #12 /root/aspect/build/aspect: void aspect::Simulator<2>::get_artificial_viscosity(dealii::Vector&, aspect::Simulator<2>::AdvectionField const&, bool) const #13 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_advection_system(aspect::Simulator<2>::AdvectionField const&) #14 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_and_solve_temperature(double const&, double*) #15 /root/aspect/build/aspect: aspect::Simulator<2>::solve_single_advection_iterated_defect_correction_stokes() #16 /root/aspect/build/aspect: aspect::Simulator<2>::solve_timestep() #17 /root/aspect/build/aspect: aspect::Simulator<2>::run() #18 /root/aspect/build/aspect: void run_simulator<2>(std::__cxx11::basic_string, std::allocator > const&, std::__cxx11::basic_string, std::allocator > const&, bool, bool, bool) #19 /root/aspect/build/aspect: main -------------------------------------------------------- Calling MPI_Abort now. To break execution in a GDB session, execute 'break MPI_Abort' before running. You can also put the following into your ~/.gdbinit: set breakpoint pending on break MPI_Abort set breakpoint pending auto const #10 /root/aspect/build/aspect: aspect::MaterialModel::Rheology::StrainDependent<2>::fill_reaction_outputs(aspect::MaterialModel::MaterialModelInputs<2> const&, int, double, bool, aspect::MaterialModel::MaterialModelOutputs<2>&) const #11 /root/aspect/build/aspect: aspect::MaterialModel::ViscoPlastic<2>::evaluate(aspect::MaterialModel::MaterialModelInputs<2> const&, aspect::MaterialModel::MaterialModelOutputs<2>&) const #12 /root/aspect/build/aspect: void aspect::Simulator<2>::get_artificial_viscosity(dealii::Vector&, aspect::Simulator<2>::AdvectionField const&, bool) const #13 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_advection_system(aspect::Simulator<2>::AdvectionField const&) #14 /root/aspect/build/aspect: aspect::Simulator<2>::assemble_and_solve_temperature(double const&, double*) #15 /root/aspect/build/aspect: aspect::Simulator<2>::solve_single_advection_iterated_defect_correction_stokes() #16 /root/aspect/build/aspect: aspect::Simulator<2>::solve_timestep() #17 /root/aspect/build/aspect: aspect::Simulator<2>::run() #18 /root/aspect/build/aspect: void run_simulator<2>(std::__cxx11::basic_string, std::allocator > const&, std::__cxx11::basic_string, std::allocator > const&, bool, bool, bool) #19 /root/aspect/build/aspect: main -------------------------------------------------------- Calling MPI_Abort now. To break execution in a GDB session, execute 'break MPI_Abort' before running. You can also put the following into your ~/.gdbinit: set breakpoint pending on break MPI_Abort set breakpoint pending auto -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 3 in communicator MPI_COMM_WORLD with errorcode 255. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them.