[cig-commits] r16582 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205

brad at geodynamics.org brad at geodynamics.org
Mon Apr 26 15:42:21 PDT 2010


Author: brad
Date: 2010-04-26 15:42:20 -0700 (Mon, 26 Apr 2010)
New Revision: 16582

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/test.cfg
Log:
Updated settings.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/test.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/test.cfg	2010-04-26 16:02:03 UTC (rev 16581)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/test.cfg	2010-04-26 22:42:20 UTC (rev 16582)
@@ -28,7 +28,7 @@
 reader = pylith.meshio.MeshIOCubit
 
 [pylithapp.mesh_generator.reader]
-filename = SCEC2004.exo
+filename = test.exo
 coordsys.space_dim = 3
 
 # ----------------------------------------------------------------------
@@ -50,9 +50,6 @@
 
 # Nondimensionalize problem using wave propagation parameters.
 normalizer = spatialdata.units.NondimElasticDynamic
-#normalizer.mass_density = 1000.0*kg/m**3
-#normalizer.shear_wave_speed = 1.0*m/s
-#normalizer.wave_period = 1.0*s
  
 bc = [x_neg,x_pos,z_neg,y_neg,y_pos]
 
@@ -75,12 +72,11 @@
 
 materials = [elastic]
 materials.elastic = pylith.materials.ElasticIsotropic3D
-#materials.elastic = pylith.materials.Homogeneous
 
 [pylithapp.timedependent.materials.elastic]
 label = Elastic material
 id = 1
-db_properties.iohandler.filename = SCEC2004matprops.spatialdb
+db_properties.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
 
@@ -91,14 +87,14 @@
 [pylithapp.timedependent.bc.x_pos]
 label = 26
 db.label = Absorbing BC +x
-db.iohandler.filename = SCEC2004matprops.spatialdb
+db.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
 [pylithapp.timedependent.bc.x_neg]
 label = 27
 db.label = Absorbing BC -x
-db.iohandler.filename = SCEC2004matprops.spatialdb
+db.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
@@ -106,21 +102,21 @@
 label = 28
 up_dir = [0, 1, 0]
 db.label = Absorbing BC -z
-db.iohandler.filename = SCEC2004matprops.spatialdb
+db.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
 [pylithapp.timedependent.bc.y_pos]
 label = 22
 db.label = Absorbing BC +y
-db.iohandler.filename = SCEC2004matprops.spatialdb
+db.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
 [pylithapp.timedependent.bc.y_neg]
 label = 23
 db.label = Absorbing BC -y
-db.iohandler.filename = SCEC2004matprops.spatialdb
+db.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
@@ -149,7 +145,7 @@
 quadrature.cell.quad_order = 2
 
 db_initial_tractions = spatialdata.spatialdb.SimpleDB
-db_initial_tractions.iohandler.filename = SCEC2004initTract.spatialdb
+db_initial_tractions.iohandler.filename = tractions.spatialdb
 db_initial_tractions.label = "Initial fault tractions"
 
 friction = pylith.friction.SlipWeakening
@@ -164,7 +160,6 @@
 # PETSc
 # ----------------------------------------------------------------------
 [pylithapp.petsc]
-ksp_type = gmres
 pc_type = asm
 
 # Change the preconditioner settings.
@@ -172,35 +167,34 @@
 
 ksp_rtol = 1.0e-8
 ksp_atol = 1.0e-15
-ksp_max_it = 200
-ksp_gmres_restart = 250
-snes_max_it = 300
+ksp_max_it = 100
+ksp_gmres_restart = 50
+snes_max_it = 200
 
-#ksp_monitor = true
-#ksp_view = true
+ksp_monitor = true
+ksp_view = true
 ksp_converged_reason = true
 
-snes_rtol = 1.0e-8
-snes_atol = 1.0e-15
 snes_monitor = true
-#snes_view = true
+snes_view = true
 snes_converged_reason = true
 
 #log_summary = true
+#start_in_debugger = true
 
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
 # Give basename for VTK domain output of solution over domain.
 [pylithapp.problem.formulation.output.output]
-writer.filename = SCEC2004.vtk
+writer.filename = output/test.vtk
 
 # Give basename for VTK fault output.
 [pylithapp.timedependent.interfaces.fault.output]
 vertex_info_fields = [strike_dir,normal_dir,initial_traction]
-writer.filename = SCEC2004-fault.vtk
+writer.filename = output/test-fault.vtk
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.elastic.output]
 cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = SCEC2004-elastic.vtk
+writer.filename = output/test-elastic.vtk



More information about the CIG-COMMITS mailing list