Merge pull request #6003 from MFraters/switch_continental_extension_t…

…o_amg Switch continental extension to AMG solver.
geodynamics · Aug 20, 2024 · 5148ea0 · 5148ea0
2 parents a15f78a + 90b0c2e
commit 5148ea0
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 8 deletions.
diff --git a/cookbooks/continental_extension/continental_extension.prm b/cookbooks/continental_extension/continental_extension.prm
@@ -72,6 +72,7 @@ end
 subsection Solver parameters
   subsection Stokes solver parameters
     set Number of cheap Stokes solver steps = 0
+    set Stokes solver type = block AMG
   end
 
   subsection Newton solver parameters

diff --git a/tests/continental_extension/screen-output b/tests/continental_extension/screen-output
@@ -12,10 +12,11 @@ Number of mesh deformation degrees of freedom: 462
    Solving upper_crust system ... 0 iterations.
    Solving lower_crust system ... 0 iterations.
    Solving lithospheric_mantle system ... 0 iterations.
-   Initial Newton Stokes residual = 2.88429e+14, v = 2.79119e+14, p = 7.26906e+13
+   Initial Newton Stokes residual = 3.04034e+14, v = 2.95216e+14, p = 7.26906e+13
 
-   Solving Stokes system... 0+56 iterations.
-      Newton system information: Norm of the rhs: 2.88429e+14
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 0+20 iterations.
+      Newton system information: Norm of the rhs: 3.04034e+14
       Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1
 
 
@@ -27,18 +28,19 @@ Continuing to the next timestep even though solution is not fully converged.
 
 *** Timestep 1:  t=10000 years, dt=10000 years
    Solving mesh surface diffusion
-   Solving mesh displacement system... 4 iterations.
+   Solving mesh displacement system... 1 iterations.
    Solving temperature system... 7 iterations.
    Solving noninitial_plastic_strain system ... 12 iterations.
    Solving plastic_strain system ... 9 iterations.
    Solving upper_crust system ... 9 iterations.
    Solving lower_crust system ... 9 iterations.
    Solving lithospheric_mantle system ... 9 iterations.
-   Initial Newton Stokes residual = 7.74397e+16, v = 7.74397e+16, p = 1.21967e+13
+   Initial Newton Stokes residual = 8.15787e+16, v = 8.15787e+16, p = 1.24975e+13
 
-   Solving Stokes system... 0+21 iterations.
-      Newton system information: Norm of the rhs: 2.35285e+13
-      Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 0.00030383
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 0+18 iterations.
+      Newton system information: Norm of the rhs: 2.11532e+13
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 0.000259297
 
 
 WARNING: The nonlinear solver in the current timestep failed to converge.