[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