From c85dd7be78f3e83912adc0a99ed4e7cf8e749df7 Mon Sep 17 00:00:00 2001 From: Rene Gassmoeller Date: Fri, 22 Nov 2024 13:25:47 +0100 Subject: [PATCH] Simplify Spiegelman run script --- .../spiegelman_et_al_2016/bash.sh | 73 ++++++------ .../spiegelman_et_al_2016/plot.gnuplot | 108 +++++++++--------- 2 files changed, 89 insertions(+), 92 deletions(-) diff --git a/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/bash.sh b/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/bash.sh index 2996b8da4dc..1382f607b62 100755 --- a/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/bash.sh +++ b/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/bash.sh @@ -28,48 +28,45 @@ for i_BV in "1e23" "1e24" "5e24"; do # Background viscosity # Now run Newton solver models with varying solver paramters - for i_NSP in "SPD" "none"; do # Newton stabilization preconditioner - for i_NSA in "SPD" "none"; do # Newton stabilization A block of matrix - for i_UFS in "false"; do # Use Newton failsafe - for i_P in 0 5 15 25 150; do # Minimum Picard iterations, note 150 represents pure defect-correction Picard - for i_MLT in "9e-1" "1e-8"; do # Maximum linear tolerance - for i_RSM in "false"; do # Residual scaling method + for i_stabilization in "SPD" "none"; do # Newton stabilization + for i_UFS in "false"; do # Use Newton failsafe + for i_P in 0 5 15 25 150; do # Minimum Picard iterations, note 150 represents pure defect-correction Picard + for i_MLT in "9e-1" "1e-8"; do # Maximum linear tolerance + for i_RSM in "false"; do # Residual scaling method + U="7.92219116e-11" + if [ $i_vel == 25 ]; then 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 + elif [ $i_vel == 50 ]; then + U="1.58443823e-10" + elif [ $i_vel == 125 ]; then + U="3.96109558e-10" + fi - dirname_base="singleAdvectionandNewtonStokes""_UFS_${i_UFS}_NSP_${i_NSP}_NSA_${i_NSA}_mLT_${i_MLT}_P${i_P}_RSM_${i_RSM}_vel_${i_vel}_BV_${i_BV}" - dirname="results/$dirname_base" - infilename="${dirname}/input.prm" - outfilename="${dirname}/output.log" - errorfilename="${dirname}/error.log" - outplotfilename="${dirname}/plot.dat" - - echo "Starting $dirname" - mkdir -p $dirname - cp "$infile" "$dirname/$infile" + dirname_base="singleAdvectionandNewtonStokes""_UFS_${i_UFS}_NS_${i_stabilization}_mLT_${i_MLT}_P${i_P}_RSM_${i_RSM}_vel_${i_vel}_BV_${i_BV}" + dirname="results/$dirname_base" + infilename="${dirname}/input.prm" + outfilename="${dirname}/output.log" + errorfilename="${dirname}/error.log" + outplotfilename="${dirname}/plot.dat" + + echo "Starting $dirname" + mkdir -p $dirname + cp "$infile" "$dirname/$infile" - sed \ - -e "s/set Function expression = if(x<60e3,.*/ set Function expression = if(x<60e3,$U,-$U);0/g" \ - -e "s/set Reference viscosity .*/ set Reference viscosity = $i_BV/g" \ - -e "s/set Output directory .*/set Output directory = results\/$dirname_base/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 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.prm > "$infilename" + sed \ + -e "s/set Function expression = if(x<60e3,.*/ set Function expression = if(x<60e3,$U,-$U);0/g" \ + -e "s/set Reference viscosity .*/ set Reference viscosity = $i_BV/g" \ + -e "s/set Output directory .*/set Output directory = results\/$dirname_base/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_stabilization/g" \ + -e "s/set Stabilization velocity block .*/set Stabilization velocity block = $i_stabilization/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.prm > "$infilename" - nohup mpirun -np $processes ./aspect-release $infilename > $outfilename 2>$errorfilename - grep "Relative nonlinear residual " $outfilename > $outplotfilename - - done + nohup mpirun -np $processes ./aspect-release $infilename > $outfilename 2>$errorfilename + grep "Relative nonlinear residual " $outfilename > $outplotfilename done done done diff --git a/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/plot.gnuplot b/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/plot.gnuplot index cd5bda6480d..c321d786814 100644 --- a/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/plot.gnuplot +++ b/benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/plot.gnuplot @@ -11,15 +11,15 @@ set key set title "Vel: 2.5cm/yr, \\eta_{ref}: 1e23, mLT: 9e-1" plot \ 'results/singleAdvectionanditeratedStokes_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'purple' t 'Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P150_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P150_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' unset key @@ -27,44 +27,44 @@ set origin 0.33,0.48 set title "Vel: 5cm/yr, \\eta_{ref}: 1e24, mLT: 9e-1" plot \ 'results/singleAdvectionanditeratedStokes_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'purple' t 'Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P150_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P150_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' set origin 0.66,0.48 set title "Vel: 12.5cm/yr, \\eta_{ref}: 5e24, mLT: 9e-1" plot \ 'results/singleAdvectionanditeratedStokes_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'purple' t 'Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P150_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_9e-1_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_9e-1_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P150_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_9e-1_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_9e-1_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' set origin 0.0,0.0 set title "Vel: 2.5cm/yr, \\eta_{ref}: 1e23, mLT: 1e-8" plot \ 'results/singleAdvectionanditeratedStokes_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'purple' t 'Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P150_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P150_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P0_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P5_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P15_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P25_RSM_false_vel_25_BV_1e23/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' unset key @@ -72,26 +72,26 @@ set origin 0.33,0.0 set title "Vel: 5cm/yr, \\eta_{ref}: 1e24, mLT: 1e-8" plot \ 'results/singleAdvectionanditeratedStokes_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'purple' t 'Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P150_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P150_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P0_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P5_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P15_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P25_RSM_false_vel_50_BV_1e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' set origin 0.66,0.0 set title "Vel: 12.5cm/yr, \\eta_{ref}: 5e24, mLT: 1e-8" plot \ 'results/singleAdvectionanditeratedStokes_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'purple' t 'Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P150_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_none_NSA_none_mLT_1e-8_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ -'results/singleAdvectionandNewtonStokes_UFS_false_NSP_SPD_NSA_SPD_mLT_1e-8_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied' +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P150_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'purple' t 'DC Picard', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'black' t '0 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P0_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'black' t '0 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'blue' t '5 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P5_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'blue' t '5 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'orange' t '15 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P15_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'orange' t '15 Picard stabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_none_mLT_1e-8_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w l lw 4 lc rgb 'red' t '25 Picard unstabilizied', \ +'results/singleAdvectionandNewtonStokes_UFS_false_NS_SPD_mLT_1e-8_P25_RSM_false_vel_125_BV_5e24/plot.dat' u 9:10 w p pt 1 lw 4 lc rgb 'red' t '25 Picard stabilizied'