diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/README.md b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/README.md index 7acc5876a84..09abf4470c3 100644 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/README.md +++ b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/README.md @@ -5,11 +5,10 @@ can be used to recreate the nonlinear channel flow figures from {cite:t}`fraters:etal:2019`. Start by compiling the `simple_nonlinear` plugin in this directory as described (at the example of another benchmark) in {ref}`sec:benchmark-run` . -After compiling the plugin adjust `metabash.sh` and/or `bash.sh` to reflect -the parameter search you are interested in and run the script. -Running the `bash.sh` script without changes will recreate the model runs of Fig. 1 and Fig. 2 in {cite:t}`fraters:etal:2019` +After compiling the plugin you can run (or adjust) the script `run.sh`, which will run all the models necessary +to recreate Fig. 1 and Fig. 2 of {cite:t}`fraters:etal:2019` and the other figures on this page (by default the script will assume your processor has 4 compute cores, but you can adjust the number in the script). -After running the model series you can use the gnuplot files in `plotfiles/` to recreate the following figures. +After running the model series you can use the gnuplot files `plot_v.gnuplot` and `plot_t.gnuplot` to recreate the following figures. ```{figure-md} fig:benchmark-newton-nonlinear-channel-flow-tractions @@ -24,3 +23,31 @@ Convergence history for several methods for a rheology with n = 3 where in- and ``` The nonlinear convergence behavior is shown in {numref}`fig:benchmark-newton-nonlinear-channel-flow-tractions` and {numref}`fig:benchmark-newton-nonlinear-channel-flow-velocities`. + +Note, that the linear residual is stalling in some cases. This is related to the choice of boundary velocity (in {numref}`fig:benchmark-newton-nonlinear-channel-flow-velocities`) +and boundary traction (in {numref}`fig:benchmark-newton-nonlinear-channel-flow-tractions`). You can adjust the base input files `input_t.prm` and `input_v.prm` to +test what happens when varying the boundary tractions or velocities. + +Also note that compared to the original paper we have added a new model series to the figure that uses the default Picard solver `single Advection, iterated Stokes` (purple lines). You can see that this solver converges at about the same rate as the defect correction Picard solver, however its residual eventually stagnates at a fixed value. This is the expected behavior of a solver that computes the actual solution vector instead of an update to the solution vector. The accuracy of the linear solver in these cases limits the accuracy that the nonlinear solver can reach and changing the linear solver tolerance controls the exact residual at which the nonlinear solver stagnates. At the same time this illustrates the advantage of defect-correction solvers (using either Picard or Newton iterations). Even with a coarse linear solver tolerance these nonlinear solvers continue to converge to the actual solution. + +Since the publication of {cite:t}`fraters:etal:2019` we have continued to work on improving and optimizing the Newton solver implementation and parameters. The result of this work are default parameters that can solve many models faster than using the parameters originally used. In particular two choices influence the solver convergence and are discussed here: + +- We can make use of a coarser maximum linear solver tolerance than the one used for the publications (here `1e-20`). See the parameter [](parameters:Solver_20parameters/Newton_20solver_20parameters/Maximum_20linear_20Stokes_20solver_20tolerance). For the example with velocity boundary conditions this change cuts the cost for the linear Stokes solver in half, while not affecting the number of nonlinear iterations significantly. + +- Choosing less (but not 0) line search iterations actually benefits the stability of the Newton solver. This is because too many line search iterations will prevent the Newton solver from modifying the solution to a value that has a higher residual. In other words if the solver gets stuck in a local minimum there is no way for it to leave this minimum behind. Therefore by default we now limit the number of line-search iterations. See the parameter [](parameters:Solver_20parameters/Newton_20solver_20parameters/Max_20Newton_20line_20search_20iterations) for more details. + +Applying these modifications here is a version of the figures above that was run with the current (as of version 3.0) default parameters of ASPECT: + +```{figure-md} fig:benchmark-newton-nonlinear-channel-flow-tractions-default-parameters + + +As above for the model with pressure boundary conditions, but using the updated default solver parameters. +``` + +```{figure-md} fig:benchmark-newton-nonlinear-channel-flow-velocities-default-parameters + + +As above for the model with velocity boundary conditions, but using the updated default solver parameters. +``` + +You can easily recreate these figures by replacing `mLT_1e-8` with `mLT_1e-2`, and by replacing `LS_100` with `LS_5` in the plotting scripts. diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/auto_logs/.gitignore b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/auto_logs/.gitignore deleted file mode 100644 index 5e7d2734cfc..00000000000 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/auto_logs/.gitignore +++ /dev/null @@ -1,4 +0,0 @@ -# Ignore everything in this directory -* -# Except this file -!.gitignore diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash.sh b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash.sh deleted file mode 100755 index ce03fd92dc0..00000000000 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash.sh +++ /dev/null @@ -1,171 +0,0 @@ -#!/bin/bash - -version="b1" -build_dir="../../../build/" -declare -a boundary_type=("v" "t") # velocity or traction boundary conditions -declare -a grid=(4) # global refinement levels -declare -a agrid=(0) # adaptive refinement levels -declare -a UFS=("false") # use Newton failsafe -declare -a NSP=("SPD") # ("Symmetric" "PD" "none") # Use Newton stabilization preconditioner -declare -a NSA=("SPD" "none") # ("Symmetric" "PD" "none") # Use Newton stabilization A block -declare -a n=(3) #1 2 5 10 25 50 100) -I=150 # maximum nonlinear iterations -declare -a P=(5 10 15 150) # (0 5 10 15 20 25 30 50 150) # Picard iterations -declare -a LS=(0 100) # (0 5 10 100) # Line search iterations -declare -a ST=(1e-20) # Nonlinear newton solver switch tolerance -declare -a LT=("1e-5") # "1e-8" # Linear solver tolerance -declare -a NLT=("1e-14") # "1e-9" "1e-10" "1e-14" "1e-20" # Nonlinear solver tolerance -declare -a ABT=("1e-2") # A block solver tolerance -declare -a RSM=("false" "true") # "true" #"false" # Residual scaling method -declare -a SF=("9e-1") # SPD safety factor -declare -a UDS=("false") # ("false" "true") # Use deviator of strain-rate -declare -a AEW=("false") # always use Eisenstat Walker, even for picard -AV=-1 # Viscosity averaging -COMP="0" # "4e-12" # compressibility -declare -a MLT=("9e-1" "1e-8") # ("9e-1" "1e-1" "1e-2" "1e-8") # maximum linear Stokes solver tolerance -materialmodelnameShort="SNL" -materialmodelname="SNL" -phi=0 -declare -a vel=(0) #(25 50 125) # 50 125) #25 50 125) -declare -a BV=(1e19) #("1e23" "1e24" "5e24") #"1e23" "1e24" "5e24") -processes=4 - -CohesionLine="s/set List of cohesions.*/ set List of cohesions = 1e8,0/g" -PhiLine="s/set List of angles of internal friction.*/ set List of angles of internal friction = $phi,0/g" - -if [ $materialmodelnameShort == "DP" ]; then - materialmodelname="drucker prager compositions" -elif [ $materialmodelnameShort == "vM" ]; then - materialmodelname="drucker prager compositions" -elif [ $materialmodelnameShort == "SNL" ]; then - materialmodelname="simple nonlinear compositions" - CohesionLine="s/ set List of cohesion of fields.*/ /g" - PhiLine="s/set List of angles of internal friction of fields.*/ /g" -elif [ $materialmodelnameShort == "VP2" ]; then - materialmodelname="viscoplastic2" -fi - -SOLVER_SHORT="NS" #"itAdandSt" #"NS" #NS" -SOLVER="NS" -if [ $SOLVER_SHORT == "NS" ]; then - SOLVER="iterated Advection and Newton Stokes" -elif [ $SOLVER_SHORT == "itAdandSt" ]; then - SOLVER="iterated Advection and Stokes" -fi - - -for i_grid in "${grid[@]}" -do - for i_agrid in "${agrid[@]}" - do - for i_NSP in "${NSP[@]}" - do - for i_NSA in "${NSA[@]}" - do - for i_NLT in "${NLT[@]}" - do - for i_ABT in "${ABT[@]}" - do - for i_ST in "${ST[@]}" - do - for i_UDS in "${UDS[@]}" - do - for i_SF in "${SF[@]}" - do - for i_LT in "${LT[@]}" - do - for i_n in "${n[@]}" - do - for i_UFS in "${UFS[@]}" - do - for i_P in "${P[@]}" - do - for i_boundary_type in "${boundary_type[@]}" - do - for i_LS in "${LS[@]}" - do - for i_AEW in "${AEW[@]}" - do - for i_MLT in "${MLT[@]}" - do - for i_RSM in "${RSM[@]}" - do - for i_BV in "${BV[@]}" - do - for i_vel in "${vel[@]}" - do - U="7.92219116e-11" - if [ $i_vel == 25 ]; then - U="7.92219116e-11" - elif [ $i_vel == 50 ]; then - U="1.58443823e-10" - elif [ $i_vel == 125 ]; then - U="3.96109558e-10" - fi - - dirname_clean="$version""$materialmodelnameShort""_BT""$i_boundary_type""_""$SOLVER_SHORT""_ST""$i_ST""_UFS""$i_UFS""_NSP-""$i_NSP""_NSA-""$i_NSA""_C""$COMP""_g""$i_grid""_ag""$i_agrid""_AEW""$i_AEW""_UDS""$i_UDS""_SF""$i_SF""_NLT""$i_NLT""_ABT""$i_ABT""_LT""$i_LT""_mLT""$i_MLT""_I""$I""_P""$i_P""_EW1""_theta1""_LS""$i_LS""_RSM""$i_RSM""_AV""$AV""_phi""$phi""_vel""$i_vel""_BV""$i_BV""_n""$i_n" - dirname="results/""$dirname_clean" - infilename="$dirname""/input.prm" - outfilename="$dirname""/output.log" - errorfilename="$dirname""/error.log" - outplotfilename="$dirname""/plot.dat" - - input_name="input_""$i_boundary_type"".prm" - - mkdir -p $dirname - - echo "$dirname" - sed \ - -e "$PhiLine" \ - -e "$CohesionLine" \ - -e "s/ set Function expression = if(x<60e3,.*/ set Function expression = if(x<60e3,$U,-$U);0/g" \ - -e "s/ set Nonlinear Newton solver switch tolerance.*/ set Nonlinear Newton solver switch tolerance = $i_ST/g" \ - -e "s/set Output directory .*/set Output directory = results\/$dirname_clean/g" \ - -e "s/ set Model name = .*/ set Model name = $materialmodelname/g" \ - -e "s/set Nonlinear solver scheme.*/set Nonlinear solver scheme = $SOLVER/g" \ - -e "s/ set Initial global refinement = .*/ set Initial global refinement = $i_grid/g" \ - -e "s/ set Initial adaptive refinement = .*/ set Initial adaptive refinement = $i_agrid/g" \ - -e "s/ set List of stress exponents of fields = .*/ set List of stress exponents of fields = $i_n, 1/g" \ - -e "s/set Stress exponent .*/ set Stress exponent = $i_n/g" \ - -e "s/set Function constants = n .*/ set Function constants = n = $i_n/g" \ - -e "s/ set Viscosity averaging p = .*/ set Viscosity averaging p = $AV/g" \ - -e "s/set Max nonlinear iterations = .*/set Max nonlinear iterations = $I /g" \ - -e "s/set Linear solver tolerance =.*/set Linear solver tolerance = $i_LT/g" \ - -e "s/set Nonlinear solver tolerance.*/set Nonlinear solver tolerance = $i_NLT/g" \ - -e "s/set Linear solver A block tolerance.*/set Linear solver A block tolerance = $i_ABT/g" \ - -e "s/ set Reference compressibility .*/ set Reference compressibility = $COMP/g" \ - -e "s/set Max pre-Newton nonlinear iterations = .*/ set Max pre-Newton nonlinear iterations = $i_P/g" \ - -e "s/set Use Newton failsafe = .*/set Use Newton failsafe = $i_UFS/g" \ - -e "s/set Stabilization preconditioner = .*/set Stabilization preconditioner = $i_NSP/g" \ - -e "s/set Stabilization velocity block = .*/set Stabilization velocity block = $i_NSA/g" \ - -e "s/set SPD safety factor = .*/set SPD safety factor = $i_SF/g" \ - -e "s/set Use deviator of strain-rate = .*/set Use deviator of strain-rate = $i_UDS/g" \ - -e "s/set Max Newton line search iterations = .*/ set Max Newton line search iterations = $i_LS/g" \ - -e "s/set Maximum linear Stokes solver tolerance =.*/set Maximum linear Stokes solver tolerance = $i_MLT/g" \ - -e "s/set Use Newton residual scaling method .*/ set Use Newton residual scaling method = $i_RSM/g" \ - $input_name > "$infilename" - - nohup mpirun -np $processes $build_dir./aspect $infilename > $outfilename 2>$errorfilename - - grep "Relative nonlinear residual " $outfilename > $outplotfilename - - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done -done diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash_autos/.gitignore b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash_autos/.gitignore deleted file mode 100644 index 5e7d2734cfc..00000000000 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash_autos/.gitignore +++ /dev/null @@ -1,4 +0,0 @@ -# Ignore everything in this directory -* -# Except this file -!.gitignore diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash_itIMPES.sh b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash_itIMPES.sh deleted file mode 100755 index 7092d31956b..00000000000 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/bash_itIMPES.sh +++ /dev/null @@ -1,168 +0,0 @@ -!/bin/bash -version="a5" -declare -a boundary_type=("v") # "t") -declare -a grid=(5) -declare -a UFS=("false") # "true") # use failsafe -declare -a NSP=("SPD") # "none") # "Symmetric" "PD" "none") #Use Newton stabilization preconditioner -declare -a NSA=("SPD") # "symmetric" "none") # "Symmetric" "PD" "none") #Use Newton stabilization A block -declare -a agrid=(0) -declare -a n=(3) #1 2 5 10 25 50 100) -I=150 -declare -a P=(150) #0 5 10 50 15 150 20 25 30) #0 5 10 15 20 150) # 1 2 3 4 5 10 15 25 50 100) -declare -a LS=(0) # 5 10) -declare -a ST=(1e-20) -declare -a LT=("1e-6") #"5e-13" "1e-13" "1e-14" "1e-15" "1e-16") # "1e-8") # "1e-8") #"1e-5") # "1e-8") -declare -a NLT=("1e-14") # "1e-9" "1e-10" "1e-14" "1e-20") -declare -a ABT=("1e-2") -declare -a RSM=("false") # "true") #"true" #"false" -declare -a SF=("9e-1") -declare -a UDS=("false" "true") -declare -a AEW=("false") # always use Eisenstat Walker, even for picard -AV=-1 -COMP="0" #"4e-12" -declare -a OS=("9e-1") # "1e-1" "1e-2" "1e-8") # "1e-3" "1e-4" "1e-5" "1e-6" "1e-6" "1e-7" "1e-8") -materialmodelnameShort="SNL" #"VP2" -materialmodelname="SNL" -phi=0 # 0 -declare -a vel=(0) #(25 50 125) # 50 125) #25 50 125) -declare -a BV=(1e19) #("1e23" "1e24" "5e24") #"1e23" "1e24" "5e24") -processes=8 - -CohesionLine="s/set List of cohesions.*/ set List of cohesions = 1e8,0/g" -PhiLine="s/set List of angles of internal friction.*/ set List of angles of internal friction = $phi,0/g" -if [ $materialmodelnameShort == "DP" ]; then - materialmodelname="drucker prager compositions" -elif [ $materialmodelnameShort == "vM" ]; then - materialmodelname="drucker prager compositions" -elif [ $materialmodelnameShort == "SNL" ]; then - materialmodelname="simple nonlinear compositions" - CohesionLine="s/ set List of cohesion of fields.*/ /g" - PhiLine="s/set List of angles of internal friction of fields.*/ /g" -elif [ $materialmodelnameShort == "VP2" ]; then - materialmodelname="viscoplastic2" -fi - -SOLVER_SHORT="itAdandSt" #"NS" #NS" -SOLVER="NS" -if [ $SOLVER_SHORT == "NS" ]; then - SOLVER="iterated Advection and Newton Stokes" -elif [ $SOLVER_SHORT == "itAdandSt" ]; then - SOLVER="iterated Advection and Stokes" -fi - - -for i_grid in "${grid[@]}" -do - for i_agrid in "${agrid[@]}" - do - for i_NSP in "${NSP[@]}" - do - for i_NSA in "${NSA[@]}" - do - for i_NLT in "${NLT[@]}" - do - for i_ABT in "${ABT[@]}" - do - for i_ST in "${ST[@]}" - do - for i_UDS in "${UDS[@]}" - do - for i_SF in "${SF[@]}" - do - for i_LT in "${LT[@]}" - do - for i_n in "${n[@]}" - do - for i_UFS in "${UFS[@]}" - do - for i_P in "${P[@]}" - do - for i_boundary_type in "${boundary_type[@]}" - do - for i_LS in "${LS[@]}" - do - for i_AEW in "${AEW[@]}" - do - for i_OS in "${OS[@]}" - do - for i_RSM in "${RSM[@]}" - do - for i_BV in "${BV[@]}" - do - for i_vel in "${vel[@]}" - do - U="7.92219116e-11" - if [ $i_vel == 25 ]; then - U="7.92219116e-11" - elif [ $i_vel == 50 ]; then - U="1.58443823e-10" - elif [ $i_vel == 125 ]; then - U="3.96109558e-10" - fi - - dirname_clean="$version""$materialmodelnameShort""_BT""$i_boundary_type""_""$SOLVER_SHORT""_ST""$i_ST""_UFS""$i_UFS""_NSP-""$i_NSP""_NSA-""$i_NSA""_C""$COMP""_g""$i_grid""_ag""$i_agrid""_AEW""$i_AEW""_UDS""$i_UDS""_SF""$i_SF""_NLT""$i_NLT""_ABT""$i_ABT""_LT""$i_LT""_mLT""$i_OS""_I""$I""_P""$i_P""_EW1""_theta1""_LS""$i_LS""_RSM""$i_RSM""_AV""$AV""_phi""$phi""_vel""$i_vel""_BV""$i_BV""_n""$i_n" - dirname="results/""$dirname_clean" - infilename="$dirname""/input.prm" - outfilename="$dirname""/output.log" - errorfilename="$dirname""/error.log" - outplotfilename="$dirname""/plot.dat" - - input_name="input_""$i_boundary_type"".prm" - - mkdir -p $dirname - - echo "$dirname" - sed \ - -e "$PhiLine" \ - -e "$CohesionLine" \ - -e "s/ set Function expression = if(x<60e3,.*/ set Function expression = if(x<60e3,$U,-$U);0/g" \ - -e "s/ set Nonlinear Newton solver switch tolerance.*/ set Nonlinear Newton solver switch tolerance = $i_ST/g" \ - -e "s/set Output directory .*/set Output directory = results\/$dirname_clean/g" \ - -e "s/ set Model name = .*/ set Model name = $materialmodelname/g" \ - -e "s/set Nonlinear solver scheme.*/set Nonlinear solver scheme = $SOLVER/g" \ - -e "s/ set Initial global refinement = .*/ set Initial global refinement = $i_grid/g" \ - -e "s/ set Initial adaptive refinement = .*/ set Initial adaptive refinement = $i_agrid/g" \ - -e "s/ set List of stress exponents of fields = .*/ set List of stress exponents of fields = $i_n, 1/g" \ - -e "s/set Stress exponent .*/ set Stress exponent = $i_n/g" \ - -e "s/set Function constants = n .*/ set Function constants = n = $i_n/g" \ - -e "s/ set Viscosity averaging p = .*/ set Viscosity averaging p = $AV/g" \ - -e "s/set Max nonlinear iterations = .*/set Max nonlinear iterations = $I /g" \ - -e "s/set Linear solver tolerance =.*/set Linear solver tolerance = $i_LT/g" \ - -e "s/set Nonlinear solver tolerance.*/set Nonlinear solver tolerance = $i_NLT/g" \ - -e "s/set Linear solver A block tolerance.*/set Linear solver A block tolerance = $i_ABT/g" \ - -e "s/ set Reference compressibility .*/ set Reference compressibility = $COMP/g" \ - -e "s/set Max pre-Newton nonlinear iterations = .*/ set Max pre-Newton nonlinear iterations = $i_P/g" \ - -e "s/set Use Newton failsafe = .*/set Use Newton failsafe = $i_UFS/g" \ - -e "s/set Stabilization preconditioner = .*/set Stabilization preconditioner = $i_NSP/g" \ - -e "s/set Stabilization velocity block = .*/set Stabilization velocity block = $i_NSA/g" \ - -e "s/set SPD safety factor = .*/set SPD safety factor = $i_SF/g" \ - -e "s/set Use deviator of strain-rate = .*/set Use deviator of strain-rate = $i_UDS/g" \ - -e "s/set Max Newton line search iterations = .*/ set Max Newton line search iterations = $i_LS/g" \ - -e "s/set Maximum linear Stokes solver tolerance =.*/set Maximum linear Stokes solver tolerance = $i_OS/g" \ - -e "s/set Use Newton residual scaling method .*/ set Use Newton residual scaling method = $i_RSM/g" \ - $input_name > "$infilename" - - nohup mpirun -np $processes ../../.././aspect $infilename > $outfilename 2>$errorfilename - - grep "Relative nonlinear residual (total" $outfilename > $outplotfilename - - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done - done -done diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t.png b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t.png index 5c466200cee..b908eeab83a 100644 Binary files a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t.png and b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t.png differ diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t_default_parameters.png b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t_default_parameters.png new file mode 100644 index 00000000000..f8d522faa0b Binary files /dev/null and b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_t_default_parameters.png differ diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v.png b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v.png index 9ab3c7f883e..3287c7b0a7a 100644 Binary files a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v.png and b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v.png differ diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v_default_parameters.png b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v_default_parameters.png new file mode 100644 index 00000000000..e54677404f7 Binary files /dev/null and b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/doc/figure_v_default_parameters.png differ diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_t.prm b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_t.prm index ca7b64af468..e1ea81d66e5 100644 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_t.prm +++ b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_t.prm @@ -1,8 +1,8 @@ # Like the poiseuille_2d.prm test and based on the nonlinear channel -# flow benchmark. This is used to test the velocity boundary conditions +# flow benchmark. This is used to test the traction boundary conditions # of the Newton Stokes solver. -set Output directory = results/b1SNL_BTt_NS_ST1e-20_UFSfalse_NSP-SPD_NSA-SPD_C0_g6_ag0_AEWfalse_UDStrue_SF9e-1_NLT1e-14_ABT1e-2_LT1e-5_mLT9e-1_I150_P5_EW1_theta1_LS0_RSMtrue_AV-1_phi0_vel0_BV1e19_n5 +set Output directory = output_t set Dimension = 2 set CFL number = 1.0 set Maximum time step = 1 @@ -11,22 +11,23 @@ set Start time = 0 set Adiabatic surface temperature = 0 set Surface pressure = 0 set Use years in output instead of seconds = false -set Nonlinear solver scheme = iterated Advection and Newton Stokes +set Nonlinear solver scheme = single Advection, iterated Newton Stokes set Max nonlinear iterations = 150 set Nonlinear solver tolerance = 1e-14 set Additional shared libraries = libsimple_nonlinear.so subsection Solver parameters subsection Stokes solver parameters - set Linear solver tolerance = 1e-5 + set Linear solver tolerance = 1e-7 + set Stokes solver type = block AMG end subsection Newton solver parameters - set Max pre-Newton nonlinear iterations = 5 + set Max pre-Newton nonlinear iterations = 3 set Nonlinear Newton solver switch tolerance = 1e-20 - set Max Newton line search iterations = 0 - set Maximum linear Stokes solver tolerance = 9e-1 - set Use Newton residual scaling method = true + set Max Newton line search iterations = 5 + set Maximum linear Stokes solver tolerance = 1e-2 + set Use Newton residual scaling method = false set Use Newton failsafe = false set Stabilization preconditioner = SPD set Stabilization velocity block = SPD @@ -35,7 +36,7 @@ end subsection Boundary temperature model set List of model names = box - set Fixed temperature boundary indicators = 2, 3 + set Fixed temperature boundary indicators = bottom, top subsection Box set Left temperature = 0 @@ -73,62 +74,58 @@ subsection Material model subsection Simple nonlinear set Minimum viscosity = 1e19 - set Maximum viscosity = 1e25 - set Stress exponent = 5 - set Viscosity averaging p = 10000 + set Maximum viscosity = 1e24 + set Stress exponent = 3 + set Viscosity averaging p = -1 set Viscosity prefactor = 1e-37 + set Use deviator of strain-rate = false end end subsection Mesh refinement set Initial adaptive refinement = 0 - set Initial global refinement = 6 + set Initial global refinement = 4 end -subsection Boundary traction model - set Prescribed traction boundary indicators = 2 y: function, 3 y: function # pressure bc: Prescribe a horizontal traction on the vertical boundaries +subsection Boundary velocity model + set Zero velocity boundary indicators = left, right + # pressure bc: Prescribe a zero vertical velocity component on the vertical boundaries + set Prescribed velocity boundary indicators = bottom x: function, top x: function subsection Function - set Variable names = x,z + set Function constants = n = 3 + set Variable names = x,y - # We want to prescribe a pressure of 2 at the left boundary - # and -2 at the right boundary. - # The traction in this case is defined as: - # tau = - pressure * normal_vector. - # On the left boundary, the outward pointing normal vector is - # (-1;0). On the right (1;0). - # Therefore: - # Left boundary: tau = - pressure(left) (-1;0) = - (2) (-1;0) = (2;0). - # Right boundary: tau = - pressure(right) (1;0) = - (-2) (1;0) = (2;0). - # Conveniently, the traction is the same on both boundaries. - set Function expression = 0;if(z>0,0,1e9) + # For velocity boundary conditions both are used, for pressure boundary conditions only the first (x) component + set Function expression = 0;(1e-37/(n+1))*((1e9/8e3)^n)*(((5e3)^(n+1))-(abs(x-(5e3))^(n+1))); end end -subsection Boundary velocity model - set Zero velocity boundary indicators = 0, 1 - set Prescribed velocity boundary indicators = 2 x: function, 3 x: function # pressure bc: Prescribe a zero vertical velocity component on the vertical boundaries - +subsection Boundary traction model + # pressure bc: Prescribe a horizontal traction on the vertical boundaries + set Prescribed traction boundary indicators = bottom y: function, top y: function subsection Function - set Function constants = n = 5 - set Variable names = x,z + set Variable names = x,y - # For velocity boundary conditions both are used, for pressure boundary conditions only the first (x) component - set Function expression = 0;(1e-37/(n+1))*((1e9/8e3)^n)*(((5e3)^(n+1))-((x-(5e3))^(n+1))); + # We want to prescribe a pressure of 1e9 at the bottom boundary + # and 0 at the top boundary. + # The traction in this case is defined as: + # tau = - pressure * normal_vector. + # On the bottom boundary, the outward pointing normal vector is + # (0;-1). On the top (0;1). + # Therefore: + # Top boundary: tau = - pressure(top) (0;1) = - (0) (0;1) = (0;0). + # Bottom boundary: tau = - pressure(bottom) (0;-1) = - (1e9) (0;-1) = (0;1e9). + # Conveniently, the traction is the same on both boundaries. + set Function expression = 0; y>4e3 ? 0 : 1e9 end end subsection Postprocess - set List of postprocessors = velocity statistics, pressure statistics, mass flux statistics, visualization #, depth average - - subsection Depth average - set List of output variables = viscosity - end + set List of postprocessors = velocity statistics, pressure statistics, mass flux statistics, visualization subsection Visualization - set List of output variables = material properties,strain rate, spd factor - - # set Output format = gnuplot + set List of output variables = material properties, strain rate set Time between graphical output = 2 end end diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_v.prm b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_v.prm index 1ceef140ad8..600b5ad0219 100644 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_v.prm +++ b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_v.prm @@ -2,7 +2,7 @@ # flow benchmark. This is used to test the velocity boundary conditions # of the Newton Stokes solver. -set Output directory = a2_itAdandSt_ST1e-20_UFSfalse_NSP-SPD_NSA-none_C0_g6_ag0_AEWfalse_SF9e-1_NLT1e-14_ABT1e-2_LT1e-5_mLT9e-1_I150_P3_EW1_theta1_LS5_RSMfalse_AV-1_n5 +set Output directory = output_v set Dimension = 2 set CFL number = 1.0 set Maximum time step = 1 @@ -11,31 +11,32 @@ set Start time = 0 set Adiabatic surface temperature = 0 set Surface pressure = 0 set Use years in output instead of seconds = false -set Nonlinear solver scheme = iterated Advection and Stokes #Newton Stokes +set Nonlinear solver scheme = single Advection, iterated Newton Stokes set Max nonlinear iterations = 150 set Nonlinear solver tolerance = 1e-14 set Additional shared libraries = libsimple_nonlinear.so subsection Solver parameters subsection Stokes solver parameters - set Linear solver tolerance = 1e-15 + set Linear solver tolerance = 1e-7 + set Stokes solver type = block AMG end subsection Newton solver parameters set Max pre-Newton nonlinear iterations = 3 set Nonlinear Newton solver switch tolerance = 1e-20 set Max Newton line search iterations = 5 - set Maximum linear Stokes solver tolerance = 9e-1 + set Maximum linear Stokes solver tolerance = 1e-2 set Use Newton residual scaling method = false set Use Newton failsafe = false set Stabilization preconditioner = SPD - set Stabilization velocity block = none + set Stabilization velocity block = SPD end end subsection Boundary temperature model set List of model names = box - set Fixed temperature boundary indicators = 2, 3 + set Fixed temperature boundary indicators = bottom, top subsection Box set Left temperature = 0 @@ -74,8 +75,8 @@ subsection Material model subsection Simple nonlinear set Minimum viscosity = 1e19 set Maximum viscosity = 1e24 - set Stress exponent = 5 - set Viscosity averaging p = 10000 + set Stress exponent = 3 + set Viscosity averaging p = -1 set Viscosity prefactor = 1e-37 set Use deviator of strain-rate = false end @@ -83,16 +84,16 @@ end subsection Mesh refinement set Initial adaptive refinement = 0 - set Initial global refinement = 6 + set Initial global refinement = 4 end subsection Boundary velocity model - set Zero velocity boundary indicators = 0, 1 - set Prescribed velocity boundary indicators = 2: function, 3: function # velocity bc + set Zero velocity boundary indicators = left, right + set Prescribed velocity boundary indicators = bottom: function, top: function subsection Function - set Function constants = n = 5 - set Variable names = x,z + set Function constants = n = 3 + set Variable names = x,y # For velocity boundary conditions both are used, for pressure boundary conditions only the first (x) component set Function expression = 0;(1e-37/(n+1))*((1e9/8e3)^n)*(((5e3)^(n+1))-(abs(x-(5e3))^(n+1))); @@ -100,34 +101,30 @@ subsection Boundary velocity model end subsection Boundary traction model + set Prescribed traction boundary indicators = + subsection Function - set Variable names = x,z + set Variable names = x,y - # We want to prescribe a pressure of 2 at the left boundary - # and -2 at the right boundary. + # We want to prescribe a pressure of 1e9 at the bottom boundary + # and 0 at the top boundary. # The traction in this case is defined as: # tau = - pressure * normal_vector. - # On the left boundary, the outward pointing normal vector is - # (-1;0). On the right (1;0). + # On the bottom boundary, the outward pointing normal vector is + # (0;-1). On the top (0;1). # Therefore: - # Left boundary: tau = - pressure(left) (-1;0) = - (2) (-1;0) = (2;0). - # Right boundary: tau = - pressure(right) (1;0) = - (-2) (1;0) = (2;0). + # Top boundary: tau = - pressure(top) (0;1) = - (0) (0;1) = (0;0). + # Bottom boundary: tau = - pressure(bottom) (0;-1) = - (1e9) (0;-1) = (0;1e9). # Conveniently, the traction is the same on both boundaries. - set Function expression = 0;if(z>0,0,1e9) + set Function expression = 0; y>4e3 ? 0 : 1e9 end end subsection Postprocess - set List of postprocessors = velocity statistics, pressure statistics, mass flux statistics, visualization #, depth average - - subsection Depth average - set List of output variables = viscosity - end + set List of postprocessors = velocity statistics, pressure statistics, mass flux statistics, visualization subsection Visualization - set List of output variables = material properties,strain rate#, spd factor - - # set Output format = gnuplot + set List of output variables = material properties, strain rate set Time between graphical output = 2 end end diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/metabash.sh b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/metabash.sh deleted file mode 100755 index de6a6250772..00000000000 --- a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/metabash.sh +++ /dev/null @@ -1,81 +0,0 @@ -#!/bin/bash - -version="b1" -declare -a grid=(4) -declare -a NSP=("SPD") # "symmetric" "none") #Use Newton stabilization preconditioner -declare -a NSA=("SPD" "none") # "symmetric" "none") #Use Newton stabilization A block -declare -a UFS=("truefalse") #"false") # use failsafe -declare -a RSM=("true" "false") # use residual scaling -declare -a LS=(0-100) #0-5-10) #(10) # line search -declare -a COMP=(0) # copmressibiltiy -declare -a phi=(0) -declare -a ST=("1e-20") # Nonlinear Newton solver switch tolerance -declare -a UDS=("falsetrue") # Use deviator of strain-rate -declare -a SF=("9e-1") # "9.9e-1") # "9.99e-1" "1") -declare -a AEW=("false") # "true") # always use Eisenstat Walker, even for picard -n="3" -boundary_type="t" -processes=4 - -for i_grid in "${grid[@]}" -do - for i_NSP in "${NSP[@]}" - do - for i_NSA in "${NSA[@]}" - do - for i_LS in "${LS[@]}" - do - for i_RSM in "${RSM[@]}" - do - for i_UFS in "${UFS[@]}" - do - for i_COMP in "${COMP[@]}" - do - for i_phi in "${phi[@]}" - do - for i_ST in "${ST[@]}" - do - for i_UDS in "${UDS[@]}" - do - for i_SF in "${SF[@]}" - do - for i_AEW in "${AEW[@]}" - do - output_name="$version""_g""$i_grid""_UFS""$i_UFS""_NSP-""$i_NSP""_NSA-""$i_NSA""_LS""$i_LS""_RSM""$i_RSM""_phi""$i_phi""_COMP""$i_COMP""_ST""$i_ST""_UDS""$i_UDS""_SF""$i_SF""_AEW""$i_AEW""_n""$n""_BT""$boundary_type" - log_file="auto_logs/""$output_name"".log" - bash_auto_file="bash_autos/bash_""$output_name"".sh" - echo "$output_name" - sed \ - -e "s/declare -a grid.*/declare -a grid=($i_grid)/g" \ - -e "s/declare -a UFS.*/declare -a UFS=(\"false\" \"true\")/g" \ - -e "s/declare -a NSP.*/declare -a NSP=($i_NSP)/g" \ - -e "s/declare -a NSA.*/declare -a NSA=($i_NSA)/g" \ - -e "s/declare -a ST.*/declare -a ST=($i_ST)/g" \ - -e "s/declare -a UDS.*/declare -a UDS=(\"false\" \"true\")/g" \ - -e "s/declare -a SF.*/declare -a SF=($i_SF)/g" \ - -e "s/declare -a boundary type.*/declare -a boundary_type=($i_boundary_type)/g" \ - -e "s/declare -a AEW.*/declare -a AEW=($i_AEW)/g" \ - -e "s/declare -a LS.*/declare -a LS=(0 100)/g" \ - -e "s/version=.*/declare -a version=\"$version\"/g" \ - -e "s/declare -a RSM.*/declare -a RSM=(\"$i_RSM\")/g" \ - -e "s/phi=.*/phi=\"$i_phi\"/g" \ - -e "s/declare -a n=.*/declare -a n=\"$n\"/g" \ - -e "s/COMP=.*/COMP=\"$i_COMP\"/g" \ - -e "s/processes=.*/processes=\"$processes\"/g" \ - bash.sh > $bash_auto_file - - nohup bash ./$bash_auto_file > $log_file 2>&1 & - sleep 2 - - done - done - done - done - done - done - done - done - done - done - done -done diff --git a/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/plot_t.gnuplot b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/plot_t.gnuplot new file mode 100644 index 00000000000..6f0627eaeca --- /dev/null +++ b/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/plot_t.gnuplot @@ -0,0 +1,85 @@ +set terminal png size 1500,1000 enhanced font 'Verdana,14' +set output 'figure_t.png' +set yrange [1e-14:1] +set xrange [0:60] +set format y '%g' +set logscale y +set multiplot title "With a minimum linear tolerance of 1e-8, a stress exponent of 3 and pressure boundary conditions" font 'Verdana, 20' +set size 0.51,0.49 +set origin 0,0.48 +unset key +set title "Without stabilization of the velocity block and without the RSM" +plot \ +'results/BT_t_singleAdvectioniteratedStokes/plot.dat' u 9:10 w l lw 2 lc rgb 'purple' t 'Picard', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_150_LS_0_RSM_false/plot.dat' u 9:10 w l lw 2 lc rgb 'green' t 'DC Picard', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_5_LS_0_RSM_false/plot.dat' u 9:10 w l dt 1 lw 2 lc rgb 'blue' t 'Newton solver, 5 Picard, no line search', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_5_LS_100_RSM_false/plot.dat' u 9:10 w p dt 3 lw 2 lc rgb 'blue' pt 7 t 'Newton solver, 5 Picard, line search', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_10_LS_0_RSM_false/plot.dat' u 9:10 w l dt 1 lw 2 lc rgb 'orange' t 'Newton solver, 10 Picard, no line search', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_10_LS_100_RSM_false/plot.dat' u 9:10 w p dt 3 lw 2 lc rgb 'orange' pt 7 ps 1 t 'Newton solver, 10 Picard, line search', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_15_LS_0_RSM_false/plot.dat' u 9:10 w l dt 1 lw 2 lc rgb 'red' t 'Newton solver, 15 Picard, no line search', \ +'results/BT_t_singleAdvectioniteratedNewtonStokes_stabilization_none_mLT_1e-8_P_15_LS_100_RSM_false/plot.dat' u 9:10 w p dt 3 lw 2 lc rgb 'red' pt 7 ps 1 t 'Newton solver, 15 Picard, line search', \ +' $dirname/$infile + + nohup mpirun -np $processes ./aspect $dirname/$infile > $dirname/$outfile 2>$dirname/$errorfile + grep "Relative nonlinear residual " $dirname/$outfile > $dirname/$plotfile + + # Then run all the Newton solver models + for stabilization in "SPD" "none"; do + for picard_iterations in 5 10 15 150; do + for line_search_iterations in 0 5 100; do + for max_linear_tolerance in "1e-8" "1e-2"; do + for residual_scaling_method in "false" "true"; do + dirname_clean="BT_${boundary_type}_singleAdvectioniteratedNewtonStokes_stabilization_${stabilization}_mLT_${max_linear_tolerance}_P_${picard_iterations}_LS_${line_search_iterations}_RSM_${residual_scaling_method}" + dirname="results/${dirname_clean}" + mkdir -p $dirname + echo "Starting $dirname" + + sed \ + -e "s/set Output directory .*/set Output directory = results\/$dirname_clean/g" \ + -e "s/set Max pre-Newton nonlinear iterations = .*/set Max pre-Newton nonlinear iterations = $picard_iterations/g" \ + -e "s/set Stabilization preconditioner = .*/set Stabilization preconditioner = $stabilization/g" \ + -e "s/set Stabilization velocity block = .*/set Stabilization velocity block = $stabilization/g" \ + -e "s/set Max Newton line search iterations = .*/set Max Newton line search iterations = $line_search_iterations/g" \ + -e "s/set Maximum linear Stokes solver tolerance =.*/set Maximum linear Stokes solver tolerance = $max_linear_tolerance/g" \ + -e "s/set Use Newton residual scaling method .*/set Use Newton residual scaling method = $residual_scaling_method/g" \ + $infile > $dirname/$infile + + nohup mpirun -np $processes ./aspect $dirname/$infile > $dirname/$outfile 2>$dirname/$errorfile + grep "Relative nonlinear residual " $dirname/$outfile > $dirname/$plotfile + done + done + done + done + done +done diff --git a/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/screen-output b/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/screen-output index d8b72fe9c9a..96fe4b5ec3c 100644 --- a/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/screen-output +++ b/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/screen-output @@ -1,86 +1,69 @@ Loading shared library <./libnonlinear_channel_flow_tractions_Newton_Stokes_GMG.debug.so> -Number of active cells: 4,096 (on 7 levels) -Number of degrees of freedom: 54,148 (33,282+4,225+16,641) +Number of active cells: 256 (on 5 levels) +Number of degrees of freedom: 3,556 (2,178+289+1,089) *** Timestep 0: t=0 seconds, dt=0 seconds Skipping temperature solve because RHS is zero. - Initial Newton Stokes residual = 9.30238e+11, v = 9.30238e+11, p = 0 + Initial Newton Stokes residual = 1.85171e+12, v = 1.85171e+12, p = 0 - Solving Stokes system (GMG)... 8+0 iterations. - Newton system information: Norm of the rhs: 9.30238e+11 - Relative nonlinear residuals (temperature, Stokes system): 0, 1 - Relative nonlinear residual (total system) after nonlinear iteration 1: 1 + Solving Stokes system (GMG)... 13+0 iterations. + Newton system information: Norm of the rhs: 1.85171e+12 + Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1 - Skipping temperature solve because RHS is zero. - Solving Stokes system (GMG)... 10+0 iterations. - Newton system information: Norm of the rhs: 8.35883e+10, Derivative scaling factor: 0 - Relative nonlinear residuals (temperature, Stokes system): 0, 0.0898569 - Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0898569 + Solving Stokes system (GMG)... 12+0 iterations. + Newton system information: Norm of the rhs: 3.04031e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.16419 - Skipping temperature solve because RHS is zero. - Solving Stokes system (GMG)... 7+0 iterations. - Newton system information: Norm of the rhs: 1.08226e+10, Derivative scaling factor: 0 - Relative nonlinear residuals (temperature, Stokes system): 0, 0.0116342 - Relative nonlinear residual (total system) after nonlinear iteration 3: 0.0116342 + Solving Stokes system (GMG)... 14+0 iterations. + Newton system information: Norm of the rhs: 2.64066e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 0.142607 - Skipping temperature solve because RHS is zero. - Solving Stokes system (GMG)... 4+0 iterations. - Newton system information: Norm of the rhs: 105099, Derivative scaling factor: 0 - Relative nonlinear residuals (temperature, Stokes system): 0, 1.12981e-07 - Relative nonlinear residual (total system) after nonlinear iteration 4: 1.12981e-07 + Switching from defect correction form of Picard to the Newton solver scheme. + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Solving Stokes system (GMG)... 10+0 iterations. + Newton system information: Norm of the rhs: 2.06072e+11, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 0.111288 - Skipping temperature solve because RHS is zero. - Solving Stokes system (GMG)... 8+0 iterations. - Newton system information: Norm of the rhs: 0.672766, Derivative scaling factor: 0 - Relative nonlinear residuals (temperature, Stokes system): 0, 7.23219e-13 - Relative nonlinear residual (total system) after nonlinear iteration 5: 7.23219e-13 + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Solving Stokes system (GMG)... 11+0 iterations. + Newton system information: Norm of the rhs: 1.49894e+11, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 0.080949 - Skipping temperature solve because RHS is zero. - Switching from defect correction form of Picard to the Newton solver scheme. - The linear solver tolerance is set to 1.12984e-06. Stabilization Preconditioner is SPD and A block is SPD. + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. Solving Stokes system (GMG)... 11+0 iterations. - Newton system information: Norm of the rhs: 0.0951387, Derivative scaling factor: 0 - Relative nonlinear residuals (temperature, Stokes system): 0, 1.02274e-13 - Relative nonlinear residual (total system) after nonlinear iteration 6: 1.02274e-13 + Newton system information: Norm of the rhs: 7.71743e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 6: 0.0416774 - Skipping temperature solve because RHS is zero. - The linear solver tolerance is set to 0.121142. Stabilization Preconditioner is SPD and A block is SPD. - Solving Stokes system (GMG)... 1+0 iterations. - Newton system information: Norm of the rhs: 0.0946666, Derivative scaling factor: 0.858586 - Relative nonlinear residuals (temperature, Stokes system): 0, 1.01766e-13 - Relative nonlinear residual (total system) after nonlinear iteration 7: 1.01766e-13 + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Solving Stokes system (GMG)... 12+0 iterations. + Newton system information: Norm of the rhs: 1.81606e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 7: 0.00980748 - Skipping temperature solve because RHS is zero. - The linear solver tolerance is set to 0.9. Stabilization Preconditioner is SPD and A block is SPD. - Solving Stokes system (GMG)... 1+0 iterations. - Newton system information: Norm of the rhs: 0.0887962, Derivative scaling factor: 0.859288 - Relative nonlinear residuals (temperature, Stokes system): 0, 9.54553e-14 - Relative nonlinear residual (total system) after nonlinear iteration 8: 9.54553e-14 + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Solving Stokes system (GMG)... 12+0 iterations. + Newton system information: Norm of the rhs: 4.89728e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 8: 0.00264474 - Skipping temperature solve because RHS is zero. - The linear solver tolerance is set to 0.89746. Stabilization Preconditioner is SPD and A block is SPD. - Solving Stokes system (GMG)... 1+0 iterations. - Newton system information: Norm of the rhs: 0.0890302, Derivative scaling factor: 0.868013 - Relative nonlinear residuals (temperature, Stokes system): 0, 9.57069e-14 - Relative nonlinear residual (total system) after nonlinear iteration 9: 9.57069e-14 + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Solving Stokes system (GMG)... 5+0 iterations. + Newton system information: Norm of the rhs: 2.64989e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 9: 0.00143105 - Skipping temperature solve because RHS is zero. - The linear solver tolerance is set to 0.9. Stabilization Preconditioner is SPD and A block is SPD. - Solving Stokes system (GMG)... 1+0 iterations. - Newton system information: Norm of the rhs: 0.0982373, Derivative scaling factor: 0.867665 - Relative nonlinear residuals (temperature, Stokes system): 0, 1.05604e-13 - Relative nonlinear residual (total system) after nonlinear iteration 10: 1.05604e-13 + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Solving Stokes system (GMG)... 7+0 iterations. + Newton system information: Norm of the rhs: 2.42373e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 10: 0.00130892 WARNING: The nonlinear solver in the current timestep failed to converge. Acting according to the parameter 'Nonlinear solver failure strategy'... Continuing to the next timestep even though solution is not fully converged. Postprocessing: - RMS, max velocity: 1.14e-07 m/s, 1.56e-07 m/s - Pressure min/avg/max: -7.077e-06 Pa, 5e+08 Pa, 1e+09 Pa - Mass fluxes through boundary parts: 0 kg/s, 0 kg/s, -3.473 kg/s, 3.473 kg/s + RMS, max velocity: 2.55e-08 m/s, 3.02e-08 m/s + Pressure min/avg/max: -296.2 Pa, 5e+08 Pa, 1e+09 Pa + Mass fluxes through boundary parts: 0 kg/s, 0 kg/s, -0.8043 kg/s, 0.8043 kg/s Writing graphical output: output-nonlinear_channel_flow_tractions_Newton_Stokes_GMG/solution/solution-00000 Termination requested by criterion: end time diff --git a/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/statistics b/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/statistics index d35c196824c..cb9fd5f4c3e 100644 --- a/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/statistics +++ b/tests/nonlinear_channel_flow_tractions_Newton_Stokes_GMG/statistics @@ -19,4 +19,4 @@ # 19: Outward mass flux through boundary with indicator 2 ("bottom") (kg/s) # 20: Outward mass flux through boundary with indicator 3 ("top") (kg/s) # 21: Visualization file name -0 0.000000000000e+00 0.000000000000e+00 4096 37507 16641 10 0 42 62 62 1.14108866e-07 1.56248062e-07 -7.07741689e-06 5.00000000e+08 1.00000000e+09 0.00000000e+00 0.00000000e+00 -3.47275156e+00 3.47275156e+00 output-nonlinear_channel_flow_tractions_Newton_Stokes_GMG/solution/solution-00000 +0 0.000000000000e+00 0.000000000000e+00 256 2467 1089 10 0 97 117 117 2.54638657e-08 3.02425152e-08 -2.96236519e+02 4.99999462e+08 9.99999881e+08 0.00000000e+00 0.00000000e+00 -8.04290283e-01 8.04290290e-01 output-nonlinear_channel_flow_tractions_Newton_Stokes_GMG/solution/solution-00000 diff --git a/tests/nonlinear_channel_flow_tractions_benchmark.cc b/tests/nonlinear_channel_flow_tractions_benchmark.cc new file mode 100644 index 00000000000..b0c1b8e736e --- /dev/null +++ b/tests/nonlinear_channel_flow_tractions_benchmark.cc @@ -0,0 +1 @@ +#include "../benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/simple_nonlinear.cc" diff --git a/tests/nonlinear_channel_flow_tractions_benchmark.prm b/tests/nonlinear_channel_flow_tractions_benchmark.prm new file mode 100644 index 00000000000..f5af47a7150 --- /dev/null +++ b/tests/nonlinear_channel_flow_tractions_benchmark.prm @@ -0,0 +1,4 @@ +# Test the benchmark parameter file +# benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_t.prm + +include $ASPECT_SOURCE_DIR/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_t.prm diff --git a/tests/nonlinear_channel_flow_tractions_benchmark/screen-output b/tests/nonlinear_channel_flow_tractions_benchmark/screen-output new file mode 100644 index 00000000000..2857f0b5cf3 --- /dev/null +++ b/tests/nonlinear_channel_flow_tractions_benchmark/screen-output @@ -0,0 +1,1469 @@ + +Loading shared library <./libnonlinear_channel_flow_tractions_benchmark.debug.so> + +Number of active cells: 256 (on 5 levels) +Number of degrees of freedom: 3,556 (2,178+289+1,089) + +*** Timestep 0: t=0 seconds, dt=0 seconds + Skipping temperature solve because RHS is zero. + Initial Newton Stokes residual = 1.85171e+12, v = 1.85171e+12, p = 0 + + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 38+0 iterations. + Newton system information: Norm of the rhs: 1.85171e+12 + Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1 + + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 49+0 iterations. + Newton system information: Norm of the rhs: 3.1208e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.168536 + + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 83+0 iterations. + Newton system information: Norm of the rhs: 2.72888e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 0.147371 + + Switching from defect correction form of Picard to the Newton solver scheme. + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 54+0 iterations. + Newton system information: Norm of the rhs: 1.79996e+11, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 0.0972052 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 66+0 iterations. + Newton system information: Norm of the rhs: 7.85851e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 0.0424393 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 96+0 iterations. + Newton system information: Norm of the rhs: 1.67375e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 6: 0.00903893 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 98+0 iterations. + Newton system information: Norm of the rhs: 7.74748e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 7: 0.00418397 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 82+0 iterations. + Line search iteration 0, with norm of the rhs 8.0074e+09 and going to 7.7467e+09, relative residual: 0.00432433 + Newton system information: Norm of the rhs: 5.28972e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 8: 0.00285667 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 79+0 iterations. + Newton system information: Norm of the rhs: 4.26618e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 9: 0.00230391 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 43+0 iterations. + Line search iteration 0, with norm of the rhs 5.28821e+09 and going to 4.26575e+09, relative residual: 0.00285586 + Newton system information: Norm of the rhs: 3.57208e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 10: 0.00192908 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 38+0 iterations. + Newton system information: Norm of the rhs: 2.5051e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 11: 0.00135286 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 42+0 iterations. + Line search iteration 0, with norm of the rhs 3.06636e+09 and going to 2.50484e+09, relative residual: 0.00165596 + Newton system information: Norm of the rhs: 2.11287e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 12: 0.00114104 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 49+0 iterations. + Newton system information: Norm of the rhs: 2.03291e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 13: 0.00109786 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 40+0 iterations. + Line search iteration 0, with norm of the rhs 2.56906e+09 and going to 2.03271e+09, relative residual: 0.0013874 + Newton system information: Norm of the rhs: 1.73886e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 14: 0.000939057 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 29+0 iterations. + Newton system information: Norm of the rhs: 1.18235e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 15: 0.000638522 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 39+0 iterations. + Newton system information: Norm of the rhs: 8.87341e+08, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 16: 0.000479202 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 39+0 iterations. + Newton system information: Norm of the rhs: 8.81726e+06, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 17: 4.76169e-06 + + The linear solver tolerance is set to 1.46493e-07. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 644+0 iterations. + Newton system information: Norm of the rhs: 208.524, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 18: 1.12612e-10 + + The linear solver tolerance is set to 2.35095e-05. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 354+0 iterations. + Newton system information: Norm of the rhs: 127.164, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 19: 6.86738e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 7+0 iterations. + Newton system information: Norm of the rhs: 106.258, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 20: 5.73835e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Newton system information: Norm of the rhs: 98.9143, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 21: 5.34179e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Newton system information: Norm of the rhs: 89.702, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 22: 4.84429e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 7+0 iterations. + Newton system information: Norm of the rhs: 69.2659, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 23: 3.74065e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Newton system information: Norm of the rhs: 54.7388, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 24: 2.95613e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 11+0 iterations. + Line search iteration 0, with norm of the rhs 141.149 and going to 54.7312, relative residual: 7.62263e-11 + Line search iteration 1, with norm of the rhs 106.165 and going to 54.733, relative residual: 5.73334e-11 + Line search iteration 2, with norm of the rhs 71.4401 and going to 54.733, relative residual: 3.85807e-11 + Line search iteration 3, with norm of the rhs 69.8671 and going to 54.733, relative residual: 3.77312e-11 + Line search iteration 4, with norm of the rhs 101.045 and going to 54.733, relative residual: 5.45685e-11 + Newton system information: Norm of the rhs: 54.1083, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 25: 2.92208e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 71.9683 and going to 54.1031, relative residual: 3.88659e-11 + Line search iteration 1, with norm of the rhs 77.0575 and going to 54.1049, relative residual: 4.16143e-11 + Line search iteration 2, with norm of the rhs 65.6731 and going to 54.1049, relative residual: 3.54662e-11 + Line search iteration 3, with norm of the rhs 69.3664 and going to 54.1049, relative residual: 3.74608e-11 + Line search iteration 4, with norm of the rhs 63.5005 and going to 54.1049, relative residual: 3.42929e-11 + Newton system information: Norm of the rhs: 67.1661, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 26: 3.62725e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 133.064 and going to 67.1593, relative residual: 7.186e-11 + Line search iteration 1, with norm of the rhs 97.0918 and going to 67.1615, relative residual: 5.24337e-11 + Line search iteration 2, with norm of the rhs 83.8918 and going to 67.1615, relative residual: 4.53051e-11 + Line search iteration 3, with norm of the rhs 100.683 and going to 67.1615, relative residual: 5.4373e-11 + Line search iteration 4, with norm of the rhs 99.3966 and going to 67.1615, relative residual: 5.36784e-11 + Newton system information: Norm of the rhs: 65.0426, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 27: 3.51258e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 114.288 and going to 65.0361, relative residual: 6.17204e-11 + Line search iteration 1, with norm of the rhs 125.705 and going to 65.0382, relative residual: 6.78858e-11 + Line search iteration 2, with norm of the rhs 95.5727 and going to 65.0382, relative residual: 5.16133e-11 + Line search iteration 3, with norm of the rhs 114.966 and going to 65.0382, relative residual: 6.20866e-11 + Newton system information: Norm of the rhs: 58.6722, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 28: 3.16855e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 92.2561 and going to 58.6662, relative residual: 4.98222e-11 + Line search iteration 1, with norm of the rhs 72.2866 and going to 58.6682, relative residual: 3.90378e-11 + Line search iteration 2, with norm of the rhs 102.111 and going to 58.6682, relative residual: 5.51443e-11 + Line search iteration 3, with norm of the rhs 77.399 and going to 58.6682, relative residual: 4.17987e-11 + Line search iteration 4, with norm of the rhs 60.0485 and going to 58.6682, relative residual: 3.24287e-11 + Newton system information: Norm of the rhs: 60.0719, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 29: 3.24414e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 96.94 and going to 60.0659, relative residual: 5.23517e-11 + Line search iteration 1, with norm of the rhs 74.6639 and going to 60.0679, relative residual: 4.03216e-11 + Line search iteration 2, with norm of the rhs 83.5532 and going to 60.0679, relative residual: 4.51222e-11 + Line search iteration 3, with norm of the rhs 68.6286 and going to 60.0679, relative residual: 3.70623e-11 + Newton system information: Norm of the rhs: 58.4881, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 30: 3.15861e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 109.245 and going to 58.4821, relative residual: 5.8997e-11 + Line search iteration 1, with norm of the rhs 76.6844 and going to 58.4841, relative residual: 4.14128e-11 + Line search iteration 2, with norm of the rhs 115.313 and going to 58.4841, relative residual: 6.22738e-11 + Line search iteration 3, with norm of the rhs 76.4071 and going to 58.4841, relative residual: 4.12631e-11 + Newton system information: Norm of the rhs: 57.1635, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 31: 3.08707e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 92.0602 and going to 57.1577, relative residual: 4.97164e-11 + Line search iteration 1, with norm of the rhs 70.0682 and going to 57.1596, relative residual: 3.78398e-11 + Line search iteration 2, with norm of the rhs 68.5812 and going to 57.1596, relative residual: 3.70367e-11 + Line search iteration 3, with norm of the rhs 67.5424 and going to 57.1596, relative residual: 3.64758e-11 + Line search iteration 4, with norm of the rhs 65.4275 and going to 57.1596, relative residual: 3.53336e-11 + Newton system information: Norm of the rhs: 56.0573, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 32: 3.02733e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 11+0 iterations. + Line search iteration 0, with norm of the rhs 87.6513 and going to 56.0516, relative residual: 4.73354e-11 + Line search iteration 1, with norm of the rhs 78.7423 and going to 56.0534, relative residual: 4.25242e-11 + Line search iteration 2, with norm of the rhs 81.7539 and going to 56.0534, relative residual: 4.41505e-11 + Line search iteration 3, with norm of the rhs 65.9692 and going to 56.0534, relative residual: 3.56261e-11 + Line search iteration 4, with norm of the rhs 68.4846 and going to 56.0534, relative residual: 3.69846e-11 + Newton system information: Norm of the rhs: 65.1964, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 33: 3.52088e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 84.5136 and going to 65.1899, relative residual: 4.56409e-11 + Line search iteration 1, with norm of the rhs 74.4175 and going to 65.192, relative residual: 4.01886e-11 + Line search iteration 2, with norm of the rhs 67.6422 and going to 65.192, relative residual: 3.65297e-11 + Newton system information: Norm of the rhs: 58.7119, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 34: 3.17069e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 11+0 iterations. + Newton system information: Norm of the rhs: 54.9922, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 35: 2.96981e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 81.6035 and going to 54.9863, relative residual: 4.40693e-11 + Line search iteration 1, with norm of the rhs 61.7278 and going to 54.9882, relative residual: 3.33356e-11 + Line search iteration 2, with norm of the rhs 85.8706 and going to 54.9882, relative residual: 4.63738e-11 + Line search iteration 3, with norm of the rhs 60.5786 and going to 54.9882, relative residual: 3.2715e-11 + Newton system information: Norm of the rhs: 53.1169, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 36: 2.86854e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 71.8091 and going to 53.112, relative residual: 3.878e-11 + Line search iteration 1, with norm of the rhs 86.3138 and going to 53.1137, relative residual: 4.66131e-11 + Line search iteration 2, with norm of the rhs 65.467 and going to 53.1137, relative residual: 3.5355e-11 + Line search iteration 3, with norm of the rhs 60.2022 and going to 53.1137, relative residual: 3.25117e-11 + Line search iteration 4, with norm of the rhs 66.3223 and going to 53.1137, relative residual: 3.58168e-11 + Newton system information: Norm of the rhs: 53.6058, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 37: 2.89494e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 144.169 and going to 53.6009, relative residual: 7.78575e-11 + Line search iteration 1, with norm of the rhs 73.235 and going to 53.6027, relative residual: 3.955e-11 + Line search iteration 2, with norm of the rhs 93.7335 and going to 53.6027, relative residual: 5.06201e-11 + Line search iteration 3, with norm of the rhs 54.5159 and going to 53.6027, relative residual: 2.94409e-11 + Line search iteration 4, with norm of the rhs 69.7257 and going to 53.6027, relative residual: 3.76548e-11 + Newton system information: Norm of the rhs: 60.8851, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 38: 3.28805e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 76.2832 and going to 60.8794, relative residual: 4.11961e-11 + Line search iteration 1, with norm of the rhs 74.7197 and going to 60.8814, relative residual: 4.03518e-11 + Line search iteration 2, with norm of the rhs 71.9084 and going to 60.8814, relative residual: 3.88336e-11 + Line search iteration 3, with norm of the rhs 69.446 and going to 60.8814, relative residual: 3.75037e-11 + Newton system information: Norm of the rhs: 56.8168, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 39: 3.06835e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 75.2377 and going to 56.8123, relative residual: 4.06315e-11 + Line search iteration 1, with norm of the rhs 95.2209 and going to 56.8142, relative residual: 5.14233e-11 + Line search iteration 2, with norm of the rhs 62.6654 and going to 56.8142, relative residual: 3.3842e-11 + Line search iteration 3, with norm of the rhs 67.9339 and going to 56.8142, relative residual: 3.66872e-11 + Line search iteration 4, with norm of the rhs 62.9425 and going to 56.8142, relative residual: 3.39916e-11 + Newton system information: Norm of the rhs: 57.6812, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 40: 3.11503e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 94.6227 and going to 57.6762, relative residual: 5.11003e-11 + Line search iteration 1, with norm of the rhs 69.7469 and going to 57.6781, relative residual: 3.76663e-11 + Line search iteration 2, with norm of the rhs 66.2357 and going to 57.6781, relative residual: 3.57701e-11 + Line search iteration 3, with norm of the rhs 72.7964 and going to 57.6781, relative residual: 3.93131e-11 + Line search iteration 4, with norm of the rhs 85.0942 and going to 57.6781, relative residual: 4.59545e-11 + Newton system information: Norm of the rhs: 60.0174, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 41: 3.24119e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 150.275 and going to 60.0119, relative residual: 8.1155e-11 + Line search iteration 1, with norm of the rhs 71.8798 and going to 60.0139, relative residual: 3.88181e-11 + Line search iteration 2, with norm of the rhs 116.808 and going to 60.0139, relative residual: 6.30814e-11 + Line search iteration 3, with norm of the rhs 62.3617 and going to 60.0139, relative residual: 3.3678e-11 + Line search iteration 4, with norm of the rhs 69.1471 and going to 60.0139, relative residual: 3.73423e-11 + Newton system information: Norm of the rhs: 65.1053, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 42: 3.51596e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 72.5498 and going to 65.0991, relative residual: 3.918e-11 + Line search iteration 1, with norm of the rhs 83.8274 and going to 65.1013, relative residual: 4.52704e-11 + Newton system information: Norm of the rhs: 63.1255, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 43: 3.40904e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 122.909 and going to 63.119, relative residual: 6.63763e-11 + Line search iteration 1, with norm of the rhs 66.3361 and going to 63.1211, relative residual: 3.58243e-11 + Line search iteration 2, with norm of the rhs 68.5549 and going to 63.1211, relative residual: 3.70226e-11 + Line search iteration 3, with norm of the rhs 71.5265 and going to 63.1211, relative residual: 3.86273e-11 + Line search iteration 4, with norm of the rhs 77.1731 and going to 63.1211, relative residual: 4.16767e-11 + Newton system information: Norm of the rhs: 64.3583, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 44: 3.47562e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 94.2433 and going to 64.3519, relative residual: 5.08953e-11 + Line search iteration 1, with norm of the rhs 87.6775 and going to 64.354, relative residual: 4.73495e-11 + Line search iteration 2, with norm of the rhs 77.3101 and going to 64.354, relative residual: 4.17507e-11 + Newton system information: Norm of the rhs: 62.5799, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 45: 3.37958e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 93.1778 and going to 62.5733, relative residual: 5.032e-11 + Line search iteration 1, with norm of the rhs 66.2924 and going to 62.5754, relative residual: 3.58007e-11 + Line search iteration 2, with norm of the rhs 64.9035 and going to 62.5754, relative residual: 3.50506e-11 + Line search iteration 3, with norm of the rhs 88.2454 and going to 62.5754, relative residual: 4.76563e-11 + Line search iteration 4, with norm of the rhs 74.0135 and going to 62.5754, relative residual: 3.99704e-11 + Newton system information: Norm of the rhs: 60.465, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 46: 3.26536e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 98.9979 and going to 60.4589, relative residual: 5.3463e-11 + Line search iteration 1, with norm of the rhs 77.1671 and going to 60.4609, relative residual: 4.16735e-11 + Line search iteration 2, with norm of the rhs 77.328 and going to 60.4609, relative residual: 4.17604e-11 + Line search iteration 3, with norm of the rhs 65.0095 and going to 60.4609, relative residual: 3.51079e-11 + Line search iteration 4, with norm of the rhs 62.7622 and going to 60.4609, relative residual: 3.38942e-11 + Newton system information: Norm of the rhs: 62.2316, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 47: 3.36077e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 75.29 and going to 62.2253, relative residual: 4.06598e-11 + Line search iteration 1, with norm of the rhs 84.2954 and going to 62.2274, relative residual: 4.55231e-11 + Line search iteration 2, with norm of the rhs 77.2832 and going to 62.2274, relative residual: 4.17362e-11 + Line search iteration 3, with norm of the rhs 65.3405 and going to 62.2274, relative residual: 3.52866e-11 + Newton system information: Norm of the rhs: 61.4032, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 48: 3.31603e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 65.7954 and going to 61.397, relative residual: 3.55323e-11 + Line search iteration 1, with norm of the rhs 63.462 and going to 61.3991, relative residual: 3.42721e-11 + Newton system information: Norm of the rhs: 56.4803, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 49: 3.05018e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 78.2009 and going to 56.4744, relative residual: 4.22318e-11 + Line search iteration 1, with norm of the rhs 86.2755 and going to 56.4763, relative residual: 4.65924e-11 + Line search iteration 2, with norm of the rhs 63.3983 and going to 56.4763, relative residual: 3.42378e-11 + Line search iteration 3, with norm of the rhs 64.1747 and going to 56.4763, relative residual: 3.4657e-11 + Newton system information: Norm of the rhs: 55.606, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 50: 3.00296e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 76.5426 and going to 55.6005, relative residual: 4.13363e-11 + Line search iteration 1, with norm of the rhs 65.2374 and going to 55.6023, relative residual: 3.5231e-11 + Newton system information: Norm of the rhs: 54.7265, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 51: 2.95546e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 88.6393 and going to 54.7211, relative residual: 4.78689e-11 + Line search iteration 1, with norm of the rhs 60.7199 and going to 54.723, relative residual: 3.27913e-11 + Line search iteration 2, with norm of the rhs 75.2792 and going to 54.723, relative residual: 4.06539e-11 + Newton system information: Norm of the rhs: 48.9375, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 52: 2.64283e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 59.1857 and going to 48.9326, relative residual: 3.19627e-11 + Newton system information: Norm of the rhs: 48.5501, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 53: 2.62191e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 9+0 iterations. + Line search iteration 0, with norm of the rhs 59.7764 and going to 48.5448, relative residual: 3.22818e-11 + Line search iteration 1, with norm of the rhs 49.3143 and going to 48.5465, relative residual: 2.66318e-11 + Line search iteration 2, with norm of the rhs 50.2002 and going to 48.5465, relative residual: 2.71102e-11 + Line search iteration 3, with norm of the rhs 83.9602 and going to 48.5465, relative residual: 4.53421e-11 + Newton system information: Norm of the rhs: 43.674, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 54: 2.35858e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 8+0 iterations. + Line search iteration 0, with norm of the rhs 73.7423 and going to 43.6693, relative residual: 3.9824e-11 + Line search iteration 1, with norm of the rhs 69.3376 and going to 43.6708, relative residual: 3.74452e-11 + Line search iteration 2, with norm of the rhs 84.1309 and going to 43.6708, relative residual: 4.54342e-11 + Newton system information: Norm of the rhs: 40.6249, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 55: 2.19392e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 72.3782 and going to 40.6202, relative residual: 3.90873e-11 + Line search iteration 1, with norm of the rhs 82.9002 and going to 40.6216, relative residual: 4.47696e-11 + Line search iteration 2, with norm of the rhs 52.2679 and going to 40.6216, relative residual: 2.82269e-11 + Newton system information: Norm of the rhs: 39.3707, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 56: 2.12619e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 11+0 iterations. + Line search iteration 0, with norm of the rhs 80.6077 and going to 39.3663, relative residual: 4.35316e-11 + Line search iteration 1, with norm of the rhs 98.0594 and going to 39.3676, relative residual: 5.29562e-11 + Line search iteration 2, with norm of the rhs 47.3573 and going to 39.3676, relative residual: 2.55749e-11 + Line search iteration 3, with norm of the rhs 45.8856 and going to 39.3676, relative residual: 2.47802e-11 + Line search iteration 4, with norm of the rhs 39.5445 and going to 39.3676, relative residual: 2.13557e-11 + Newton system information: Norm of the rhs: 39.2894, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 57: 2.12179e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 10+0 iterations. + Line search iteration 0, with norm of the rhs 90.9045 and going to 39.2853, relative residual: 4.90923e-11 + Line search iteration 1, with norm of the rhs 58.7041 and going to 39.2866, relative residual: 3.17027e-11 + Line search iteration 2, with norm of the rhs 93.4486 and going to 39.2866, relative residual: 5.04662e-11 + Line search iteration 3, with norm of the rhs 67.5639 and going to 39.2866, relative residual: 3.64873e-11 + Line search iteration 4, with norm of the rhs 68.5101 and going to 39.2866, relative residual: 3.69983e-11 + Newton system information: Norm of the rhs: 39.0123, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 58: 2.10683e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 75.0831 and going to 39.0085, relative residual: 4.0548e-11 + Line search iteration 1, with norm of the rhs 79.8856 and going to 39.0098, relative residual: 4.31416e-11 + Line search iteration 2, with norm of the rhs 67.9521 and going to 39.0098, relative residual: 3.6697e-11 + Line search iteration 3, with norm of the rhs 44.87 and going to 39.0098, relative residual: 2.42317e-11 + Line search iteration 4, with norm of the rhs 39.7556 and going to 39.0098, relative residual: 2.14697e-11 + Newton system information: Norm of the rhs: 39.3808, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 59: 2.12673e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 64.3868 and going to 39.3768, relative residual: 3.47716e-11 + Line search iteration 1, with norm of the rhs 63.2695 and going to 39.3781, relative residual: 3.41682e-11 + Line search iteration 2, with norm of the rhs 41.9213 and going to 39.3781, relative residual: 2.26393e-11 + Line search iteration 3, with norm of the rhs 39.7958 and going to 39.3781, relative residual: 2.14914e-11 + Line search iteration 4, with norm of the rhs 40.5822 and going to 39.3781, relative residual: 2.19161e-11 + Newton system information: Norm of the rhs: 37.6724, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 60: 2.03447e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 68.4433 and going to 37.6686, relative residual: 3.69623e-11 + Line search iteration 1, with norm of the rhs 47.0032 and going to 37.6698, relative residual: 2.53837e-11 + Line search iteration 2, with norm of the rhs 40.4853 and going to 37.6698, relative residual: 2.18638e-11 + Line search iteration 3, with norm of the rhs 43.3423 and going to 37.6698, relative residual: 2.34067e-11 + Line search iteration 4, with norm of the rhs 42.0367 and going to 37.6698, relative residual: 2.27016e-11 + Newton system information: Norm of the rhs: 41.0233, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 61: 2.21543e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 79.2213 and going to 41.0192, relative residual: 4.27828e-11 + Line search iteration 1, with norm of the rhs 69.7864 and going to 41.0205, relative residual: 3.76876e-11 + Line search iteration 2, with norm of the rhs 68.2362 and going to 41.0205, relative residual: 3.68504e-11 + Line search iteration 3, with norm of the rhs 62.0235 and going to 41.0205, relative residual: 3.34953e-11 + Newton system information: Norm of the rhs: 40.6054, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 62: 2.19287e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 84.0837 and going to 40.6013, relative residual: 4.54088e-11 + Line search iteration 1, with norm of the rhs 68.8634 and going to 40.6026, relative residual: 3.71892e-11 + Line search iteration 2, with norm of the rhs 65.0652 and going to 40.6026, relative residual: 3.51379e-11 + Line search iteration 3, with norm of the rhs 69.7101 and going to 40.6026, relative residual: 3.76464e-11 + Line search iteration 4, with norm of the rhs 44.4574 and going to 40.6026, relative residual: 2.40088e-11 + Newton system information: Norm of the rhs: 41.427, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 63: 2.23723e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 91.393 and going to 41.4228, relative residual: 4.93561e-11 + Line search iteration 1, with norm of the rhs 70.6241 and going to 41.4242, relative residual: 3.814e-11 + Line search iteration 2, with norm of the rhs 51.6393 and going to 41.4242, relative residual: 2.78874e-11 + Line search iteration 3, with norm of the rhs 42.21 and going to 41.4242, relative residual: 2.27952e-11 + Newton system information: Norm of the rhs: 40.4071, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 64: 2.18215e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 73.1141 and going to 40.4028, relative residual: 3.94847e-11 + Line search iteration 1, with norm of the rhs 71.5397 and going to 40.4042, relative residual: 3.86345e-11 + Line search iteration 2, with norm of the rhs 51.5476 and going to 40.4042, relative residual: 2.78379e-11 + Line search iteration 3, with norm of the rhs 49.1214 and going to 40.4042, relative residual: 2.65276e-11 + Line search iteration 4, with norm of the rhs 46.1603 and going to 40.4042, relative residual: 2.49285e-11 + Newton system information: Norm of the rhs: 43.1157, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 65: 2.32843e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 82.2179 and going to 43.1113, relative residual: 4.44011e-11 + Line search iteration 1, with norm of the rhs 56.2959 and going to 43.1127, relative residual: 3.04021e-11 + Newton system information: Norm of the rhs: 42.7743, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 66: 2.30999e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 67.2775 and going to 42.7711, relative residual: 3.63327e-11 + Line search iteration 1, with norm of the rhs 75.8909 and going to 42.7725, relative residual: 4.09843e-11 + Line search iteration 2, with norm of the rhs 72.0866 and going to 42.7725, relative residual: 3.89298e-11 + Line search iteration 3, with norm of the rhs 109.571 and going to 42.7725, relative residual: 5.91731e-11 + Line search iteration 4, with norm of the rhs 43.1516 and going to 42.7725, relative residual: 2.33037e-11 + Newton system information: Norm of the rhs: 43.0511, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 67: 2.32494e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 102.624 and going to 43.0471, relative residual: 5.54215e-11 + Line search iteration 1, with norm of the rhs 72.5989 and going to 43.0486, relative residual: 3.92065e-11 + Line search iteration 2, with norm of the rhs 76.3727 and going to 43.0486, relative residual: 4.12445e-11 + Line search iteration 3, with norm of the rhs 68.9536 and going to 43.0486, relative residual: 3.72379e-11 + Newton system information: Norm of the rhs: 41.6223, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 68: 2.24778e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 72.9563 and going to 41.6188, relative residual: 3.93995e-11 + Line search iteration 1, with norm of the rhs 49.0433 and going to 41.6202, relative residual: 2.64855e-11 + Newton system information: Norm of the rhs: 37.8174, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 69: 2.0423e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 88.1504 and going to 37.8151, relative residual: 4.76049e-11 + Line search iteration 1, with norm of the rhs 83.3936 and going to 37.8163, relative residual: 4.50361e-11 + Line search iteration 2, with norm of the rhs 42.4439 and going to 37.8163, relative residual: 2.29215e-11 + Line search iteration 3, with norm of the rhs 41.8248 and going to 37.8163, relative residual: 2.25871e-11 + Line search iteration 4, with norm of the rhs 39.979 and going to 37.8163, relative residual: 2.15903e-11 + Newton system information: Norm of the rhs: 39.7937, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 70: 2.14903e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 82.7086 and going to 39.7901, relative residual: 4.46661e-11 + Line search iteration 1, with norm of the rhs 60.3488 and going to 39.7914, relative residual: 3.25909e-11 + Line search iteration 2, with norm of the rhs 91.6488 and going to 39.7914, relative residual: 4.94942e-11 + Line search iteration 3, with norm of the rhs 72.9879 and going to 39.7914, relative residual: 3.94166e-11 + Line search iteration 4, with norm of the rhs 77.2484 and going to 39.7914, relative residual: 4.17174e-11 + Newton system information: Norm of the rhs: 38.357, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 71: 2.07144e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 45.3792 and going to 38.3536, relative residual: 2.45067e-11 + Line search iteration 1, with norm of the rhs 40.6032 and going to 38.3549, relative residual: 2.19274e-11 + Line search iteration 2, with norm of the rhs 39.3992 and going to 38.3549, relative residual: 2.12772e-11 + Line search iteration 3, with norm of the rhs 38.4839 and going to 38.3549, relative residual: 2.07829e-11 + Line search iteration 4, with norm of the rhs 39.4088 and going to 38.3549, relative residual: 2.12824e-11 + Newton system information: Norm of the rhs: 42.5414, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 72: 2.29741e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 57.9868 and going to 42.5374, relative residual: 3.13153e-11 + Line search iteration 1, with norm of the rhs 42.5593 and going to 42.5389, relative residual: 2.29838e-11 + Newton system information: Norm of the rhs: 36.67, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 73: 1.98033e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 81.3676 and going to 36.6682, relative residual: 4.39419e-11 + Line search iteration 1, with norm of the rhs 75.9867 and going to 36.6695, relative residual: 4.10361e-11 + Line search iteration 2, with norm of the rhs 81.3429 and going to 36.6695, relative residual: 4.39286e-11 + Line search iteration 3, with norm of the rhs 43.1671 and going to 36.6695, relative residual: 2.3312e-11 + Line search iteration 4, with norm of the rhs 38.3608 and going to 36.6695, relative residual: 2.07164e-11 + Newton system information: Norm of the rhs: 37.6486, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 74: 2.03318e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 98.714 and going to 37.6452, relative residual: 5.33097e-11 + Line search iteration 1, with norm of the rhs 74.1047 and going to 37.6464, relative residual: 4.00196e-11 + Line search iteration 2, with norm of the rhs 89.8679 and going to 37.6464, relative residual: 4.85325e-11 + Line search iteration 3, with norm of the rhs 47.6565 and going to 37.6464, relative residual: 2.57365e-11 + Line search iteration 4, with norm of the rhs 38.854 and going to 37.6464, relative residual: 2.09828e-11 + Newton system information: Norm of the rhs: 37.1852, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 75: 2.00816e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 107.098 and going to 37.1818, relative residual: 5.78375e-11 + Line search iteration 1, with norm of the rhs 64.8796 and going to 37.1831, relative residual: 3.50377e-11 + Line search iteration 2, with norm of the rhs 74.6353 and going to 37.1831, relative residual: 4.03062e-11 + Line search iteration 3, with norm of the rhs 73.4419 and going to 37.1831, relative residual: 3.96617e-11 + Line search iteration 4, with norm of the rhs 83.9079 and going to 37.1831, relative residual: 4.53138e-11 + Newton system information: Norm of the rhs: 38.5168, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 76: 2.08007e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 57.5162 and going to 38.5132, relative residual: 3.10612e-11 + Line search iteration 1, with norm of the rhs 51.9497 and going to 38.5145, relative residual: 2.80551e-11 + Line search iteration 2, with norm of the rhs 44.6212 and going to 38.5145, relative residual: 2.40974e-11 + Line search iteration 3, with norm of the rhs 42.7153 and going to 38.5145, relative residual: 2.30681e-11 + Line search iteration 4, with norm of the rhs 43.3514 and going to 38.5145, relative residual: 2.34116e-11 + Newton system information: Norm of the rhs: 42.1656, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 77: 2.27712e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 81.6363 and going to 42.1617, relative residual: 4.40871e-11 + Line search iteration 1, with norm of the rhs 61.3829 and going to 42.1631, relative residual: 3.31493e-11 + Line search iteration 2, with norm of the rhs 96.4107 and going to 42.1631, relative residual: 5.20658e-11 + Line search iteration 3, with norm of the rhs 70.8959 and going to 42.1631, relative residual: 3.82868e-11 + Line search iteration 4, with norm of the rhs 69.6808 and going to 42.1631, relative residual: 3.76305e-11 + Newton system information: Norm of the rhs: 35.4704, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 78: 1.91555e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 14+0 iterations. + Line search iteration 0, with norm of the rhs 75.8833 and going to 35.4668, relative residual: 4.09802e-11 + Line search iteration 1, with norm of the rhs 85.8802 and going to 35.468, relative residual: 4.63789e-11 + Line search iteration 2, with norm of the rhs 82.9517 and going to 35.468, relative residual: 4.47974e-11 + Line search iteration 3, with norm of the rhs 39.681 and going to 35.468, relative residual: 2.14294e-11 + Line search iteration 4, with norm of the rhs 36.1738 and going to 35.468, relative residual: 1.95354e-11 + Newton system information: Norm of the rhs: 34.8771, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 79: 1.88351e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 83.8286 and going to 34.8732, relative residual: 4.5271e-11 + Line search iteration 1, with norm of the rhs 57.4612 and going to 34.8744, relative residual: 3.10315e-11 + Line search iteration 2, with norm of the rhs 104.359 and going to 34.8744, relative residual: 5.63584e-11 + Line search iteration 3, with norm of the rhs 68.7145 and going to 34.8744, relative residual: 3.71088e-11 + Line search iteration 4, with norm of the rhs 69.2729 and going to 34.8744, relative residual: 3.74103e-11 + Newton system information: Norm of the rhs: 37.3304, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 80: 2.016e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 106.203 and going to 37.3266, relative residual: 5.7354e-11 + Line search iteration 1, with norm of the rhs 71.939 and going to 37.3278, relative residual: 3.88501e-11 + Line search iteration 2, with norm of the rhs 95.1808 and going to 37.3278, relative residual: 5.14017e-11 + Line search iteration 3, with norm of the rhs 47.9516 and going to 37.3278, relative residual: 2.58959e-11 + Newton system information: Norm of the rhs: 36.4437, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 81: 1.96812e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 71.0836 and going to 36.4399, relative residual: 3.83881e-11 + Line search iteration 1, with norm of the rhs 66.8347 and going to 36.4411, relative residual: 3.60935e-11 + Line search iteration 2, with norm of the rhs 67.2466 and going to 36.4411, relative residual: 3.6316e-11 + Line search iteration 3, with norm of the rhs 66.2435 and going to 36.4411, relative residual: 3.57743e-11 + Newton system information: Norm of the rhs: 35.2602, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 82: 1.9042e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 62.3608 and going to 35.2561, relative residual: 3.36775e-11 + Line search iteration 1, with norm of the rhs 109.63 and going to 35.2572, relative residual: 5.92049e-11 + Line search iteration 2, with norm of the rhs 69.319 and going to 35.2572, relative residual: 3.74352e-11 + Line search iteration 3, with norm of the rhs 75.3861 and going to 35.2572, relative residual: 4.07117e-11 + Line search iteration 4, with norm of the rhs 42.6247 and going to 35.2572, relative residual: 2.30191e-11 + Newton system information: Norm of the rhs: 34.5363, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 83: 1.86511e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 73.3972 and going to 34.5329, relative residual: 3.96376e-11 + Line search iteration 1, with norm of the rhs 61.7538 and going to 34.5341, relative residual: 3.33497e-11 + Line search iteration 2, with norm of the rhs 44.5177 and going to 34.5341, relative residual: 2.40415e-11 + Line search iteration 3, with norm of the rhs 41.7849 and going to 34.5341, relative residual: 2.25656e-11 + Line search iteration 4, with norm of the rhs 40.8201 and going to 34.5341, relative residual: 2.20446e-11 + Newton system information: Norm of the rhs: 37.997, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 84: 2.052e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 64.4429 and going to 37.993, relative residual: 3.48019e-11 + Line search iteration 1, with norm of the rhs 87.2256 and going to 37.9943, relative residual: 4.71055e-11 + Line search iteration 2, with norm of the rhs 68.1327 and going to 37.9943, relative residual: 3.67945e-11 + Line search iteration 3, with norm of the rhs 87.1168 and going to 37.9943, relative residual: 4.70468e-11 + Newton system information: Norm of the rhs: 36.7504, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 85: 1.98468e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 62.9142 and going to 36.7467, relative residual: 3.39763e-11 + Line search iteration 1, with norm of the rhs 37.1457 and going to 36.7479, relative residual: 2.00602e-11 + Line search iteration 2, with norm of the rhs 39.508 and going to 36.7479, relative residual: 2.1336e-11 + Line search iteration 3, with norm of the rhs 39.7138 and going to 36.7479, relative residual: 2.14471e-11 + Newton system information: Norm of the rhs: 34.8699, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 86: 1.88312e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 90.2589 and going to 34.8661, relative residual: 4.87436e-11 + Line search iteration 1, with norm of the rhs 66.4401 and going to 34.8673, relative residual: 3.58804e-11 + Line search iteration 2, with norm of the rhs 64.9266 and going to 34.8673, relative residual: 3.50631e-11 + Line search iteration 3, with norm of the rhs 68.4206 and going to 34.8673, relative residual: 3.695e-11 + Line search iteration 4, with norm of the rhs 51.5351 and going to 34.8673, relative residual: 2.78311e-11 + Newton system information: Norm of the rhs: 39.5031, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 87: 2.13333e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 71.3004 and going to 39.4991, relative residual: 3.85052e-11 + Line search iteration 1, with norm of the rhs 69.4359 and going to 39.5004, relative residual: 3.74983e-11 + Line search iteration 2, with norm of the rhs 67.5279 and going to 39.5004, relative residual: 3.64679e-11 + Line search iteration 3, with norm of the rhs 67.7788 and going to 39.5004, relative residual: 3.66034e-11 + Line search iteration 4, with norm of the rhs 46.1019 and going to 39.5004, relative residual: 2.4897e-11 + Newton system information: Norm of the rhs: 35.8257, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 88: 1.93474e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 56.6352 and going to 35.8218, relative residual: 3.05854e-11 + Line search iteration 1, with norm of the rhs 46.232 and going to 35.823, relative residual: 2.49673e-11 + Line search iteration 2, with norm of the rhs 38.0422 and going to 35.823, relative residual: 2.05444e-11 + Newton system information: Norm of the rhs: 35.1107, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 89: 1.89613e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 75.091 and going to 35.1072, relative residual: 4.05523e-11 + Line search iteration 1, with norm of the rhs 67.3964 and going to 35.1083, relative residual: 3.63969e-11 + Line search iteration 2, with norm of the rhs 71.6995 and going to 35.1083, relative residual: 3.87207e-11 + Line search iteration 3, with norm of the rhs 88.7585 and going to 35.1083, relative residual: 4.79334e-11 + Line search iteration 4, with norm of the rhs 39.95 and going to 35.1083, relative residual: 2.15747e-11 + Newton system information: Norm of the rhs: 35.2, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 90: 1.90095e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 70.7744 and going to 35.1965, relative residual: 3.82212e-11 + Line search iteration 1, with norm of the rhs 65.7289 and going to 35.1977, relative residual: 3.54964e-11 + Line search iteration 2, with norm of the rhs 44.4796 and going to 35.1977, relative residual: 2.40209e-11 + Line search iteration 3, with norm of the rhs 42.5438 and going to 35.1977, relative residual: 2.29755e-11 + Line search iteration 4, with norm of the rhs 44.6914 and going to 35.1977, relative residual: 2.41352e-11 + Newton system information: Norm of the rhs: 39.5889, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 91: 2.13797e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 76.8182 and going to 39.5849, relative residual: 4.14851e-11 + Line search iteration 1, with norm of the rhs 89.809 and going to 39.5863, relative residual: 4.85006e-11 + Line search iteration 2, with norm of the rhs 73.7514 and going to 39.5863, relative residual: 3.98289e-11 + Line search iteration 3, with norm of the rhs 39.9213 and going to 39.5863, relative residual: 2.15592e-11 + Line search iteration 4, with norm of the rhs 41.1427 and going to 39.5863, relative residual: 2.22188e-11 + Newton system information: Norm of the rhs: 40.2836, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 92: 2.17549e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 54.306 and going to 40.2795, relative residual: 2.93275e-11 + Line search iteration 1, with norm of the rhs 78.2058 and going to 40.2808, relative residual: 4.22344e-11 + Line search iteration 2, with norm of the rhs 68.593 and going to 40.2808, relative residual: 3.70431e-11 + Line search iteration 3, with norm of the rhs 76.7391 and going to 40.2808, relative residual: 4.14423e-11 + Newton system information: Norm of the rhs: 38.1026, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 93: 2.0577e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 74.4648 and going to 38.0997, relative residual: 4.02141e-11 + Line search iteration 1, with norm of the rhs 74.896 and going to 38.101, relative residual: 4.0447e-11 + Line search iteration 2, with norm of the rhs 65.6456 and going to 38.101, relative residual: 3.54514e-11 + Line search iteration 3, with norm of the rhs 88.5321 and going to 38.101, relative residual: 4.78111e-11 + Newton system information: Norm of the rhs: 36.8331, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 94: 1.98914e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 66.0896 and going to 36.8305, relative residual: 3.56912e-11 + Line search iteration 1, with norm of the rhs 77.7208 and going to 36.8318, relative residual: 4.19725e-11 + Line search iteration 2, with norm of the rhs 66.3358 and going to 36.8318, relative residual: 3.58241e-11 + Line search iteration 3, with norm of the rhs 77.5337 and going to 36.8318, relative residual: 4.18715e-11 + Line search iteration 4, with norm of the rhs 39.3178 and going to 36.8318, relative residual: 2.12333e-11 + Newton system information: Norm of the rhs: 38.9547, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 95: 2.10372e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 85.0978 and going to 38.9515, relative residual: 4.59564e-11 + Line search iteration 1, with norm of the rhs 63.6179 and going to 38.9528, relative residual: 3.43563e-11 + Line search iteration 2, with norm of the rhs 73.161 and going to 38.9528, relative residual: 3.951e-11 + Line search iteration 3, with norm of the rhs 68.2586 and going to 38.9528, relative residual: 3.68625e-11 + Line search iteration 4, with norm of the rhs 39.8399 and going to 38.9528, relative residual: 2.15152e-11 + Newton system information: Norm of the rhs: 40.0938, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 96: 2.16523e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 84.1499 and going to 40.0905, relative residual: 4.54445e-11 + Line search iteration 1, with norm of the rhs 42.8209 and going to 40.0919, relative residual: 2.31251e-11 + Line search iteration 2, with norm of the rhs 43.1362 and going to 40.0919, relative residual: 2.32954e-11 + Newton system information: Norm of the rhs: 39.5961, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 97: 2.13836e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 44.5671 and going to 39.5936, relative residual: 2.40681e-11 + Line search iteration 1, with norm of the rhs 40.2651 and going to 39.595, relative residual: 2.17448e-11 + Line search iteration 2, with norm of the rhs 44.0597 and going to 39.595, relative residual: 2.37941e-11 + Line search iteration 3, with norm of the rhs 41.8524 and going to 39.595, relative residual: 2.2602e-11 + Line search iteration 4, with norm of the rhs 41.2988 and going to 39.595, relative residual: 2.23031e-11 + Newton system information: Norm of the rhs: 39.0284, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 98: 2.1077e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 53.489 and going to 39.0252, relative residual: 2.88863e-11 + Line search iteration 1, with norm of the rhs 48.1275 and going to 39.0265, relative residual: 2.59909e-11 + Line search iteration 2, with norm of the rhs 55.1171 and going to 39.0265, relative residual: 2.97655e-11 + Newton system information: Norm of the rhs: 37.3156, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 99: 2.0152e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 51.2138 and going to 37.3113, relative residual: 2.76576e-11 + Line search iteration 1, with norm of the rhs 72.8551 and going to 37.3126, relative residual: 3.93448e-11 + Line search iteration 2, with norm of the rhs 71.1698 and going to 37.3126, relative residual: 3.84347e-11 + Line search iteration 3, with norm of the rhs 102.885 and going to 37.3126, relative residual: 5.55623e-11 + Line search iteration 4, with norm of the rhs 39.4346 and going to 37.3126, relative residual: 2.12963e-11 + Newton system information: Norm of the rhs: 36.105, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 100: 1.94982e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 100.326 and going to 36.1012, relative residual: 5.41801e-11 + Line search iteration 1, with norm of the rhs 71.4637 and going to 36.1024, relative residual: 3.85934e-11 + Line search iteration 2, with norm of the rhs 73.319 and going to 36.1024, relative residual: 3.95954e-11 + Line search iteration 3, with norm of the rhs 60.1541 and going to 36.1024, relative residual: 3.24857e-11 + Newton system information: Norm of the rhs: 34.9283, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 101: 1.88627e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 15+0 iterations. + Line search iteration 0, with norm of the rhs 52.7335 and going to 34.9245, relative residual: 2.84783e-11 + Line search iteration 1, with norm of the rhs 86.7657 and going to 34.9256, relative residual: 4.68571e-11 + Line search iteration 2, with norm of the rhs 59.3485 and going to 34.9256, relative residual: 3.20507e-11 + Line search iteration 3, with norm of the rhs 36.0839 and going to 34.9256, relative residual: 1.94868e-11 + Line search iteration 4, with norm of the rhs 35.8131 and going to 34.9256, relative residual: 1.93406e-11 + Newton system information: Norm of the rhs: 37.7581, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 102: 2.0391e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 90.7695 and going to 37.7541, relative residual: 4.90193e-11 + Line search iteration 1, with norm of the rhs 45.2294 and going to 37.7554, relative residual: 2.44258e-11 + Line search iteration 2, with norm of the rhs 42.3479 and going to 37.7554, relative residual: 2.28697e-11 + Line search iteration 3, with norm of the rhs 39.6788 and going to 37.7554, relative residual: 2.14282e-11 + Newton system information: Norm of the rhs: 36.8394, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 103: 1.98948e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 114.812 and going to 36.8356, relative residual: 6.20033e-11 + Line search iteration 1, with norm of the rhs 48.45 and going to 36.8368, relative residual: 2.61651e-11 + Line search iteration 2, with norm of the rhs 66.6491 and going to 36.8368, relative residual: 3.59933e-11 + Line search iteration 3, with norm of the rhs 66.2091 and going to 36.8368, relative residual: 3.57557e-11 + Line search iteration 4, with norm of the rhs 97.1006 and going to 36.8368, relative residual: 5.24384e-11 + Newton system information: Norm of the rhs: 35.2928, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 104: 1.90596e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 82.53 and going to 35.289, relative residual: 4.45697e-11 + Line search iteration 1, with norm of the rhs 112.881 and going to 35.2902, relative residual: 6.09604e-11 + Line search iteration 2, with norm of the rhs 69.9667 and going to 35.2902, relative residual: 3.7785e-11 + Line search iteration 3, with norm of the rhs 66.7375 and going to 35.2902, relative residual: 3.6041e-11 + Line search iteration 4, with norm of the rhs 49.242 and going to 35.2902, relative residual: 2.65927e-11 + Newton system information: Norm of the rhs: 34.2684, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 105: 1.85064e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 80.6002 and going to 34.2648, relative residual: 4.35275e-11 + Line search iteration 1, with norm of the rhs 69.8051 and going to 34.2659, relative residual: 3.76977e-11 + Line search iteration 2, with norm of the rhs 84.8134 and going to 34.2659, relative residual: 4.58028e-11 + Line search iteration 3, with norm of the rhs 37.1202 and going to 34.2659, relative residual: 2.00465e-11 + Newton system information: Norm of the rhs: 34.0558, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 106: 1.83916e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 55.8874 and going to 34.0522, relative residual: 3.01816e-11 + Line search iteration 1, with norm of the rhs 47.2654 and going to 34.0533, relative residual: 2.55253e-11 + Line search iteration 2, with norm of the rhs 43.0006 and going to 34.0533, relative residual: 2.32221e-11 + Line search iteration 3, with norm of the rhs 39.9245 and going to 34.0533, relative residual: 2.15609e-11 + Line search iteration 4, with norm of the rhs 34.2763 and going to 34.0533, relative residual: 1.85107e-11 + Newton system information: Norm of the rhs: 38.0155, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 107: 2.053e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 81.2628 and going to 38.0115, relative residual: 4.38854e-11 + Line search iteration 1, with norm of the rhs 72.9474 and going to 38.0128, relative residual: 3.93946e-11 + Line search iteration 2, with norm of the rhs 91.3199 and going to 38.0128, relative residual: 4.93166e-11 + Line search iteration 3, with norm of the rhs 40.4797 and going to 38.0128, relative residual: 2.18608e-11 + Line search iteration 4, with norm of the rhs 39.7675 and going to 38.0128, relative residual: 2.14761e-11 + Newton system information: Norm of the rhs: 35.9312, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 108: 1.94044e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 57.4398 and going to 35.9274, relative residual: 3.10199e-11 + Line search iteration 1, with norm of the rhs 69.2502 and going to 35.9286, relative residual: 3.7398e-11 + Line search iteration 2, with norm of the rhs 71.7658 and going to 35.9286, relative residual: 3.87565e-11 + Line search iteration 3, with norm of the rhs 100.548 and going to 35.9286, relative residual: 5.43e-11 + Line search iteration 4, with norm of the rhs 37.0107 and going to 35.9286, relative residual: 1.99873e-11 + Newton system information: Norm of the rhs: 33.8617, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 109: 1.82867e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 15+0 iterations. + Line search iteration 0, with norm of the rhs 47.8819 and going to 33.8582, relative residual: 2.58582e-11 + Line search iteration 1, with norm of the rhs 48.9744 and going to 33.8593, relative residual: 2.64482e-11 + Line search iteration 2, with norm of the rhs 42.8747 and going to 33.8593, relative residual: 2.31541e-11 + Line search iteration 3, with norm of the rhs 38.9096 and going to 33.8593, relative residual: 2.10128e-11 + Line search iteration 4, with norm of the rhs 36.3422 and going to 33.8593, relative residual: 1.96263e-11 + Newton system information: Norm of the rhs: 35.2788, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 110: 1.9052e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 101.138 and going to 35.275, relative residual: 5.46188e-11 + Line search iteration 1, with norm of the rhs 70.1928 and going to 35.2762, relative residual: 3.79071e-11 + Line search iteration 2, with norm of the rhs 72.2708 and going to 35.2762, relative residual: 3.90293e-11 + Line search iteration 3, with norm of the rhs 54.3793 and going to 35.2762, relative residual: 2.93671e-11 + Line search iteration 4, with norm of the rhs 35.2933 and going to 35.2762, relative residual: 1.90599e-11 + Newton system information: Norm of the rhs: 37.3706, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 111: 2.01817e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 14+0 iterations. + Line search iteration 0, with norm of the rhs 69.4558 and going to 37.3667, relative residual: 3.75091e-11 + Line search iteration 1, with norm of the rhs 78.6517 and going to 37.368, relative residual: 4.24752e-11 + Line search iteration 2, with norm of the rhs 54.0108 and going to 37.368, relative residual: 2.91681e-11 + Newton system information: Norm of the rhs: 37.321, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 112: 2.01549e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 120.808 and going to 37.3167, relative residual: 6.52415e-11 + Line search iteration 1, with norm of the rhs 84.0117 and going to 37.318, relative residual: 4.53699e-11 + Line search iteration 2, with norm of the rhs 68.3736 and going to 37.318, relative residual: 3.69246e-11 + Line search iteration 3, with norm of the rhs 68.3453 and going to 37.318, relative residual: 3.69093e-11 + Line search iteration 4, with norm of the rhs 62.925 and going to 37.318, relative residual: 3.39821e-11 + Newton system information: Norm of the rhs: 34.1257, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 113: 1.84293e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 73.8333 and going to 34.1221, relative residual: 3.98731e-11 + Line search iteration 1, with norm of the rhs 68.2812 and going to 34.1232, relative residual: 3.68747e-11 + Line search iteration 2, with norm of the rhs 68.9694 and going to 34.1232, relative residual: 3.72464e-11 + Line search iteration 3, with norm of the rhs 95.6304 and going to 34.1232, relative residual: 5.16444e-11 + Line search iteration 4, with norm of the rhs 36.4743 and going to 34.1232, relative residual: 1.96977e-11 + Newton system information: Norm of the rhs: 36.0034, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 114: 1.94434e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 14+0 iterations. + Line search iteration 0, with norm of the rhs 49.0154 and going to 35.9996, relative residual: 2.64704e-11 + Line search iteration 1, with norm of the rhs 41.9851 and going to 36.0008, relative residual: 2.26737e-11 + Line search iteration 2, with norm of the rhs 40.3134 and going to 36.0008, relative residual: 2.17709e-11 + Line search iteration 3, with norm of the rhs 39.6716 and going to 36.0008, relative residual: 2.14243e-11 + Line search iteration 4, with norm of the rhs 36.1947 and going to 36.0008, relative residual: 1.95467e-11 + Newton system information: Norm of the rhs: 33.8937, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 115: 1.8304e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 68.8794 and going to 33.8901, relative residual: 3.71978e-11 + Line search iteration 1, with norm of the rhs 114.655 and going to 33.8912, relative residual: 6.19183e-11 + Line search iteration 2, with norm of the rhs 65.6379 and going to 33.8912, relative residual: 3.54473e-11 + Line search iteration 3, with norm of the rhs 70.5936 and going to 33.8912, relative residual: 3.81235e-11 + Line search iteration 4, with norm of the rhs 42.2377 and going to 33.8912, relative residual: 2.28101e-11 + Newton system information: Norm of the rhs: 39.0339, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 116: 2.10799e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 80.8561 and going to 39.0299, relative residual: 4.36657e-11 + Line search iteration 1, with norm of the rhs 91.2027 and going to 39.0312, relative residual: 4.92533e-11 + Line search iteration 2, with norm of the rhs 90.3795 and going to 39.0312, relative residual: 4.88088e-11 + Newton system information: Norm of the rhs: 35.2729, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 117: 1.90489e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 14+0 iterations. + Line search iteration 0, with norm of the rhs 71.6347 and going to 35.2686, relative residual: 3.86858e-11 + Line search iteration 1, with norm of the rhs 90.2357 and going to 35.2698, relative residual: 4.87311e-11 + Line search iteration 2, with norm of the rhs 68.544 and going to 35.2698, relative residual: 3.70166e-11 + Line search iteration 3, with norm of the rhs 66.051 and going to 35.2698, relative residual: 3.56703e-11 + Line search iteration 4, with norm of the rhs 47.4205 and going to 35.2698, relative residual: 2.56091e-11 + Newton system information: Norm of the rhs: 36.3758, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 118: 1.96444e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 50.0703 and going to 36.3721, relative residual: 2.70401e-11 + Line search iteration 1, with norm of the rhs 39.6714 and going to 36.3733, relative residual: 2.14243e-11 + Line search iteration 2, with norm of the rhs 43.3118 and going to 36.3733, relative residual: 2.33902e-11 + Line search iteration 3, with norm of the rhs 41.7323 and going to 36.3733, relative residual: 2.25372e-11 + Line search iteration 4, with norm of the rhs 39.0486 and going to 36.3733, relative residual: 2.10879e-11 + Newton system information: Norm of the rhs: 37.9837, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 119: 2.05128e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 100.172 and going to 37.9797, relative residual: 5.40969e-11 + Line search iteration 1, with norm of the rhs 56.9397 and going to 37.9809, relative residual: 3.07498e-11 + Line search iteration 2, with norm of the rhs 48.964 and going to 37.9809, relative residual: 2.64426e-11 + Line search iteration 3, with norm of the rhs 42.6019 and going to 37.9809, relative residual: 2.30068e-11 + Line search iteration 4, with norm of the rhs 39.3125 and going to 37.9809, relative residual: 2.12304e-11 + Newton system information: Norm of the rhs: 38.7369, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 120: 2.09196e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 61.7248 and going to 38.733, relative residual: 3.3334e-11 + Line search iteration 1, with norm of the rhs 78.9823 and going to 38.7342, relative residual: 4.26538e-11 + Line search iteration 2, with norm of the rhs 69.8775 and going to 38.7342, relative residual: 3.77368e-11 + Line search iteration 3, with norm of the rhs 86.1692 and going to 38.7342, relative residual: 4.6535e-11 + Newton system information: Norm of the rhs: 38.1881, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 121: 2.06232e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 68.6949 and going to 38.1846, relative residual: 3.70982e-11 + Line search iteration 1, with norm of the rhs 77.5149 and going to 38.1858, relative residual: 4.18613e-11 + Line search iteration 2, with norm of the rhs 77.577 and going to 38.1858, relative residual: 4.18948e-11 + Line search iteration 3, with norm of the rhs 39.4857 and going to 38.1858, relative residual: 2.13239e-11 + Line search iteration 4, with norm of the rhs 39.6012 and going to 38.1858, relative residual: 2.13863e-11 + Newton system information: Norm of the rhs: 42.2532, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 122: 2.28185e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 71.1122 and going to 42.2496, relative residual: 3.84036e-11 + Line search iteration 1, with norm of the rhs 87.0236 and going to 42.251, relative residual: 4.69964e-11 + Line search iteration 2, with norm of the rhs 48.2561 and going to 42.251, relative residual: 2.60603e-11 + Line search iteration 3, with norm of the rhs 45.5428 and going to 42.251, relative residual: 2.4595e-11 + Newton system information: Norm of the rhs: 40.1094, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 123: 2.16608e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 79.2321 and going to 40.1063, relative residual: 4.27887e-11 + Line search iteration 1, with norm of the rhs 45.6605 and going to 40.1076, relative residual: 2.46586e-11 + Line search iteration 2, with norm of the rhs 44.2952 and going to 40.1076, relative residual: 2.39213e-11 + Line search iteration 3, with norm of the rhs 44.1685 and going to 40.1076, relative residual: 2.38529e-11 + Line search iteration 4, with norm of the rhs 41.7235 and going to 40.1076, relative residual: 2.25324e-11 + Newton system information: Norm of the rhs: 41.3686, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 124: 2.23408e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 53.7098 and going to 41.365, relative residual: 2.90056e-11 + Line search iteration 1, with norm of the rhs 42.2716 and going to 41.3664, relative residual: 2.28284e-11 + Newton system information: Norm of the rhs: 39.8916, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 125: 2.15431e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 71.9235 and going to 39.8896, relative residual: 3.88417e-11 + Line search iteration 1, with norm of the rhs 71.0707 and going to 39.891, relative residual: 3.83812e-11 + Line search iteration 2, with norm of the rhs 65.941 and going to 39.891, relative residual: 3.56109e-11 + Newton system information: Norm of the rhs: 35.998, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 126: 1.94405e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 69.314 and going to 35.9955, relative residual: 3.74325e-11 + Line search iteration 1, with norm of the rhs 85.6074 and going to 35.9967, relative residual: 4.62316e-11 + Line search iteration 2, with norm of the rhs 75.3569 and going to 35.9967, relative residual: 4.06959e-11 + Line search iteration 3, with norm of the rhs 71.4413 and going to 35.9967, relative residual: 3.85813e-11 + Line search iteration 4, with norm of the rhs 40.2121 and going to 35.9967, relative residual: 2.17162e-11 + Newton system information: Norm of the rhs: 37.136, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 127: 2.0055e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 115.868 and going to 37.1327, relative residual: 6.25737e-11 + Line search iteration 1, with norm of the rhs 71.8526 and going to 37.1339, relative residual: 3.88034e-11 + Line search iteration 2, with norm of the rhs 77.9159 and going to 37.1339, relative residual: 4.20779e-11 + Line search iteration 3, with norm of the rhs 75.4443 and going to 37.1339, relative residual: 4.07431e-11 + Line search iteration 4, with norm of the rhs 86.3033 and going to 37.1339, relative residual: 4.66074e-11 + Newton system information: Norm of the rhs: 39.927, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 128: 2.15623e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 82.1772 and going to 39.9236, relative residual: 4.43791e-11 + Line search iteration 1, with norm of the rhs 72.3866 and going to 39.9249, relative residual: 3.90918e-11 + Line search iteration 2, with norm of the rhs 82.2322 and going to 39.9249, relative residual: 4.44088e-11 + Line search iteration 3, with norm of the rhs 53.2375 and going to 39.9249, relative residual: 2.87505e-11 + Newton system information: Norm of the rhs: 38.6037, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 129: 2.08476e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 67.0169 and going to 38.6007, relative residual: 3.61919e-11 + Line search iteration 1, with norm of the rhs 73.4455 and going to 38.602, relative residual: 3.96637e-11 + Line search iteration 2, with norm of the rhs 59.2646 and going to 38.602, relative residual: 3.20054e-11 + Newton system information: Norm of the rhs: 38.0264, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 130: 2.05359e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 68.7701 and going to 38.024, relative residual: 3.71387e-11 + Line search iteration 1, with norm of the rhs 80.8949 and going to 38.0253, relative residual: 4.36867e-11 + Line search iteration 2, with norm of the rhs 73.52 and going to 38.0253, relative residual: 3.97039e-11 + Line search iteration 3, with norm of the rhs 83.1227 and going to 38.0253, relative residual: 4.48898e-11 + Line search iteration 4, with norm of the rhs 41.5531 and going to 38.0253, relative residual: 2.24404e-11 + Newton system information: Norm of the rhs: 38.8271, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 131: 2.09683e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 87.5419 and going to 38.8239, relative residual: 4.72763e-11 + Line search iteration 1, with norm of the rhs 73.3114 and going to 38.8252, relative residual: 3.95913e-11 + Line search iteration 2, with norm of the rhs 63.2955 and going to 38.8252, relative residual: 3.41823e-11 + Line search iteration 3, with norm of the rhs 40.9057 and going to 38.8252, relative residual: 2.20908e-11 + Line search iteration 4, with norm of the rhs 39.5691 and going to 38.8252, relative residual: 2.1369e-11 + Newton system information: Norm of the rhs: 38.5323, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 132: 2.08091e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 68.3124 and going to 38.5291, relative residual: 3.68916e-11 + Line search iteration 1, with norm of the rhs 49.1458 and going to 38.5304, relative residual: 2.65408e-11 + Line search iteration 2, with norm of the rhs 39.1008 and going to 38.5304, relative residual: 2.11161e-11 + Newton system information: Norm of the rhs: 37.8596, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 133: 2.04458e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 75.2774 and going to 37.8575, relative residual: 4.0653e-11 + Line search iteration 1, with norm of the rhs 78.2847 and going to 37.8588, relative residual: 4.2277e-11 + Line search iteration 2, with norm of the rhs 73.7094 and going to 37.8588, relative residual: 3.98062e-11 + Line search iteration 3, with norm of the rhs 40.7084 and going to 37.8588, relative residual: 2.19843e-11 + Line search iteration 4, with norm of the rhs 39.3615 and going to 37.8588, relative residual: 2.12569e-11 + Newton system information: Norm of the rhs: 40.242, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 134: 2.17324e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 75.8896 and going to 40.2386, relative residual: 4.09836e-11 + Line search iteration 1, with norm of the rhs 92.3553 and going to 40.24, relative residual: 4.98758e-11 + Line search iteration 2, with norm of the rhs 72.556 and going to 40.24, relative residual: 3.91833e-11 + Line search iteration 3, with norm of the rhs 69.6229 and going to 40.24, relative residual: 3.75993e-11 + Line search iteration 4, with norm of the rhs 49.7596 and going to 40.24, relative residual: 2.68723e-11 + Newton system information: Norm of the rhs: 38.0003, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 135: 2.05217e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 90.8209 and going to 37.9972, relative residual: 4.90471e-11 + Line search iteration 1, with norm of the rhs 72.8523 and going to 37.9984, relative residual: 3.93433e-11 + Line search iteration 2, with norm of the rhs 81.1421 and going to 37.9984, relative residual: 4.38202e-11 + Line search iteration 3, with norm of the rhs 54.9077 and going to 37.9984, relative residual: 2.96525e-11 + Line search iteration 4, with norm of the rhs 38.4629 and going to 37.9984, relative residual: 2.07716e-11 + Newton system information: Norm of the rhs: 37.2794, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 136: 2.01325e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 116.928 and going to 37.2764, relative residual: 6.3146e-11 + Line search iteration 1, with norm of the rhs 80.2857 and going to 37.2776, relative residual: 4.33577e-11 + Line search iteration 2, with norm of the rhs 70.5161 and going to 37.2776, relative residual: 3.80817e-11 + Line search iteration 3, with norm of the rhs 69.4598 and going to 37.2776, relative residual: 3.75112e-11 + Line search iteration 4, with norm of the rhs 77.8235 and going to 37.2776, relative residual: 4.2028e-11 + Newton system information: Norm of the rhs: 37.2352, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 137: 2.01086e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 77.7721 and going to 37.2322, relative residual: 4.20002e-11 + Line search iteration 1, with norm of the rhs 85.673 and going to 37.2335, relative residual: 4.62671e-11 + Line search iteration 2, with norm of the rhs 72.0915 and going to 37.2335, relative residual: 3.89324e-11 + Line search iteration 3, with norm of the rhs 77.2954 and going to 37.2335, relative residual: 4.17428e-11 + Line search iteration 4, with norm of the rhs 41.3326 and going to 37.2335, relative residual: 2.23214e-11 + Newton system information: Norm of the rhs: 40.659, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 138: 2.19576e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 71.3688 and going to 40.6555, relative residual: 3.85422e-11 + Line search iteration 1, with norm of the rhs 68.9098 and going to 40.6569, relative residual: 3.72142e-11 + Newton system information: Norm of the rhs: 39.5767, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 139: 2.13731e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 85.4278 and going to 39.5756, relative residual: 4.61346e-11 + Line search iteration 1, with norm of the rhs 70.4154 and going to 39.5769, relative residual: 3.80273e-11 + Line search iteration 2, with norm of the rhs 73.7134 and going to 39.5769, relative residual: 3.98084e-11 + Line search iteration 3, with norm of the rhs 73.0193 and going to 39.5769, relative residual: 3.94335e-11 + Line search iteration 4, with norm of the rhs 66.4546 and going to 39.5769, relative residual: 3.58883e-11 + Newton system information: Norm of the rhs: 38.5438, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 140: 2.08153e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 91.4187 and going to 38.5407, relative residual: 4.937e-11 + Line search iteration 1, with norm of the rhs 86.6055 and going to 38.542, relative residual: 4.67706e-11 + Line search iteration 2, with norm of the rhs 75.9166 and going to 38.542, relative residual: 4.09982e-11 + Line search iteration 3, with norm of the rhs 43.0914 and going to 38.542, relative residual: 2.32712e-11 + Line search iteration 4, with norm of the rhs 43.8307 and going to 38.542, relative residual: 2.36704e-11 + Newton system information: Norm of the rhs: 39.8522, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 141: 2.15219e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 90.1228 and going to 39.849, relative residual: 4.86701e-11 + Line search iteration 1, with norm of the rhs 91.7147 and going to 39.8504, relative residual: 4.95298e-11 + Line search iteration 2, with norm of the rhs 74.0256 and going to 39.8504, relative residual: 3.99769e-11 + Line search iteration 3, with norm of the rhs 69.8442 and going to 39.8504, relative residual: 3.77188e-11 + Line search iteration 4, with norm of the rhs 69.0208 and going to 39.8504, relative residual: 3.72741e-11 + Newton system information: Norm of the rhs: 37.4987, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 142: 2.02509e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 67.0635 and going to 37.4958, relative residual: 3.62171e-11 + Line search iteration 1, with norm of the rhs 72.7607 and going to 37.497, relative residual: 3.92938e-11 + Line search iteration 2, with norm of the rhs 43.7419 and going to 37.497, relative residual: 2.36225e-11 + Line search iteration 3, with norm of the rhs 44.0869 and going to 37.497, relative residual: 2.38088e-11 + Line search iteration 4, with norm of the rhs 43.926 and going to 37.497, relative residual: 2.37219e-11 + Newton system information: Norm of the rhs: 42.4208, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 143: 2.2909e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 59.7291 and going to 42.4172, relative residual: 3.22562e-11 + Newton system information: Norm of the rhs: 40.8012, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 144: 2.20344e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 59.3285 and going to 40.8001, relative residual: 3.20399e-11 + Line search iteration 1, with norm of the rhs 49.0719 and going to 40.8014, relative residual: 2.65009e-11 + Line search iteration 2, with norm of the rhs 43.7013 and going to 40.8014, relative residual: 2.36005e-11 + Line search iteration 3, with norm of the rhs 43.2965 and going to 40.8014, relative residual: 2.33819e-11 + Line search iteration 4, with norm of the rhs 43.4236 and going to 40.8014, relative residual: 2.34506e-11 + Newton system information: Norm of the rhs: 40.5587, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 145: 2.19034e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 78.4526 and going to 40.5552, relative residual: 4.23677e-11 + Line search iteration 1, with norm of the rhs 78.0563 and going to 40.5566, relative residual: 4.21537e-11 + Line search iteration 2, with norm of the rhs 65.4517 and going to 40.5566, relative residual: 3.53467e-11 + Line search iteration 3, with norm of the rhs 42.8692 and going to 40.5566, relative residual: 2.31512e-11 + Line search iteration 4, with norm of the rhs 42.6239 and going to 40.5566, relative residual: 2.30187e-11 + Newton system information: Norm of the rhs: 40.9117, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 146: 2.20941e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 53.2962 and going to 40.9081, relative residual: 2.87822e-11 + Line search iteration 1, with norm of the rhs 44.4739 and going to 40.9095, relative residual: 2.40178e-11 + Line search iteration 2, with norm of the rhs 41.7716 and going to 40.9095, relative residual: 2.25584e-11 + Line search iteration 3, with norm of the rhs 41.1415 and going to 40.9095, relative residual: 2.22181e-11 + Newton system information: Norm of the rhs: 40.5898, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 147: 2.19202e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 93.1928 and going to 40.5865, relative residual: 5.03281e-11 + Line search iteration 1, with norm of the rhs 74.2232 and going to 40.5879, relative residual: 4.00836e-11 + Line search iteration 2, with norm of the rhs 52.203 and going to 40.5879, relative residual: 2.81918e-11 + Newton system information: Norm of the rhs: 40.0188, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 148: 2.16119e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 13+0 iterations. + Line search iteration 0, with norm of the rhs 86.6157 and going to 40.016, relative residual: 4.67761e-11 + Line search iteration 1, with norm of the rhs 68.7344 and going to 40.0173, relative residual: 3.71195e-11 + Line search iteration 2, with norm of the rhs 52.8158 and going to 40.0173, relative residual: 2.85228e-11 + Line search iteration 3, with norm of the rhs 40.3493 and going to 40.0173, relative residual: 2.17903e-11 + Line search iteration 4, with norm of the rhs 40.8877 and going to 40.0173, relative residual: 2.20811e-11 + Newton system information: Norm of the rhs: 43.3464, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 149: 2.34089e-11 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 12+0 iterations. + Line search iteration 0, with norm of the rhs 99.4983 and going to 43.3427, relative residual: 5.37333e-11 + Line search iteration 1, with norm of the rhs 70.9317 and going to 43.3441, relative residual: 3.83061e-11 + Line search iteration 2, with norm of the rhs 75.4052 and going to 43.3441, relative residual: 4.0722e-11 + Line search iteration 3, with norm of the rhs 75.9567 and going to 43.3441, relative residual: 4.10198e-11 + Line search iteration 4, with norm of the rhs 58.647 and going to 43.3441, relative residual: 3.16719e-11 + Newton system information: Norm of the rhs: 37.8881, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 150: 2.04612e-11 + + +WARNING: The nonlinear solver in the current timestep failed to converge. +Acting according to the parameter 'Nonlinear solver failure strategy'... +Continuing to the next timestep even though solution is not fully converged. + Postprocessing: + RMS, max velocity: 2.57e-08 m/s, 3.05e-08 m/s + Pressure min/avg/max: -0.0002411 Pa, 5e+08 Pa, 1e+09 Pa + Mass fluxes through boundary parts: 0 kg/s, 0 kg/s, -0.8139 kg/s, 0.8139 kg/s + Writing graphical output: output-nonlinear_channel_flow_tractions_benchmark/solution/solution-00000 + +Termination requested by criterion: end time + + + + +WARNING: During this computation 1 nonlinear solver failures occurred! diff --git a/tests/nonlinear_channel_flow_tractions_benchmark/statistics b/tests/nonlinear_channel_flow_tractions_benchmark/statistics new file mode 100644 index 00000000000..7eee9e34247 --- /dev/null +++ b/tests/nonlinear_channel_flow_tractions_benchmark/statistics @@ -0,0 +1,22 @@ +# 1: Time step number +# 2: Time (seconds) +# 3: Time step size (seconds) +# 4: Number of mesh cells +# 5: Number of Stokes degrees of freedom +# 6: Number of temperature degrees of freedom +# 7: Number of nonlinear iterations +# 8: Iterations for temperature solver +# 9: Iterations for Stokes solver +# 10: Velocity iterations in Stokes preconditioner +# 11: Schur complement iterations in Stokes preconditioner +# 12: RMS velocity (m/s) +# 13: Max. velocity (m/s) +# 14: Minimal pressure (Pa) +# 15: Average pressure (Pa) +# 16: Maximal pressure (Pa) +# 17: Outward mass flux through boundary with indicator 0 ("left") (kg/s) +# 18: Outward mass flux through boundary with indicator 1 ("right") (kg/s) +# 19: Outward mass flux through boundary with indicator 2 ("bottom") (kg/s) +# 20: Outward mass flux through boundary with indicator 3 ("top") (kg/s) +# 21: Visualization file name +0 0.000000000000e+00 0.000000000000e+00 256 2467 1089 150 0 3351 3660 4034 2.57340148e-08 3.05175054e-08 -2.41114358e-04 5.00000000e+08 1.00000000e+09 0.00000000e+00 0.00000000e+00 -8.13894262e-01 8.13894262e-01 output-nonlinear_channel_flow_tractions_benchmark/solution/solution-00000 diff --git a/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/screen-output b/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/screen-output index 37e13b7ac3f..5ff099536aa 100644 --- a/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/screen-output +++ b/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/screen-output @@ -1,61 +1,61 @@ Loading shared library <./libnonlinear_channel_flow_tractions_dc_picard_Stokes_GMG.debug.so> -Number of active cells: 4,096 (on 7 levels) -Number of degrees of freedom: 54,148 (33,282+4,225+16,641) +Number of active cells: 256 (on 5 levels) +Number of degrees of freedom: 3,556 (2,178+289+1,089) *** Timestep 0: t=0 seconds, dt=0 seconds Skipping temperature solve because RHS is zero. - Initial Newton Stokes residual = 9.30238e+11, v = 9.30238e+11, p = 0 + Initial Newton Stokes residual = 1.85171e+12, v = 1.85171e+12, p = 0 - Solving Stokes system (GMG)... 8+0 iterations. - Newton system information: Norm of the rhs: 9.30238e+11 + Solving Stokes system (GMG)... 13+0 iterations. + Newton system information: Norm of the rhs: 1.85171e+12 Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1 - Solving Stokes system (GMG)... 10+0 iterations. - Newton system information: Norm of the rhs: 8.35883e+10, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.0898569 + Solving Stokes system (GMG)... 12+0 iterations. + Newton system information: Norm of the rhs: 3.04031e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.16419 - Solving Stokes system (GMG)... 7+0 iterations. - Newton system information: Norm of the rhs: 1.08226e+10, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 0.0116342 + Solving Stokes system (GMG)... 14+0 iterations. + Newton system information: Norm of the rhs: 2.64066e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 0.142607 - Solving Stokes system (GMG)... 4+0 iterations. - Newton system information: Norm of the rhs: 105099, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 1.12981e-07 + Solving Stokes system (GMG)... 16+0 iterations. + Newton system information: Norm of the rhs: 2.15319e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 0.116281 - Solving Stokes system (GMG)... 8+0 iterations. - Newton system information: Norm of the rhs: 0.672766, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 7.23219e-13 + Solving Stokes system (GMG)... 16+0 iterations. + Newton system information: Norm of the rhs: 1.662e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 0.0897551 - Solving Stokes system (GMG)... 9+0 iterations. - Newton system information: Norm of the rhs: 0.0951174, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 6: 1.02251e-13 + Solving Stokes system (GMG)... 17+0 iterations. + Newton system information: Norm of the rhs: 1.22876e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 6: 0.0663582 - Solving Stokes system (GMG)... 7+0 iterations. - Newton system information: Norm of the rhs: 0.0939245, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 7: 1.00968e-13 + Solving Stokes system (GMG)... 17+0 iterations. + Newton system information: Norm of the rhs: 8.79976e+10, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 7: 0.0475224 - Solving Stokes system (GMG)... 7+0 iterations. - Newton system information: Norm of the rhs: 0.0839196, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 8: 9.0213e-14 + Solving Stokes system (GMG)... 18+0 iterations. + Newton system information: Norm of the rhs: 6.16064e+10, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 8: 0.0332701 - Solving Stokes system (GMG)... 7+0 iterations. - Newton system information: Norm of the rhs: 0.0914938, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 9: 9.83552e-14 + Solving Stokes system (GMG)... 18+0 iterations. + Newton system information: Norm of the rhs: 4.24552e+10, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 9: 0.0229276 - Solving Stokes system (GMG)... 7+0 iterations. - Newton system information: Norm of the rhs: 0.0996941, Derivative scaling factor: 0 - Relative nonlinear residual (Stokes system) after nonlinear iteration 10: 1.0717e-13 + Solving Stokes system (GMG)... 18+0 iterations. + Newton system information: Norm of the rhs: 2.89428e+10, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 10: 0.0156303 WARNING: The nonlinear solver in the current timestep failed to converge. Acting according to the parameter 'Nonlinear solver failure strategy'... Continuing to the next timestep even though solution is not fully converged. Postprocessing: - RMS, max velocity: 1.14e-07 m/s, 1.56e-07 m/s - Pressure min/avg/max: -2.689e-06 Pa, 5e+08 Pa, 1e+09 Pa - Mass fluxes through boundary parts: 0 kg/s, 0 kg/s, -3.473 kg/s, 3.473 kg/s + RMS, max velocity: 1.96e-08 m/s, 2.33e-08 m/s + Pressure min/avg/max: -2.116 Pa, 5e+08 Pa, 1e+09 Pa + Mass fluxes through boundary parts: 0 kg/s, 0 kg/s, -0.6176 kg/s, 0.6176 kg/s Writing graphical output: output-nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/solution/solution-00000 Termination requested by criterion: end time diff --git a/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/statistics b/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/statistics index fa9992a8cab..6a48c1b286c 100644 --- a/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/statistics +++ b/tests/nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/statistics @@ -19,4 +19,4 @@ # 19: Outward mass flux through boundary with indicator 2 ("bottom") (kg/s) # 20: Outward mass flux through boundary with indicator 3 ("top") (kg/s) # 21: Visualization file name -0 0.000000000000e+00 0.000000000000e+00 4096 37507 16641 10 0 64 84 84 1.14108866e-07 1.56248062e-07 -2.68875150e-06 5.00000000e+08 1.00000000e+09 0.00000000e+00 0.00000000e+00 -3.47275156e+00 3.47275156e+00 output-nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/solution/solution-00000 +0 0.000000000000e+00 0.000000000000e+00 256 2467 1089 10 0 149 169 169 1.95647143e-08 2.32828598e-08 -2.11569031e+00 5.00000000e+08 1.00000000e+09 0.00000000e+00 0.00000000e+00 -6.17622785e-01 6.17622785e-01 output-nonlinear_channel_flow_tractions_dc_picard_Stokes_GMG/solution/solution-00000 diff --git a/tests/nonlinear_channel_flow_velocities_benchmark.cc b/tests/nonlinear_channel_flow_velocities_benchmark.cc new file mode 100644 index 00000000000..b0c1b8e736e --- /dev/null +++ b/tests/nonlinear_channel_flow_velocities_benchmark.cc @@ -0,0 +1 @@ +#include "../benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/simple_nonlinear.cc" diff --git a/tests/nonlinear_channel_flow_velocities_benchmark.prm b/tests/nonlinear_channel_flow_velocities_benchmark.prm new file mode 100644 index 00000000000..faf912fb9aa --- /dev/null +++ b/tests/nonlinear_channel_flow_velocities_benchmark.prm @@ -0,0 +1,4 @@ +# Test the benchmark parameter file +# benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_v.prm + +include $ASPECT_SOURCE_DIR/benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow/input_v.prm diff --git a/tests/nonlinear_channel_flow_velocities_benchmark/screen-output b/tests/nonlinear_channel_flow_velocities_benchmark/screen-output new file mode 100644 index 00000000000..3cede1e78b1 --- /dev/null +++ b/tests/nonlinear_channel_flow_velocities_benchmark/screen-output @@ -0,0 +1,240 @@ + +Loading shared library <./libnonlinear_channel_flow_velocities_benchmark.debug.so> + +Number of active cells: 256 (on 5 levels) +Number of degrees of freedom: 3,556 (2,178+289+1,089) + +*** Timestep 0: t=0 seconds, dt=0 seconds + Skipping temperature solve because RHS is zero. + Initial Newton Stokes residual = 1.06033e+18, v = 7.28522e+17, p = 7.70426e+17 + + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 33+0 iterations. + Newton system information: Norm of the rhs: 1.06033e+18 + Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1 + + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 45+0 iterations. + Newton system information: Norm of the rhs: 2.19191e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 2.06719e-07 + + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 48+0 iterations. + Newton system information: Norm of the rhs: 1.47477e+11, Derivative scaling factor: 0 + Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 1.39086e-07 + + Switching from defect correction form of Picard to the Newton solver scheme. + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 26+0 iterations. + Line search iteration 0, with norm of the rhs 3.72586e+11 and going to 1.47463e+11, relative residual: 3.51387e-07 + Line search iteration 1, with norm of the rhs 2.88847e+11 and going to 1.47468e+11, relative residual: 2.72412e-07 + Newton system information: Norm of the rhs: 8.02499e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 7.56839e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 26+0 iterations. + Line search iteration 0, with norm of the rhs 1.84411e+11 and going to 8.02425e+10, relative residual: 1.73918e-07 + Line search iteration 1, with norm of the rhs 1.15262e+11 and going to 8.02452e+10, relative residual: 1.08704e-07 + Newton system information: Norm of the rhs: 4.85681e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 4.58046e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 23+0 iterations. + Line search iteration 0, with norm of the rhs 9.09999e+10 and going to 4.85634e+10, relative residual: 8.58222e-08 + Line search iteration 1, with norm of the rhs 6.72566e+10 and going to 4.85651e+10, relative residual: 6.34298e-08 + Newton system information: Norm of the rhs: 2.78022e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 6: 2.62203e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 25+0 iterations. + Line search iteration 0, with norm of the rhs 6.1696e+10 and going to 2.77995e+10, relative residual: 5.81856e-08 + Line search iteration 1, with norm of the rhs 5.20151e+10 and going to 2.78004e+10, relative residual: 4.90555e-08 + Newton system information: Norm of the rhs: 1.70634e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 7: 1.60925e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 32+0 iterations. + Line search iteration 0, with norm of the rhs 3.57069e+10 and going to 1.70617e+10, relative residual: 3.36753e-08 + Line search iteration 1, with norm of the rhs 3.3656e+10 and going to 1.70623e+10, relative residual: 3.17411e-08 + Line search iteration 2, with norm of the rhs 1.99766e+10 and going to 1.70623e+10, relative residual: 1.884e-08 + Newton system information: Norm of the rhs: 1.38215e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 8: 1.30351e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 38+0 iterations. + Line search iteration 0, with norm of the rhs 2.31006e+10 and going to 1.38201e+10, relative residual: 2.17862e-08 + Line search iteration 1, with norm of the rhs 2.38371e+10 and going to 1.38206e+10, relative residual: 2.24808e-08 + Line search iteration 2, with norm of the rhs 1.89865e+10 and going to 1.38206e+10, relative residual: 1.79062e-08 + Newton system information: Norm of the rhs: 1.26085e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 9: 1.18911e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 47+0 iterations. + Line search iteration 0, with norm of the rhs 1.28735e+10 and going to 1.26072e+10, relative residual: 1.2141e-08 + Line search iteration 1, with norm of the rhs 1.56504e+10 and going to 1.26076e+10, relative residual: 1.47599e-08 + Line search iteration 2, with norm of the rhs 1.32102e+10 and going to 1.26076e+10, relative residual: 1.24585e-08 + Newton system information: Norm of the rhs: 1.12354e+10, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 10: 1.05961e-08 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 46+0 iterations. + Newton system information: Norm of the rhs: 8.52691e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 11: 8.04174e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 38+0 iterations. + Line search iteration 0, with norm of the rhs 1.07576e+10 and going to 8.52605e+09, relative residual: 1.01455e-08 + Newton system information: Norm of the rhs: 8.31518e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 12: 7.84207e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 26+0 iterations. + Line search iteration 0, with norm of the rhs 8.88087e+09 and going to 8.31435e+09, relative residual: 8.37557e-09 + Newton system information: Norm of the rhs: 6.49023e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 13: 6.12095e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 21+0 iterations. + Line search iteration 0, with norm of the rhs 7.8591e+09 and going to 6.48958e+09, relative residual: 7.41194e-09 + Newton system information: Norm of the rhs: 6.15863e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 14: 5.80821e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 26+0 iterations. + Line search iteration 0, with norm of the rhs 7.30912e+09 and going to 6.15801e+09, relative residual: 6.89325e-09 + Newton system information: Norm of the rhs: 5.71093e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 15: 5.38599e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 27+0 iterations. + Line search iteration 0, with norm of the rhs 6.92805e+09 and going to 5.71036e+09, relative residual: 6.53386e-09 + Newton system information: Norm of the rhs: 5.35661e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 16: 5.05183e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 24+0 iterations. + Line search iteration 0, with norm of the rhs 6.38186e+09 and going to 5.35607e+09, relative residual: 6.01875e-09 + Newton system information: Norm of the rhs: 5.00374e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 17: 4.71903e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 27+0 iterations. + Line search iteration 0, with norm of the rhs 6.10245e+09 and going to 5.00323e+09, relative residual: 5.75524e-09 + Newton system information: Norm of the rhs: 4.72526e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 18: 4.4564e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 20+0 iterations. + Line search iteration 0, with norm of the rhs 5.61402e+09 and going to 4.72479e+09, relative residual: 5.29459e-09 + Newton system information: Norm of the rhs: 4.33869e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 19: 4.09183e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 25+0 iterations. + Line search iteration 0, with norm of the rhs 4.88694e+09 and going to 4.33826e+09, relative residual: 4.60888e-09 + Newton system information: Norm of the rhs: 3.76733e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 20: 3.55298e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 22+0 iterations. + Line search iteration 0, with norm of the rhs 4.54727e+09 and going to 3.76695e+09, relative residual: 4.28854e-09 + Newton system information: Norm of the rhs: 3.43811e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 21: 3.24249e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 19+0 iterations. + Line search iteration 0, with norm of the rhs 3.54679e+09 and going to 3.43777e+09, relative residual: 3.34498e-09 + Newton system information: Norm of the rhs: 2.78833e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 22: 2.62968e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 28+0 iterations. + Line search iteration 0, with norm of the rhs 3.27037e+09 and going to 2.78805e+09, relative residual: 3.08429e-09 + Newton system information: Norm of the rhs: 2.12249e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 23: 2.00173e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 30+0 iterations. + Newton system information: Norm of the rhs: 1.44888e+09, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 24: 1.36644e-09 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 24+0 iterations. + Newton system information: Norm of the rhs: 1.5275e+07, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 25: 1.44059e-11 + + The linear solver tolerance is set to 0.00112825. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 144+0 iterations. + Newton system information: Norm of the rhs: 3.44205e+06, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 26: 3.24621e-12 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 29+0 iterations. + Newton system information: Norm of the rhs: 920411, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 27: 8.68042e-13 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 40+0 iterations. + Newton system information: Norm of the rhs: 282220, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 28: 2.66162e-13 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 51+0 iterations. + Newton system information: Norm of the rhs: 105708, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 29: 9.96939e-14 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 52+0 iterations. + Newton system information: Norm of the rhs: 41579, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 30: 3.92132e-14 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 48+0 iterations. + Newton system information: Norm of the rhs: 16988.7, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 31: 1.60221e-14 + + The linear solver tolerance is set to 0.01. Stabilization Preconditioner is SPD and A block is SPD. + Rebuilding Stokes preconditioner... + Solving Stokes system (AMG)... 45+0 iterations. + Newton system information: Norm of the rhs: 6964.7, Derivative scaling factor: 1 + Relative nonlinear residual (Stokes system) after nonlinear iteration 32: 6.56842e-15 + + + Postprocessing: + RMS, max velocity: 2.57e-08 m/s, 3.05e-08 m/s + Pressure min/avg/max: -1.074e+07 Pa, 5.008e+08 Pa, 1.012e+09 Pa + Mass fluxes through boundary parts: 0 kg/s, 0 kg/s, -0.8139 kg/s, 0.8139 kg/s + Writing graphical output: output-nonlinear_channel_flow_velocities_benchmark/solution/solution-00000 + +Termination requested by criterion: end time + + + diff --git a/tests/nonlinear_channel_flow_velocities_benchmark/statistics b/tests/nonlinear_channel_flow_velocities_benchmark/statistics new file mode 100644 index 00000000000..a85071fa8a4 --- /dev/null +++ b/tests/nonlinear_channel_flow_velocities_benchmark/statistics @@ -0,0 +1,22 @@ +# 1: Time step number +# 2: Time (seconds) +# 3: Time step size (seconds) +# 4: Number of mesh cells +# 5: Number of Stokes degrees of freedom +# 6: Number of temperature degrees of freedom +# 7: Number of nonlinear iterations +# 8: Iterations for temperature solver +# 9: Iterations for Stokes solver +# 10: Velocity iterations in Stokes preconditioner +# 11: Schur complement iterations in Stokes preconditioner +# 12: RMS velocity (m/s) +# 13: Max. velocity (m/s) +# 14: Minimal pressure (Pa) +# 15: Average pressure (Pa) +# 16: Maximal pressure (Pa) +# 17: Outward mass flux through boundary with indicator 0 ("left") (kg/s) +# 18: Outward mass flux through boundary with indicator 1 ("right") (kg/s) +# 19: Outward mass flux through boundary with indicator 2 ("bottom") (kg/s) +# 20: Outward mass flux through boundary with indicator 3 ("top") (kg/s) +# 21: Visualization file name +0 0.000000000000e+00 0.000000000000e+00 256 2467 1089 32 0 1123 1188 2501 2.57348854e-08 3.05182567e-08 -1.07436330e+07 5.00753885e+08 1.01225140e+09 0.00000000e+00 0.00000000e+00 -8.13924078e-01 8.13924078e-01 output-nonlinear_channel_flow_velocities_benchmark/solution/solution-00000