[cig-commits] r16926 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters

brad at geodynamics.org brad at geodynamics.org
Mon Jun 7 21:57:30 PDT 2010 

Author: brad
Date: 2010-06-07 21:57:30 -0700 (Mon, 07 Jun 2010)
New Revision: 16926

Added:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fieldsplit.cfg
Modified:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_unif_20km.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_graded_20km.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_unif_20km.cfg
Log:
Updated parameters.

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fieldsplit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fieldsplit.cfg	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/fieldsplit.cfg	2010-06-08 04:57:30 UTC (rev 16926)
@@ -0,0 +1,21 @@
+# -*- Python -*-
+[pylithapp]
+
+# Field split
+[pylithapp.timedependent.formulation]
+split_fields = True
+use_custom_constraint_pc = True
+matrix_type = aij
+
+[pylithapp.petsc]
+fs_pc_type = fieldsplit
+fs_pc_fieldsplit_real_diagonal = 
+fs_pc_fieldsplit_type = multiplicative
+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 = preonly
+fs_fieldsplit_1_ksp_type = preonly
+fs_fieldsplit_2_ksp_type = preonly
+fs_fieldsplit_3_ksp_type = preonly

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg	2010-06-08 04:31:28 UTC (rev 16925)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg	2010-06-08 04:57:30 UTC (rev 16926)
@@ -38,8 +38,8 @@
 [pylithapp.mesh_generator]
 # Change the mesh reader to the CUBIT reader.
 reader = pylith.meshio.MeshIOCubit
-use_nodeset_names = False
-#debug = 1   ; uncomment to get very verbose mesh information
+reader.use_nodeset_names = False
+reorder_mesh = true
 
 [pylithapp.mesh_generator.reader]
 coordsys.space_dim = 3
@@ -48,15 +48,11 @@
 # problem
 # ----------------------------------------------------------------------
 [pylithapp.timedependent]
-# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'z_neg',
-# 'y_pos', and 'y_neg'.
 bc = [x_pos,x_neg,z_neg,y_pos,y_neg]
 
-# Set interfaces to an array of 1 fault: 'fault'.
 interfaces = [fault]
 
 [pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
 # We will output the solution over the domain and the ground surface.
 output = [domain,subdomain]
 
@@ -144,27 +140,55 @@
 label = 10
 
 # ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+vertex_data_fields = [displacement,velocity]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# ----------------------------------------------------------------------
 # PETSc
 # ----------------------------------------------------------------------
 # Set the solver options.
 
 [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 =
+sub_pc_factor_shift_type = nonzero
 
+ksp_rtol = 1.0e-8
+ksp_max_it = 1000
+ksp_gmres_restart = 100
+ksp_monitor_singular_value = true
+ksp_converged_reason = true
+
 ksp_monitor = true
 ksp_view = true
+
 log_summary = true
-ksp_max_it = 1000
-ksp_gmres_restart = 100
-ksp_monitor_singular_value = true
 # start_in_debugger = true
 
-snes_monitor = true
-snes_view = true
-ksp_converged_reason = true
-snes_converged_reason = true
+# End of file

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg	2010-06-08 04:31:28 UTC (rev 16925)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_graded_20km.cfg	2010-06-08 04:57:30 UTC (rev 16926)
@@ -36,35 +36,19 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for VTK files.
 [pylithapp.problem.formulation.output.domain.writer]
 filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK domain output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
 writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
-# Give basename for vtk fault rupture output.
 [pylithapp.timedependent.interfaces.fault.output.writer]
 filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
 [pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_unif_20km.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_unif_20km.cfg	2010-06-08 04:31:28 UTC (rev 16925)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_hex8_unif_20km.cfg	2010-06-08 04:57:30 UTC (rev 16926)
@@ -36,35 +36,19 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for VTK files.
 [pylithapp.problem.formulation.output.domain.writer]
 filename = ../results/spbm_hex8_unif_20km/spbm_hex8_unif_20km.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK domain output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
 writer.filename = ../results/spbm_hex8_unif_20km/spbm_hex8_unif_20km-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
-# Give basename for vtk fault rupture output.
 [pylithapp.timedependent.interfaces.fault.output.writer]
 filename = ../results/spbm_hex8_unif_20km/spbm_hex8_unif_20km-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_hex8_unif_20km/spbm_hex8_unif_20km-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
 [pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_hex8_unif_20km/spbm_hex8_unif_20km-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_graded_20km.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_graded_20km.cfg	2010-06-08 04:31:28 UTC (rev 16925)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_graded_20km.cfg	2010-06-08 04:57:30 UTC (rev 16926)
@@ -36,35 +36,20 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for VTK files.
 [pylithapp.problem.formulation.output.domain.writer]
 filename = ../results/spbm_tet4_graded_20km/spbm_tet4_graded_20km.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK domain output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
 label = 17 ; nodeset for subdomain
 writer.filename = ../results/spbm_tet4_graded_20km/spbm_tet4_graded_20km-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
-# Give basename for vtk fault rupture output.
 [pylithapp.timedependent.interfaces.fault.output.writer]
 filename = ../results/spbm_tet4_graded_20km/spbm_tet4_graded_20km-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_tet4_graded_20km/spbm_tet4_graded_20km-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
 [pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_tet4_graded_20km/spbm_tet4_graded_20km-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_unif_20km.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_unif_20km.cfg	2010-06-08 04:31:28 UTC (rev 16925)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/spbm_tet4_unif_20km.cfg	2010-06-08 04:57:30 UTC (rev 16926)
@@ -36,35 +36,19 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for VTK files.
 [pylithapp.problem.formulation.output.domain.writer]
 filename = ../results/spbm_tet4_unif_20km/spbm_tet4_unif_20km.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK domain output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
 writer.filename = ../results/spbm_tet4_unif_20km/spbm_tet4_unif_20km-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
-# Give basename for vtk fault rupture output.
 [pylithapp.timedependent.interfaces.fault.output.writer]
 filename = ../results/spbm_tet4_unif_20km/spbm_tet4_unif_20km-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
 
-# Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_tet4_unif_20km/spbm_tet4_unif_20km-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
 [pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = ../results/spbm_tet4_unif_20km/spbm_tet4_unif_20km-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year

More information about the CIG-COMMITS mailing list