[cig-commits] r20065 - short/3D/PyLith/benchmarks/trunk/static/performance/solvertest

[knepley at geodynamics.org]

Author: knepley
Date: 2012-05-10 12:21:33 -0700 (Thu, 10 May 2012)
New Revision: 20065

Modified:
   short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/fieldsplit_mult.cfg
   short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_full.cfg
Log:
More optimization of Schur parameters

Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/fieldsplit_mult.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/fieldsplit_mult.cfg	2012-05-10 18:19:04 UTC (rev 20064)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/fieldsplit_mult.cfg	2012-05-10 19:21:33 UTC (rev 20065)
@@ -11,3 +11,11 @@
 fs_fieldsplit_1_pc_type = jacobi
 fs_fieldsplit_0_ksp_type = preonly
 fs_fieldsplit_1_ksp_type = preonly
+
+# Need to put in near null vectors
+#pc_type gamg 
+#pc_gamg_type agg 
+#pc_gamg_agg_nsmooths 1
+#pc_gamg_verbose 2
+#pc_gamg_threshold .05
+#pc_gamg_coarse_eq_limit 50

Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_full.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_full.cfg	2012-05-10 18:19:04 UTC (rev 20064)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_full.cfg	2012-05-10 19:21:33 UTC (rev 20065)
@@ -6,7 +6,7 @@
 [pylithapp.petsc]
 fs_pc_type = fieldsplit
 fs_pc_fieldsplit_type = schur
-fs_pc_fieldsplit_schur_factorization_type = full
+fs_pc_fieldsplit_real_diagonal = true
 fs_pc_fieldsplit_schur_precondition = user
 
 fs_fieldsplit_1_ksp_monitor = true
@@ -17,6 +17,7 @@
 # 39 iterations fs_fieldsplit_1_ solve.
 # 1 iteration fs_fieldsplit_0_ solve.
 
+#fs_pc_fieldsplit_schur_factorization_type = full
 #fs_fieldsplit_0_pc_type = lu
 #fs_fieldsplit_0_ksp_type = gmres
 #fs_fieldsplit_0_ksp_rtol = 1.0e-10
@@ -25,13 +26,18 @@
 #fs_fieldsplit_1_ksp_rtol = 1.0e-10
 
 # ----------------------------------------------------------------------
-# Tuned parameters (13 iterations)
-# 22 iterations fs_fieldsplit_1_ solve
+# Tuned parameters (8 iterations) (15 iterations with upper factorization)
+# 18-20 iterations fs_fieldsplit_1_ solve
 # 0 iterations fs_fieldsplit_0_ solve (preonly)
 
+fs_pc_fieldsplit_schur_factorization_type = full
 fs_fieldsplit_0_pc_type = ml
+fs_fieldsplit_0_pc_mg_type = full
+fs_fieldsplit_0_pc_mg_smoothdown = 3
+fs_fieldsplit_0_pc_mg_smoothup   = 2
+fs_fieldsplit_0_mg_levels_3_ksp_type = richardson
+fs_fieldsplit_0_mg_levels_3_pc_type  = sor
 fs_fieldsplit_0_ksp_type = preonly
-fs_fieldsplit_0_ksp_rtol = 1.0e-6
 fs_fieldsplit_1_pc_type = jacobi
 fs_fieldsplit_1_ksp_type = gmres
-fs_fieldsplit_1_ksp_rtol = 1.0e-6
+fs_fieldsplit_1_ksp_rtol = 1.0e-5