[cig-commits] r16866 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210

[brad at geodynamics.org]

Author: brad
Date: 2010-06-02 10:08:19 -0700 (Wed, 02 Jun 2010)
New Revision: 16866

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg
Log:
Use optimized tet4 lumped formulation.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg	2010-06-02 17:07:52 UTC (rev 16865)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg	2010-06-02 17:08:19 UTC (rev 16866)
@@ -4,43 +4,44 @@
 # ----------------------------------------------------------------------
 # materials
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.elastic]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
+[pylithapp.timedependent.materials.elastic.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 3
+cell.quad_order = 2
 
 
 # ----------------------------------------------------------------------
 # boundary conditions
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+[pylithapp.timedependent.bc.x_pos.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 2
+cell.quad_order = 2
 
-[pylithapp.timedependent.bc.x_neg]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+[pylithapp.timedependent.bc.x_neg.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 2
+cell.quad_order = 2
 
-[pylithapp.timedependent.bc.y_pos]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+[pylithapp.timedependent.bc.y_pos.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 2
+cell.quad_order = 2
 
-[pylithapp.timedependent.bc.y_neg]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+[pylithapp.timedependent.bc.y_neg.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 2
+cell.quad_order = 2
 
-[pylithapp.timedependent.bc.z_neg]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+[pylithapp.timedependent.bc.z_neg.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 2
+cell.quad_order = 2
 
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+[pylithapp.timedependent.interfaces.fault.quadrature]
+cell = pylith.feassemble.FIATLagrange
+cell.dimension = 2
+cell.quad_order = 2

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg	2010-06-02 17:07:52 UTC (rev 16865)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg	2010-06-02 17:08:19 UTC (rev 16866)
@@ -1,5 +1,4 @@
 # -*- Python -*-
-
 [pylithapp]
 
 # ----------------------------------------------------------------------
@@ -12,8 +11,6 @@
 meshiocubit = 1
 elasticityexplicit = 1
 meshimporter = 1
-#quadrature2d = 1
-#fiatlagrange = 1
 faultcohesivedyn = 1
 
 [pylithapp.journal.debug]
@@ -44,8 +41,6 @@
 # Set materials to an array with 1 material 'elastic'.
 materials = [elastic]
 
-#formulation = pylith.problems.Explicit
-#formulation.solver = pylith.problems.SolverNonlinear
 formulation = pylith.problems.ExplicitLumped
 
 # Nondimensionalize problem using wave propagation parameters.
@@ -169,40 +164,11 @@
 output_freq = time_step
 time_step = 30.0*s
 writer.time_format = %06.3f
+cell_info_fields = []
+cell_vertex_fields = []
 
 # ----------------------------------------------------------------------
 # PETSc
 # ----------------------------------------------------------------------
 [pylithapp.petsc]
-pc_type = asm
-
-# Change the preconditioner settings.
-sub_pc_factor_shift_type = nonzero
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-15
-ksp_max_it = 100
-ksp_gmres_restart = 50
-snes_max_it = 200
-
-ksp_monitor = true
-ksp_view = true
-ksp_converged_reason = true
-
-snes_monitor = true
-snes_view = true
-snes_converged_reason = true
-
 log_summary = true
-
-# Friction sensitivity solve
-friction_pc_type = asm
-friction_sub_pc_factor_shift_type = nonzero
-
-friction_ksp_max_it = 25
-friction_ksp_gmres_restart = 30
-
-friction_ksp_monitor = true
-friction_ksp_view = true
-friction_ksp_converged_reason = true
-

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg	2010-06-02 17:07:52 UTC (rev 16865)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg	2010-06-02 17:08:19 UTC (rev 16866)
@@ -2,39 +2,59 @@
 [pylithapp]
 
 # ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation = pylith.problems.ExplicitLumpedTet4
+
+# ----------------------------------------------------------------------
 # materials
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.elastic]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = tetrahedron
+[pylithapp.timedependent.materials.elastic.quadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = tetrahedron
+cell.degree = 1
+cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # boundary conditions
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+[pylithapp.timedependent.bc.x_pos.quadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = triangle
+cell.degree = 1
+cell.quad_order = 1
 
-[pylithapp.timedependent.bc.x_neg]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+[pylithapp.timedependent.bc.x_neg.quadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = triangle
+cell.degree = 1
+cell.quad_order = 1
 
-[pylithapp.timedependent.bc.y_pos]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+[pylithapp.timedependent.bc.y_pos.quadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = triangle
+cell.degree = 1
+cell.quad_order = 1
 
-[pylithapp.timedependent.bc.y_neg]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+[pylithapp.timedependent.bc.y_neg.qadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = triangle
+cell.degree = 1
+cell.quad_order = 1
 
-[pylithapp.timedependent.bc.z_neg]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+[pylithapp.timedependent.bc.z_neg.quadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = triangle
+cell.degree = 1
+cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+[pylithapp.timedependent.interfaces.fault.quadrature]
+cell = pylith.feassemble.FIATSimplex
+cell.shape = triangle
+cell.degree = 1
+cell.quad_order = 1