[cig-commits] r19179 - short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8

[brad at geodynamics.org]

Author: brad
Date: 2011-11-10 14:05:34 -0800 (Thu, 10 Nov 2011)
New Revision: 19179

Modified:
   short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg
   short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step12.cfg
   short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step13.cfg
   short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step14.cfg
Log:
Update of solver tolerances for friction examples.

Modified: short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg	2011-11-10 21:31:08 UTC (rev 19178)
+++ short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg	2011-11-10 22:05:34 UTC (rev 19179)
@@ -78,7 +78,7 @@
 sub_pc_factor_shift_type = nonzero
 
 # Convergence parameters.
-ksp_rtol = 1.0e-8
+ksp_rtol = 1.0e-10
 ksp_atol = 1.0e-12
 ksp_max_it = 100
 ksp_gmres_restart = 50
@@ -89,10 +89,11 @@
 ksp_converged_reason = true
 
 # Nonlinear solver monitoring options.
-snes_rtol = 1.0e-8
-snes_atol = 1.0e-18
+snes_rtol = 1.0e-10
+snes_atol = 1.0e-9
 snes_max_it = 100
 snes_monitor = true
+snes_ls_monitor = true
 snes_view = true
 snes_converged_reason = true
 

Modified: short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step12.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step12.cfg	2011-11-10 21:31:08 UTC (rev 19178)
+++ short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step12.cfg	2011-11-10 22:05:34 UTC (rev 19179)
@@ -88,7 +88,7 @@
 db_initial = spatialdata.spatialdb.UniformDB
 db_initial.label = Dirichlet BC on +x
 db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [-1.0*m,0.0*m]
+db_initial.data = [-0.5*m,0.0*m]
 
 db_rate = spatialdata.spatialdb.UniformDB
 db_rate.label = Dirichlet rate BC on +x
@@ -156,10 +156,6 @@
 #friction_ksp_view = true
 friction_ksp_converged_reason = true
 
-# Reduce convergence tolerances.
-ksp_rtol = 1.0e-12
-ksp_atol = 1.0e-15
-
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------

Modified: short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step13.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step13.cfg	2011-11-10 21:31:08 UTC (rev 19178)
+++ short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step13.cfg	2011-11-10 22:05:34 UTC (rev 19179)
@@ -136,7 +136,7 @@
 friction.db_properties = spatialdata.spatialdb.UniformDB
 friction.db_properties.label = Slip weakening
 friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
-friction.db_properties.data = [0.6,0.5,0.2*m,0.0*Pa]
+friction.db_properties.data = [0.6,0.4,0.2*m,0.0*Pa]
 
 # ----------------------------------------------------------------------
 # PETSc settings
@@ -156,10 +156,6 @@
 #friction_ksp_view = true
 friction_ksp_converged_reason = true
 
-# Reduce convergence tolerances.
-ksp_rtol = 1.0e-12
-ksp_atol = 1.0e-15
-
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------

Modified: short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step14.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step14.cfg	2011-11-10 21:31:08 UTC (rev 19178)
+++ short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/step14.cfg	2011-11-10 22:05:34 UTC (rev 19179)
@@ -118,6 +118,7 @@
 [pylithapp.timedependent.interfaces.fault]
 # The label corresponds to the name of the nodeset in CUBIT.
 label = fault
+zero_tolerance = 1.0e-12
 
 # Use the rate-and-state aging friction model.
 friction = pylith.friction.RateStateAgeing