[cig-commits] [commit] knepley/fix-faults-parallel: Finish full-scale test for spontaneous rupture with embedded fault. (cdddc30)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/pylith

On branch  : knepley/fix-faults-parallel
Link       : https://github.com/geodynamics/pylith/compare/0f8af6e7a83ed15f72a991e2a2f8b66fd29f9c12...cdddc3055566e8ef7ff029fc6fab3ccb89daa25b

>---------------------------------------------------------------

commit cdddc3055566e8ef7ff029fc6fab3ccb89daa25b
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Wed May 28 10:18:43 2014 -0700

    Finish full-scale test for spontaneous rupture with embedded fault.


>---------------------------------------------------------------

cdddc3055566e8ef7ff029fc6fab3ccb89daa25b
 tests_auto/2d/tri3/TestShearDispFriction.py |  6 ++--
 tests_auto/2d/tri3/sheardispfriction.cfg    | 43 +++++++++++++++--------------
 2 files changed, 26 insertions(+), 23 deletions(-)

diff --git a/tests_auto/2d/tri3/TestShearDispFriction.py b/tests_auto/2d/tri3/TestShearDispFriction.py
index 96b08cd..5085c85 100644
--- a/tests_auto/2d/tri3/TestShearDispFriction.py
+++ b/tests_auto/2d/tri3/TestShearDispFriction.py
@@ -134,10 +134,10 @@ class TestShearDispFriction(TestTri3):
     Calculate fault info.
     """
 
-    strikeDir = (1.0, 0.0)
+    strikeDir = (-1.0, 0.0)
     normalDir = (0.0, 1.0)
-    initialTraction = (0.0, -1.0e+7)
-    frictionCoefficient = 0.6
+    initialTraction = (0.0, -100.0e+6)
+    frictionCoefficient = 1.0
 
     nvertices = self.faultMesh['nvertices']
 
diff --git a/tests_auto/2d/tri3/sheardispfriction.cfg b/tests_auto/2d/tri3/sheardispfriction.cfg
index c5e5775..abda561 100644
--- a/tests_auto/2d/tri3/sheardispfriction.cfg
+++ b/tests_auto/2d/tri3/sheardispfriction.cfg
@@ -1,5 +1,5 @@
 [sheardispfriction]
-nodes = 1
+nodes = 3
 
 [sheardispfriction.launcher] # WARNING: THIS IS NOT PORTABLE
 command = mpirun -np ${nodes}
@@ -8,17 +8,17 @@ command = mpirun -np ${nodes}
 # journal
 # ----------------------------------------------------------------------
 [sheardispfriction.journal.info]
-sheardispfriction = 1
-timedependent = 1
-implicit = 1
-petsc = 1
-solvernonlinear = 1
-meshimporter = 1
-meshiocubit = 1
-implicitelasticity = 1
-quadrature2d = 1
-faultcohesivekin = 1
-fiatsimplex = 1
+#sheardispfriction = 1
+#timedependent = 1
+#implicit = 1
+#petsc = 1
+#solvernonlinear = 1
+#meshimporter = 1
+#meshiocubit = 1
+#implicitelasticity = 1
+#quadrature2d = 1
+#faultcohesivekin = 1
+#fiatsimplex = 1
 
 # ----------------------------------------------------------------------
 # mesh_generator
@@ -37,6 +37,7 @@ coordsys.space_dim = 2
 [sheardispfriction.timedependent]
 dimension = 2
 normalizer.length_scale = 1.0*km
+normalizer.relaxation_time = 1.0*s
 formulation.solver = pylith.problems.SolverNonlinear
 
 [sheardispfriction.timedependent.formulation.time_step]
@@ -109,7 +110,8 @@ friction.label = Static friction
 friction.db_properties = spatialdata.spatialdb.UniformDB
 friction.db_properties.label = Static friction
 friction.db_properties.values = [friction-coefficient,cohesion]
-friction.db_properties.data = [0.6,0.0*Pa]
+friction.db_properties.data = [1.0,0.0*Pa]
+zero_tolerance = 1.0e-10
 
 traction_perturbation = pylith.faults.TractPerturbation
 
@@ -117,7 +119,7 @@ traction_perturbation = pylith.faults.TractPerturbation
 db_initial = spatialdata.spatialdb.UniformDB
 db_initial.label = Initial fault tractions
 db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [0.0*Pa, -10.0*MPa]
+db_initial.data = [0.0*Pa, -100.0*MPa]
 
 # ----------------------------------------------------------------------
 # PETSc
@@ -129,21 +131,21 @@ pc_type = asm
 sub_pc_factor_shift_type = none
 
 ksp_rtol = 1.0e-20
-ksp_atol = 1.0e-10
+ksp_atol = 1.0e-11
 ksp_max_it = 100
 ksp_gmres_restart = 50
 
 snes_rtol = 1.0e-20
-snes_atol = 1.0e-8
+snes_atol = 1.0e-9
 snes_max_it = 100
 
-ksp_monitor = true
+#ksp_monitor = true
 #ksp_view = true
-ksp_converged_reason = true
+#ksp_converged_reason = true
 
-snes_monitor = true
+#snes_monitor = true
 #snes_view = true
-snes_converged_reason = true
+#snes_converged_reason = true
 
 # Friction sensitivity solve used to compute the increment in slip
 # associated with changes in the Lagrange multiplier imposed by the
@@ -170,6 +172,7 @@ writer.filename = sheardispfriction.h5
 [sheardispfriction.timedependent.interfaces.fault.output]
 writer = pylith.meshio.DataWriterHDF5
 writer.filename = sheardispfriction-fault.h5
+vertex_info_fields = [strike_dir,normal_dir,traction_initial_value,friction_coefficient,cohesion]
 
 [sheardispfriction.timedependent.materials.elastic.output]
 cell_filter = pylith.meshio.CellFilterAvg