[cig-commits] r19008 - short/3D/PyLith/branches/v1.6-revisedfault/tests/3d/slipdir

[brad at geodynamics.org]

Author: brad
Date: 2011-10-04 10:02:40 -0700 (Tue, 04 Oct 2011)
New Revision: 19008

Modified:
   short/3D/PyLith/branches/v1.6-revisedfault/tests/3d/slipdir/pylithapp.cfg
Log:
Added required spatial database labels.

Modified: short/3D/PyLith/branches/v1.6-revisedfault/tests/3d/slipdir/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/tests/3d/slipdir/pylithapp.cfg	2011-10-04 16:41:52 UTC (rev 19007)
+++ short/3D/PyLith/branches/v1.6-revisedfault/tests/3d/slipdir/pylithapp.cfg	2011-10-04 17:02:40 UTC (rev 19008)
@@ -20,19 +20,12 @@
 [pylithapp.mesh_generator]
 reader = pylith.meshio.MeshIOCubit
 
-[pylithapp.mesh_generator.reader]
-
 # ----------------------------------------------------------------------
 # problem
 # ----------------------------------------------------------------------
 [pylithapp.timedependent]
 dimension = 3
 
-bc = [boundary]
-interfaces = [fault]
-
-materials = [elastic]
-
 [pylithapp.timedependent.formulation.time_step]
 total_time = 0.0*s
 dt = 1.0*year
@@ -40,6 +33,9 @@
 # ----------------------------------------------------------------------
 # materials
 # ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+
 [pylithapp.timedependent.materials]
 elastic = pylith.materials.ElasticIsotropic3D
 
@@ -48,6 +44,7 @@
 label = Elastic material
 id = 1
 
+db_properties.label = Elastic properties
 db_properties.iohandler.filename = matprops.spatialdb
 
 quadrature.cell = pylith.feassemble.FIATSimplex
@@ -56,14 +53,20 @@
 # ----------------------------------------------------------------------
 # boundary conditions
 # ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [boundary]
+
 [pylithapp.timedependent.bc.boundary]
 
-bc_dof = [0,2]
+bc_dof = [0,1,2]
 label = boundary_nofault
 
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault]
+
 [pylithapp.timedependent.interfaces.fault]
 
 id = 100
@@ -71,12 +74,16 @@
 
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
 
 eq_srcs.rupture.slip_function = pylith.faults.StepSlipFn
 
 [pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
 
+slip.label = Final slip
 slip.iohandler.filename = slip.spatialdb
+
+slip_time.label = Slip time
 slip_time.iohandler.filename = sliptime.spatialdb
 
 # ----------------------------------------------------------------------
@@ -84,20 +91,39 @@
 # ----------------------------------------------------------------------
 [pylithapp.petsc]
 ksp_type = gmres
-pc_type = asm
 
-sub_pc_factor_shift_type = nonzero
-
 ksp_rtol = 1.0e-8
 ksp_atol = 1.0e-10
-ksp_max_it = 50
-ksp_gmres_restart = 25
+ksp_max_it = 100
+ksp_gmres_restart = 100
 
 ksp_monitor = true
 #ksp_view = true
 #log_summary = true
 
+[pylithapp.timedependent.formulation]
+split_fields = True
+matrix_type = aij
 
+[pylithapp.petsc]
+ksp_gmres_restart = 500
+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
+
+[pylithapp.timedependent.formulation]
+use_custom_constraint_pc = True
+
+
+
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------