[cig-commits] r15177 - in short/3D/PyLith/benchmarks/trunk: dynamic/halfspace_strikeslip/parameters quasistatic/reversenog/parameters quasistatic/strikeslipnog/parameters static/sceccvmh
brad at geodynamics.org
brad at geodynamics.org
Wed Jun 10 12:52:36 PDT 2009
Author: brad
Date: 2009-06-10 12:52:34 -0700 (Wed, 10 Jun 2009)
New Revision: 15177
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/fixeddisp_zero.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_maxwell.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0250m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0500m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0250m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0500m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/fixeddisp_zero.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_maxwell.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_cubit.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_lagrit.cfg
short/3D/PyLith/benchmarks/trunk/static/sceccvmh/mat_homogeneous.spatialdb
short/3D/PyLith/benchmarks/trunk/static/sceccvmh/slice_dislocation.cfg
short/3D/PyLith/benchmarks/trunk/static/sceccvmh/volume_dislocation.cfg
Log:
Updated settings.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -12,9 +12,9 @@
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
#debug = 1
-importer = pylith.meshio.MeshIOCubit
+reader = pylith.meshio.MeshIOCubit
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename = ../meshes/halfspace_hex8_cubit.exo
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/mat_elastic.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/mat_elastic.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -2,7 +2,7 @@
SimpleDB {
num-values = 3
value-names = Density Vp Vs
- value-units = kg/m3 m/s m/s
+ value-units = kg/m**3 m/s m/s
num-locs = 1
data-dim = 0
cs-data = cartesian {
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -25,7 +25,6 @@
dimension = 3
formulation = pylith.problems.Explicit
-formulation.solver.initial_guess_zero = False
bc = [x_pos,x_neg,y_pos,y_neg,z_neg]
bc.x_pos = pylith.bc.AbsorbingDampers
@@ -50,8 +49,7 @@
label = Elastic material
id = 1
-db.iohandler.filename = mat_elastic.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
+db_properties.iohandler.filename = mat_elastic.spatialdb
# ----------------------------------------------------------------------
# boundary conditions
@@ -86,9 +84,7 @@
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault]
id = 100
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-mat_db.iohandler.filename = mat_elastic.spatialdb
eq_srcs.rupture.slip_function = pylith.faults.LiuCosSlipFn
[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
@@ -102,6 +98,11 @@
[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
@@ -130,4 +131,4 @@
[pylithapp.timedependent.materials.elastic.output]
skip = 999999
writer.time_format = %05.2f
-cell_filter = pylith.meshio.CellFilterAvg
+cell_filter = pylith.meshio.CellFilterAvgMesh
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 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -14,9 +14,9 @@
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
#debug = 1
-importer = pylith.meshio.MeshIOLagrit
+reader = pylith.meshio.MeshIOLagrit
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename_gmv = ../meshes/halfspace_tet4_lagrit.gmv
filename_pset = ../meshes/halfspace_tet4_lagrit.pset
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/fixeddisp_zero.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/fixeddisp_zero.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/fixeddisp_zero.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -1,7 +1,7 @@
#SPATIAL.ascii 1
SimpleDB {
num-values = 3
- value-names = dof-0 dof-1 dof-2
+ value-names = displacement-x displacement-y displacement-z
value-units = m m m
num-locs = 1
data-dim = 0
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_elastic.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_elastic.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -2,7 +2,7 @@
SimpleDB {
num-values = 3
value-names = density vs vp
- value-units = kg/m^3 m/s m/s
+ value-units = kg/m**3 m/s m/s
num-locs = 1
data-dim = 0
space-dim = 3
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_maxwell.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_maxwell.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/mat_maxwell.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -2,7 +2,7 @@
SimpleDB {
num-values = 4
value-names = density vs vp viscosity
- value-units = kg/m^3 m/s m/s Pa*s
+ value-units = kg/m**3 m/s m/s Pa*s
num-locs = 1
data-dim = 0
space-dim = 3
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -26,7 +26,7 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
-debug = 0
+#debug = 1
[pylithapp.mesh_generator.importer]
coordsys.space_dim = 3
@@ -37,14 +37,9 @@
[pylithapp.timedependent]
dimension = 3
formulation = pylith.problems.Implicit
-
-# Array of DirichletPoints (default) boundary conditions
bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
-
-# Fault
interfaces = [fault]
-# Output solution over domain + ground surface (subdomain)
[pylithapp.timedependent.implicit]
output = [domain,subdomain]
output.subdomain = pylith.meshio.OutputSolnSubset
@@ -61,50 +56,48 @@
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.y_neg]
fixed_dof = [1]
-db = spatialdata.spatialdb.SimpleDB
-db.label = Symmetry Dirichlet BC
-db.iohandler.filename = fixeddisp_zero.spatialdb
-db.query_type = nearest
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Symmetry Dirichlet BC
+db_initial.iohandler.filename = fixeddisp_zero.spatialdb
+db_initial.query_type = nearest
[pylithapp.timedependent.bc.x_pos]
fixed_dof = [0, 1, 2]
-db = spatialdata.spatialdb.SimpleDB
-db.label = Analytic displacement Dirichlet BC on +x face
-#db.iohandler.filename = fixeddisp_analytic_xpos.spatialdb
-db.iohandler.filename = fixeddisp_zero.spatialdb
-db.query_type = nearest
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Analytic displacement Dirichlet BC on +x face
+#db_initial.iohandler.filename = fixeddisp_analytic_xpos.spatialdb
+db_initial.iohandler.filename = fixeddisp_zero.spatialdb
+db_initial.query_type = nearest
[pylithapp.timedependent.bc.x_neg]
fixed_dof = [0, 1, 2]
-db = spatialdata.spatialdb.SimpleDB
-db.label = Analytic displacement Dirichlet BC on -x face
-#db.iohandler.filename = fixeddisp_analytic_xneg.spatialdb
-db.iohandler.filename = fixeddisp_zero.spatialdb
-db.query_type = nearest
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Analytic displacement Dirichlet BC on -x face
+#db_initial.iohandler.filename = fixeddisp_analytic_xneg.spatialdb
+db_initial.iohandler.filename = fixeddisp_zero.spatialdb
+db_initial.query_type = nearest
[pylithapp.timedependent.bc.y_pos]
fixed_dof = [0, 1, 2]
-db = spatialdata.spatialdb.SimpleDB
-db.label = Analytic displacement Dirichlet BC on +y face
-#db.iohandler.filename = fixeddisp_analytic_ypos.spatialdb
-db.iohandler.filename = fixeddisp_zero.spatialdb
-db.query_type = nearest
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Analytic displacement Dirichlet BC on +y face
+#db_initial.iohandler.filename = fixeddisp_analytic_ypos.spatialdb
+db_initial.iohandler.filename = fixeddisp_zero.spatialdb
+db_initial.query_type = nearest
[pylithapp.timedependent.bc.z_neg]
fixed_dof = [0, 1, 2]
-db = spatialdata.spatialdb.SimpleDB
-db.label = Analytic displacement Dirichlet BC on -z face
-#db.iohandler.filename = fixeddisp_analytic_zneg.spatialdb
-db.iohandler.filename = fixeddisp_zero.spatialdb
-db.query_type = nearest
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Analytic displacement Dirichlet BC on -z face
+#db_initial.iohandler.filename = fixeddisp_analytic_zneg.spatialdb
+db_initial.iohandler.filename = fixeddisp_zero.spatialdb
+db_initial.query_type = nearest
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault]
id = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-mat_db.iohandler.filename = mat_elastic.spatialdb
[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
slip.label = Final slip
@@ -120,6 +113,11 @@
[pylithapp.petsc]
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_gmres_restart = 70
@@ -127,7 +125,7 @@
#pc_type = lu
#mat_type = aijmumps
-ksp_max_it = 500
+ksp_max_it = 100
ksp_monitor = true
ksp_view = true
log_summary = true
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -5,7 +5,7 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
-importer = pylith.meshio.MeshIOCubit
+reader = pylith.meshio.MeshIOCubit
# ----------------------------------------------------------------------
# problem
@@ -22,16 +22,14 @@
[pylithapp.timedependent.materials.elastic]
label = Elastic material
id = 1
-db.iohandler.filename = mat_elastic.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
+db_properties.iohandler.filename = mat_elastic.spatialdb
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
[pylithapp.timedependent.materials.viscoelastic]
label = Viscoelastic material (hanging wall)
id = 2
-db.iohandler.filename = mat_maxwell.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
+db_properties.iohandler.filename = mat_maxwell.spatialdb
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
@@ -60,6 +58,7 @@
label = 10
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0250m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0250m.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0250m.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -4,7 +4,7 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename = ../meshes/reverseslip_hex8_0250m.exo
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0500m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0500m.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_0500m.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -4,7 +4,7 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename = ../meshes/reverseslip_hex8_0500m.exo
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -4,7 +4,7 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename = ../meshes/reverseslip_hex8_1000m.exo
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -5,7 +5,7 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
-importer = pylith.meshio.MeshIOLagrit
+reader = pylith.meshio.MeshIOLagrit
# ----------------------------------------------------------------------
# problem
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0250m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0250m.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0250m.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -4,7 +4,7 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename_gmv = ../meshes/reverseslip_tet4_0250m.gmv
filename_pset = ../meshes/reverseslip_tet4_0250m.pset
io_int32 = False
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0500m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0500m.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_0500m.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -4,7 +4,7 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename_gmv = ../meshes/reverseslip_tet4_0500m.gmv
filename_pset = ../meshes/reverseslip_tet4_0500m.pset
io_int32 = False
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -4,7 +4,7 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename_gmv = ../meshes/reverseslip_tet4_1000m.gmv
filename_pset = ../meshes/reverseslip_tet4_1000m.pset
io_int32 = False
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/fixeddisp_zero.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/fixeddisp_zero.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/fixeddisp_zero.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -1,7 +1,7 @@
#SPATIAL.ascii 1
SimpleDB {
num-values = 3
- value-names = dof-0 dof-1 dof-2
+ value-names = displacement-x displacement-y displacement-z
value-units = m m m
num-locs = 1
data-dim = 0
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_elastic.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_elastic.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -2,7 +2,7 @@
SimpleDB {
num-values = 3
value-names = density vs vp
- value-units = kg/m^3 m/s m/s
+ value-units = kg/m**3 m/s m/s
num-locs = 1
data-dim = 0
space-dim = 3
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_maxwell.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_maxwell.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/mat_maxwell.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -2,7 +2,7 @@
SimpleDB {
num-values = 4
value-names = density vs vp viscosity
- value-units = kg/m^3 m/s m/s Pa*s
+ value-units = kg/m**3 m/s m/s Pa*s
num-locs = 1
data-dim = 0
space-dim = 3
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/pylithapp.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/pylithapp.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -26,7 +26,7 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
-debug = 0
+#debug = 1
[pylithapp.mesh_generator.reader]
coordsys.space_dim = 3
@@ -105,6 +105,11 @@
[pylithapp.petsc]
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_gmres_restart = 70
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_cubit.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_cubit.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -22,14 +22,14 @@
[pylithapp.timedependent.materials.elastic]
label = Elastic material
id = 1
-properties_db.iohandler.filename = mat_elastic.spatialdb
+db_properties.iohandler.filename = mat_elastic.spatialdb
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
[pylithapp.timedependent.materials.viscoelastic]
label = Viscoelastic material
id = 2
-properties_db.iohandler.filename = mat_maxwell.spatialdb
+db_properties.iohandler.filename = mat_maxwell.spatialdb
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_lagrit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_lagrit.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/parameters/strikeslip_lagrit.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -23,28 +23,28 @@
[pylithapp.timedependent.materials.elastic_pos]
label = Elastic material +x
id = 1
-properties_db.iohandler.filename = mat_elastic.spatialdb
+db_properties.iohandler.filename = mat_elastic.spatialdb
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
[pylithapp.timedependent.materials.elastic_neg]
label = Elastic material -x
id = 2
-properties_db.iohandler.filename = mat_elastic.spatialdb
+db_properties.iohandler.filename = mat_elastic.spatialdb
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
[pylithapp.timedependent.materials.viscoelastic_pos]
label = Viscoelastic material +x
id = 3
-properties_db.iohandler.filename = mat_maxwell.spatialdb
+db_properties.iohandler.filename = mat_maxwell.spatialdb
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
[pylithapp.timedependent.materials.viscoelastic_neg]
label = Viscoelastic material -x
id = 4
-properties_db.iohandler.filename = mat_maxwell.spatialdb
+db_properties.iohandler.filename = mat_maxwell.spatialdb
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
Modified: short/3D/PyLith/benchmarks/trunk/static/sceccvmh/mat_homogeneous.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/sceccvmh/mat_homogeneous.spatialdb 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/static/sceccvmh/mat_homogeneous.spatialdb 2009-06-10 19:52:34 UTC (rev 15177)
@@ -2,7 +2,7 @@
SimpleDB {
num-values = 3
value-names = density vs vp
- value-units = kg/m^3 m/s m/s
+ value-units = kg/m**3 m/s m/s
num-locs = 1
data-dim = 0
space-dim = 3
Modified: short/3D/PyLith/benchmarks/trunk/static/sceccvmh/slice_dislocation.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/sceccvmh/slice_dislocation.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/static/sceccvmh/slice_dislocation.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -19,9 +19,9 @@
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
#debug = 1
-importer = pylith.meshio.MeshIOCubit
+reader = pylith.meshio.MeshIOCubit
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename = slice_hex8.exo
coordsys = spatialdata.geocoords.CSGeoProj
@@ -58,7 +58,6 @@
#db = spatialdata.spatialdb.SimpleDB
#db.iohandler.filename = mat_homogeneous.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
quadrature.cell.quad_order = 2
@@ -69,34 +68,34 @@
[pylithapp.timedependent.bc.bc_east]
fixed_dof = [0, 1]
label = 21
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC east surface
-db.values = [dof-0,dof-1]
-db.data = [0.0,-2.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC east surface
+db_initial.values = [displacement-x,displacement-y]
+db_initial.data = [0.0*m,-2.0*m]
[pylithapp.timedependent.bc.bc_west]
fixed_dof = [0, 1]
label = 22
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC west surface
-db.values = [dof-0,dof-1]
-db.data = [0.0,+2.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC west surface
+db_initial.values = [displacement-x,displacement-y]
+db_initial.data = [0.0*m,+2.0*m]
[pylithapp.timedependent.bc.bc_bottom]
fixed_dof = [2]
label = 26
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC bottom surface
-db.values = [dof-2]
-db.data = [0.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC bottom surface
+db_initial.values = [displacement-z]
+db_initial.data = [0.0*m]
[pylithapp.timedependent.bc.bc_ground]
fixed_dof = [2]
label = 25
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC ground surface
-db.values = [dof-2]
-db.data = [0.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC ground surface
+db_initial.values = [displacement-z]
+db_initial.data = [0.0*m]
# ----------------------------------------------------------------------
# faults
@@ -107,18 +106,10 @@
[pylithapp.timedependent.interfaces.fault]
id = 2
label = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
quadrature.cell.quad_order = 2
-#mat_db = spatialdata.spatialdb.SCECCVMH
-#mat_db.data_dir = /home/brad/data/sceccvm-h/vx53/bin
-#mat_db.squash = True
-#mat_db.squash_limit = -8.0*km
-mat_db = spatialdata.spatialdb.SimpleDB
-mat_db.iohandler.filename = mat_homogeneous.spatialdb
-
[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
slip.iohandler.filename = dislocation_slip.spatialdb
slip_time.iohandler.filename = dislocation_sliptime.spatialdb
@@ -130,6 +121,11 @@
[pylithapp.petsc]
ksp_rtol = 1.0e-8
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_max_it = 200
ksp_gmres_restart = 50
ksp_monitor = true
@@ -146,8 +142,8 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = slice-dislocation-statevars-elastic.vtk
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = slice-dislocation-elastic.vtk
[pylithapp.timedependent.interfaces.fault.output.writer]
filename = slice-dislocation-fault.vtk
Modified: short/3D/PyLith/benchmarks/trunk/static/sceccvmh/volume_dislocation.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/sceccvmh/volume_dislocation.cfg 2009-06-10 19:21:04 UTC (rev 15176)
+++ short/3D/PyLith/benchmarks/trunk/static/sceccvmh/volume_dislocation.cfg 2009-06-10 19:52:34 UTC (rev 15177)
@@ -19,9 +19,9 @@
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
#debug = 1
-importer = pylith.meshio.MeshIOCubit
+reader = pylith.meshio.MeshIOCubit
-[pylithapp.mesh_generator.importer]
+[pylithapp.mesh_generator.reader]
filename = volume_hex8.exo
coordsys = spatialdata.geocoords.CSGeoProj
@@ -58,7 +58,6 @@
#db = spatialdata.spatialdb.SimpleDB
#db.iohandler.filename = mat_homogeneous.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
quadrature.cell.quad_order = 2
@@ -69,26 +68,26 @@
[pylithapp.timedependent.bc.bc_east]
fixed_dof = [0, 1]
label = 21
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC east surface
-db.values = [dof-0,dof-1]
-db.data = [0.0,-2.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC east surface
+db_initial.values = [displacement-x,displacement-y]
+db_initial.data = [0.0*m,-2.0*m]
[pylithapp.timedependent.bc.bc_west]
fixed_dof = [0, 1]
label = 22
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC west surface
-db.values = [dof-0,dof-1]
-db.data = [0.0,+2.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC west surface
+db_initial.values = [displacement-x,displacement-y]
+db_initial.data = [0.0*m,+2.0*m]
[pylithapp.timedependent.bc.bc_bottom]
fixed_dof = [2]
label = 26
-db = spatialdata.spatialdb.UniformDB
-db.label = Dirichlet BC bottom surface
-db.values = [dof-2]
-db.data = [0.0]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC bottom surface
+db_initial.values = [displacement-z]
+db_initial.data = [0.0*m]
# ----------------------------------------------------------------------
# faults
@@ -99,18 +98,10 @@
[pylithapp.timedependent.interfaces.fault]
id = 2
label = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
quadrature.cell.quad_order = 2
-#mat_db = spatialdata.spatialdb.SCECCVMH
-#mat_db.data_dir = /home/brad/data/sceccvm-h/vx53/bin
-#mat_db.squash = True
-#mat_db.squash_limit = -8.0*km
-mat_db = spatialdata.spatialdb.SimpleDB
-mat_db.iohandler.filename = mat_homogeneous.spatialdb
-
[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
slip.iohandler.filename = dislocation_slip.spatialdb
slip_time.iohandler.filename = dislocation_sliptime.spatialdb
@@ -122,6 +113,11 @@
[pylithapp.petsc]
ksp_rtol = 1.0e-8
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_max_it = 200
ksp_gmres_restart = 50
ksp_monitor = true
@@ -138,8 +134,8 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = volume-dislocation-statevars-elastic.vtk
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = volume-dislocation-elastic.vtk
[pylithapp.timedependent.interfaces.fault.output.writer]
filename = volume-dislocation-fault.vtk
More information about the CIG-COMMITS
mailing list