[cig-commits] r16248 - short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters

[brad at geodynamics.org]

Author: brad
Date: 2010-02-09 16:10:20 -0800 (Tue, 09 Feb 2010)
New Revision: 16248

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg
Log:
Updated parameters.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg	2010-02-10 00:09:40 UTC (rev 16247)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg	2010-02-10 00:10:20 UTC (rev 16248)
@@ -12,6 +12,7 @@
 meshiolagrit = 1
 explicitelasticity = 1
 faultcohesivekin = 1
+fiatlagrange = 1
 
 [pylithapp.journal.debug]
 pylithapp = 1
@@ -24,7 +25,7 @@
 [pylithapp.timedependent]
 dimension = 3
 
-formulation = pylith.problems.Explicit
+formulation = pylith.problems.ExplicitLumped
 
 bc = [x_pos,x_neg,y_pos,y_neg,z_neg]
 bc.x_pos = pylith.bc.AbsorbingDampers
@@ -39,7 +40,7 @@
 formulation.output = [domain,groundsurface]
 
 [pylithapp.timedependent.formulation.time_step]
-total_time = 15.0*s
+total_time = 0.25*s
 dt = 0.05*s
 
 # ----------------------------------------------------------------------
@@ -91,20 +92,6 @@
 # PETSc
 # ----------------------------------------------------------------------
 [pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings (must turn off
-# shift_positive_definite and turn on shift_nonzero).
-sub_pc_factor_shift_positive_definite = 0
-sub_pc_factor_shift_nonzero = 
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-10
-ksp_max_it = 100
-ksp_gmres_restart = 25
-
-ksp_monitor = true
-ksp_view = true
 log_summary = true
 
 

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg	2010-02-10 00:09:40 UTC (rev 16247)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg	2010-02-10 00:10:20 UTC (rev 16248)
@@ -21,13 +21,19 @@
 filename_pset = ../meshes/halfspace_tet4_lagrit.pset
 
 # ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation = pylith.problems.ExplicitLumpedTet4
+
+# ----------------------------------------------------------------------
 # materials
 # ----------------------------------------------------------------------
 [pylithapp.timedependent.materials.elastic]
 
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = tetrahedron
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # boundary conditions
@@ -36,31 +42,31 @@
 label = elemfaces_xmax
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.x_neg]
 label = elemfaces_xmin
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.y_pos]
 label = elemfaces_ymax
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.y_neg]
 label = elemfaces_ymin
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.z_neg]
 label = elemfaces_zmin
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # faults
@@ -69,7 +75,7 @@
 label = elemfaces_fault
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # output