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

[brad at geodynamics.org]

Author: brad
Date: 2012-05-09 10:25:43 -0700 (Wed, 09 May 2012)
New Revision: 20058

Modified:
   short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_full.cfg
   short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_lower.cfg
   short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_upper.cfg
Log:
Updated schur settings.

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-09 01:19:45 UTC (rev 20057)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_full.cfg	2012-05-09 17:25:43 UTC (rev 20058)
@@ -1,22 +1,37 @@
 [pylithapp.timedependent.formulation]
 split_fields = True
 matrix_type = aij
+use_custom_constraint_pc = True
 
 [pylithapp.petsc]
 fs_pc_type = fieldsplit
 fs_pc_fieldsplit_type = schur
 fs_pc_fieldsplit_schur_factorization_type = full
-# diag -> user, turn on custom pc
-fs_pc_fieldsplit_schur_precondition = diag
-#fieldsplit_0_pc_type = ml
-fs_fieldsplit_0_pc_type = lu
-# fieldsplit_0_ksp_type = preonly
-# fieldsplit_1_ksp_type = minres
-# slowly backoff tolerance, try minres
-fs_fieldsplit_0_ksp_type = gmres
-fs_fieldsplit_0_ksp_rtol = 1.0e-10
-# jacobi or none (ilu or ml)
+fs_pc_fieldsplit_schur_precondition = user
+
+fs_fieldsplit_1_ksp_monitor = true
+fs_fieldsplit_0_ksp_monitor = true
+
+# ----------------------------------------------------------------------
+# Original settings (1 iteration)
+# 39 iterations fs_fieldsplit_1_ solve.
+# 1 iteration fs_fieldsplit_0_ solve.
+
+#fs_fieldsplit_0_pc_type = lu
+#fs_fieldsplit_0_ksp_type = gmres
+#fs_fieldsplit_0_ksp_rtol = 1.0e-10
+#fs_fieldsplit_1_pc_type = jacobi
+#fs_fieldsplit_1_ksp_type = gmres
+#fs_fieldsplit_1_ksp_rtol = 1.0e-10
+
+# ----------------------------------------------------------------------
+# Tuned parameters (13 iterations)
+# 22 iterations fs_fieldsplit_1_ solve
+# 0 iterations fs_fieldsplit_0_ solve (preonly)
+
+fs_fieldsplit_0_pc_type = ml
+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-10
-fs_fieldsplit_1_ksp_monitor = true
\ No newline at end of file
+fs_fieldsplit_1_ksp_rtol = 1.0e-6

Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_lower.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_lower.cfg	2012-05-09 01:19:45 UTC (rev 20057)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_lower.cfg	2012-05-09 17:25:43 UTC (rev 20058)
@@ -1,18 +1,25 @@
 [pylithapp.timedependent.formulation]
-split_fields = False
+split_fields = True
 matrix_type = aij
+use_custom_constraint_pc = True
 
 [pylithapp.petsc]
-ksp_gmres_restart = 100
-pc_type = fieldsplit
-pc_fieldsplit_detect_saddle_point = True
-pc_fieldsplit_type = schur
-pc_fieldsplit_schur_factorization_type = lower
-pc_fieldsplit_schur_precondition = diag
-fieldsplit_0_pc_type = ml
-fieldsplit_0_ksp_type = preonly
-fieldsplit_1_pc_type = none
-#fieldsplit_1_ksp_type = minres
-fieldsplit_1_ksp_type = gmres
-fieldsplit_1_ksp_rtol = 1e-10
+fs_pc_type = fieldsplit
+fs_pc_fieldsplit_type = schur
+fs_pc_fieldsplit_schur_factorization_type = lower
+fs_pc_fieldsplit_schur_precondition = user
 
+fs_fieldsplit_1_ksp_monitor = true
+fs_fieldsplit_0_ksp_monitor = true
+
+# ----------------------------------------------------------------------
+# Tuned parameters (24 iterations)
+# 22 iterations fs_fieldsplit_1_ solve
+# 0 iterations fs_fieldsplit_0_ solve (preonly)
+
+fs_fieldsplit_0_pc_type = ml
+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 = 1e-6

Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_upper.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_upper.cfg	2012-05-09 01:19:45 UTC (rev 20057)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/schur_upper.cfg	2012-05-09 17:25:43 UTC (rev 20058)
@@ -1,17 +1,25 @@
 [pylithapp.timedependent.formulation]
-split_fields = False
+split_fields = True
 matrix_type = aij
+use_custom_constraint_pc = True
 
 [pylithapp.petsc]
-ksp_gmres_restart = 100
-pc_type = fieldsplit
-pc_fieldsplit_detect_saddle_point = True
-pc_fieldsplit_type = schur
-pc_fieldsplit_schur_factorization_type = upper
-pc_fieldsplit_schur_precondition = diag
-fieldsplit_0_pc_type = ml
-fieldsplit_0_ksp_type = preonly
-fieldsplit_1_pc_type = none
-#fieldsplit_1_ksp_type = minres
-fieldsplit_1_ksp_type = gmres
-fieldsplit_1_ksp_rtol = 1e-10
+fs_pc_type = fieldsplit
+fs_pc_fieldsplit_type = schur
+fs_pc_fieldsplit_schur_factorization_type = upper
+fs_pc_fieldsplit_schur_precondition = user
+
+fs_fieldsplit_1_ksp_monitor = true
+fs_fieldsplit_0_ksp_monitor = true
+
+# ----------------------------------------------------------------------
+# Tuned parameters (24 iterations)
+# 22 iterations fs_fieldsplit_1_ solve
+# 0 iterations fs_fieldsplit_0_ solve (preonly)
+
+fs_fieldsplit_0_pc_type = ml
+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 = 1e-6