[cig-commits] r16635 - in short/3D/PyLith/benchmarks/trunk: dynamic/scecdynrup/tpv205-2d quasistatic static static/solvertest/parameters
brad at geodynamics.org
brad at geodynamics.org
Tue May 4 21:12:13 PDT 2010
Author: brad
Date: 2010-05-04 21:12:13 -0700 (Tue, 04 May 2010)
New Revision: 16635
Added:
short/3D/PyLith/benchmarks/trunk/static/solvertest/
short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg
Removed:
short/3D/PyLith/benchmarks/trunk/quasistatic/solvertest/
short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.jou
Log:
Moved solvertest to static.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.jou 2010-05-05 04:10:45 UTC (rev 16634)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.jou 2010-05-05 04:12:13 UTC (rev 16635)
@@ -26,5 +26,3 @@
# Export exodus file
# ----------------------------------------------------------------------
export mesh "quad4_100m.exo" dimension 2 overwrite
-
-
Copied: short/3D/PyLith/benchmarks/trunk/static/solvertest (from rev 16599, short/3D/PyLith/benchmarks/trunk/quasistatic/solvertest)
Deleted: short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/solvertest/parameters/pylithapp.cfg 2010-04-29 04:04:17 UTC (rev 16599)
+++ short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg 2010-05-05 04:12:13 UTC (rev 16635)
@@ -1,216 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 1
-implicit = 1
-petsc = 1
-solverlinear = 1
-meshimporter = 1
-meshiocubit = 1
-elasticityimplicit = 1
-faultcohesivekin = 1
-
-[pylithapp.journal.debug]
-pylithapp = 1
-problem = 1
-implicit = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-reader = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.reader]
-coordsys.space_dim = 3
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 3
-
-materials = [elastic,viscoelastic]
-bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
-interfaces = [fault_middle,fault_xpos,fault_xneg]
-
-normalizer.length_scale = 1.0*m
-
-[pylithapp.timedependent.formulation]
-output = [domain,subdomain]
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 0.0*s
-dt = 1.0*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-viscoelastic = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-
-[pylithapp.timedependent.materials.viscoelastic]
-label = Viscoelastic material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.y_neg]
-label = 13
-bc_dof = [0]
-db_initial = pylith.bc.ZeroDispDB
-db_initial.label = Dirichlet BC on -y
-
-[pylithapp.timedependent.bc.y_pos]
-label = 12
-bc_dof = [0]
-db_initial = pylith.bc.ZeroDispDB
-db_initial.label = Dirichlet BC on +y
-
-[pylithapp.timedependent.bc.x_pos]
-label = 10
-bc_dof = [1]
-db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Dirichlet BC on +x face
-db_initial.values = [displacement-y]
-db_initial.data = [-2.0*m]
-
-[pylithapp.timedependent.bc.x_neg]
-label = 11
-bc_dof = [1]
-db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Dirichlet BC on -x face
-db_initial.values = [displacement-y]
-db_initial.data = [+2.0*m]
-
-[pylithapp.timedependent.bc.z_neg]
-label = 14
-bc_dof = [2]
-db_initial = pylith.bc.ZeroDispDB
-db_initial.label = Dirichlet BC on -z face
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault_middle]
-id = 100
-label = 100
-
-[pylithapp.timedependent.interfaces.fault_middle.eq_srcs.rupture.slip_function]
-slip.label = Final slip on middle fault
-slip = spatialdata.spatialdb.UniformDB
-slip.values = [left-lateral-slip,reverse-slip,fault-opening]
-slip.data = [-1.0*m,0.0*m,0.0*m]
-slip_time = spatialdata.spatialdb.UniformDB
-slip_time.values = [slip-time]
-slip_time.data = [0.0*year]
-
-[pylithapp.timedependent.interfaces.fault_xpos]
-id = 101
-label = 101
-
-[pylithapp.timedependent.interfaces.fault_xpos.eq_srcs.rupture.slip_function]
-slip.label = Final slip on +x fault
-slip = spatialdata.spatialdb.UniformDB
-slip.values = [left-lateral-slip,reverse-slip,fault-opening]
-slip.data = [-0.5*m,0.0*m,0.0*m]
-slip_time = spatialdata.spatialdb.UniformDB
-slip_time.values = [slip-time]
-slip_time.data = [0.0*year]
-
-
-[pylithapp.timedependent.interfaces.fault_xneg]
-id = 102
-label = 102
-
-[pylithapp.timedependent.interfaces.fault_xneg.eq_srcs.rupture.slip_function]
-slip.label = Final slip on -x fault
-slip = spatialdata.spatialdb.UniformDB
-slip.values = [left-lateral-slip,reverse-slip,fault-opening]
-slip.data = [-0.5*m,0.0*m,0.0*m]
-slip_time = spatialdata.spatialdb.UniformDB
-slip_time.values = [slip-time]
-slip_time.data = [0.0*year]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-[pylithapp.problem.formulation.output.domain]
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-[pylithapp.problem.formulation.output.subdomain]
-label = 15
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-[pylithapp.problem.interfaces.fault_middle.output]
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-[pylithapp.problem.interfaces.fault_xpos.output]
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-[pylithapp.problem.interfaces.fault_xneg.output]
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-[pylithapp.problem.materials.elastic.output]
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-[pylithapp.problem.materials.viscoelastic.output]
-output_freq = time_step
-time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[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
-ksp_gmres_preallocate = true
-
-#ksp_type = preonly
-#pc_type = lu
-#mat_type = aijmumps
-
-ksp_max_it = 500
-ksp_monitor = true
-ksp_view = true
-log_summary = true
Copied: short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg (from rev 16634, short/3D/PyLith/benchmarks/trunk/quasistatic/solvertest/parameters/pylithapp.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/solvertest/parameters/pylithapp.cfg 2010-05-05 04:12:13 UTC (rev 16635)
@@ -0,0 +1,230 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshimporter = 1
+meshiocubit = 1
+elasticityimplicit = 1
+faultcohesivekin = 1
+
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+implicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+
+materials = [elastic,viscoelastic]
+bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
+interfaces = [fault_middle,fault_xpos,fault_xneg]
+
+normalizer.length_scale = 1.0*m
+
+[pylithapp.timedependent.formulation]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+dt = 1.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+viscoelastic = pylith.materials.MaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = mat_elastic.spatialdb
+
+[pylithapp.timedependent.materials.viscoelastic]
+label = Viscoelastic material
+id = 2
+db_properties.iohandler.filename = mat_maxwell.spatialdb
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.y_neg]
+label = 13
+bc_dof = [0]
+db_initial = pylith.bc.ZeroDispDB
+db_initial.label = Dirichlet BC on -y
+
+[pylithapp.timedependent.bc.y_pos]
+label = 12
+bc_dof = [0]
+db_initial = pylith.bc.ZeroDispDB
+db_initial.label = Dirichlet BC on +y
+
+[pylithapp.timedependent.bc.x_pos]
+label = 10
+bc_dof = [1]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC on +x face
+db_initial.values = [displacement-y]
+db_initial.data = [-2.0*m]
+
+[pylithapp.timedependent.bc.x_neg]
+label = 11
+bc_dof = [1]
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC on -x face
+db_initial.values = [displacement-y]
+db_initial.data = [+2.0*m]
+
+[pylithapp.timedependent.bc.z_neg]
+label = 14
+bc_dof = [2]
+db_initial = pylith.bc.ZeroDispDB
+db_initial.label = Dirichlet BC on -z face
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault_middle]
+id = 100
+label = 100
+
+[pylithapp.timedependent.interfaces.fault_middle.eq_srcs.rupture.slip_function]
+slip.label = Final slip on middle fault
+slip = spatialdata.spatialdb.UniformDB
+slip.values = [left-lateral-slip,reverse-slip,fault-opening]
+slip.data = [-1.0*m,0.0*m,0.0*m]
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.values = [slip-time]
+slip_time.data = [0.0*year]
+
+[pylithapp.timedependent.interfaces.fault_xpos]
+id = 101
+label = 101
+
+[pylithapp.timedependent.interfaces.fault_xpos.eq_srcs.rupture.slip_function]
+slip.label = Final slip on +x fault
+slip = spatialdata.spatialdb.UniformDB
+slip.values = [left-lateral-slip,reverse-slip,fault-opening]
+slip.data = [-0.5*m,0.0*m,0.0*m]
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.values = [slip-time]
+slip_time.data = [0.0*year]
+
+
+[pylithapp.timedependent.interfaces.fault_xneg]
+id = 102
+label = 102
+
+[pylithapp.timedependent.interfaces.fault_xneg.eq_srcs.rupture.slip_function]
+slip.label = Final slip on -x fault
+slip = spatialdata.spatialdb.UniformDB
+slip.values = [left-lateral-slip,reverse-slip,fault-opening]
+slip.data = [-0.5*m,0.0*m,0.0*m]
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.values = [slip-time]
+slip_time.data = [0.0*year]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+[pylithapp.problem.formulation.output.subdomain]
+label = 15
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+[pylithapp.problem.interfaces.fault_middle.output]
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+[pylithapp.problem.interfaces.fault_xpos.output]
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+[pylithapp.problem.interfaces.fault_xneg.output]
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+[pylithapp.problem.materials.elastic.output]
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+[pylithapp.problem.materials.viscoelastic.output]
+output_freq = time_step
+time_step = 0.999*year
+writer.time_format = %04.0f
+writer.time_constant = 0.1*year
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+
+pc_type = asm
+
+sub_pc_factor_shift_type = nonzero
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-20
+ksp_gmres_restart = 120
+ksp_max_it = 100
+
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+
+log_summary = true
+
+
+# Field split
+#[pylithapp.timedependent.formulation]
+#split_fields = true
+
+#[pylithapp.petsc]
+#fs_pc_fieldsplit_real_diagonal =
+#fs_pc_fieldsplit_type = additive
+#fs_fieldsplit_0_pc_type = ml
+#fs_fieldsplit_1_pc_type = ml
+#fs_fieldsplit_2_pc_type = ml
+#fs_fieldsplit_3_pc_type = ml
+#fs_fieldsplit_0_ksp_type = gmres
+#fs_fieldsplit_1_ksp_type = gmres
+#fs_fieldsplit_2_ksp_type = gmres
+#fs_fieldsplit_3_ksp_type = jacobi
+
More information about the CIG-COMMITS
mailing list