[cig-commits] [commit] master: Use optimized integrators in bar_shearwave tri3 and tet4 examples. (422d933)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : master
Link       : https://github.com/geodynamics/pylith/compare/1b3d6d3bc246edc4235d0051142d675d91e9be41...422d933ebf5e19a69537cf35cab1c4497c03fe90

>---------------------------------------------------------------

commit 422d933ebf5e19a69537cf35cab1c4497c03fe90
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Wed Apr 9 16:12:52 2014 -0700

    Use optimized integrators in bar_shearwave tri3 and tet4 examples.


>---------------------------------------------------------------

422d933ebf5e19a69537cf35cab1c4497c03fe90
 examples/bar_shearwave/tet4/pylithapp.cfg | 20 ++++++++++----------
 examples/bar_shearwave/tri3/pylithapp.cfg | 21 ++++++++++-----------
 2 files changed, 20 insertions(+), 21 deletions(-)

diff --git a/examples/bar_shearwave/tet4/pylithapp.cfg b/examples/bar_shearwave/tet4/pylithapp.cfg
index f62e858..0f9ff50 100644
--- a/examples/bar_shearwave/tet4/pylithapp.cfg
+++ b/examples/bar_shearwave/tet4/pylithapp.cfg
@@ -57,8 +57,8 @@ filename_pset = bar_tet4_200m.pset
 dimension = 3
 elastic_prestep = False
 
-# Change to an explicit time stepping formulation
-formulation = pylith.problems.Explicit
+# Change to an explicit time stepping formulation optimized for tet4 cells
+formulation = pylith.problems.ExplicitTet4
 
 # Uncomment next line to use optimized lumped formulation and change
 # the quadrature order of the bulk material to 1.
@@ -100,10 +100,10 @@ db_properties.label = Elastic properties for +x material
 db_properties.iohandler.filename = matprops.spatialdb
 
 # Set the basis functions and quadrature:
-# 3-D simplex cell with 2nd order quadrature
+# 3-D simplex cell with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 3
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.materials.neg]
 
@@ -118,10 +118,10 @@ db_properties.label = Elastic properties for -x material
 db_properties.iohandler.filename = matprops.spatialdb
 
 # Set the basis functions and quadrature:
-# 3-D simplex cell with 2nd order quadrature
+# 3-D simplex cell with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 3
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # boundary conditions
@@ -137,10 +137,10 @@ db.label = Absorbing BC +x
 db.iohandler.filename = matprops.spatialdb
 
 # Specify the basis functions and quadrature:
-# 2-D simplex cell in 3-D space with 2nd order quadrature
+# 2-D simplex cell in 3-D space with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.x_neg]
 # Absorbing boundary condition on +x face of bar
@@ -153,10 +153,10 @@ db.label = Absorbing BC -x
 db.iohandler.filename = matprops.spatialdb
 
 # Specify the basis functions and quadrature:
-# 2-D simplex cell in 3-D space with 2nd order quadrature
+# 2-D simplex cell in 3-D space with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.all_nofault]
 # Dirichlet boundary condition on all vertices except fault vertices
diff --git a/examples/bar_shearwave/tri3/pylithapp.cfg b/examples/bar_shearwave/tri3/pylithapp.cfg
index 1287128..0e0df99 100644
--- a/examples/bar_shearwave/tri3/pylithapp.cfg
+++ b/examples/bar_shearwave/tri3/pylithapp.cfg
@@ -57,9 +57,8 @@ coordsys.space_dim = 2
 dimension = 2
 elastic_prestep = False
 
-# Change to an explicit time stepping formulation
-formulation = pylith.problems.Explicit
-#formulation = pylith.problems.ExplicitTri3
+# Change to an explicit time stepping formulation optimized for tri3 cells
+formulation = pylith.problems.ExplicitTri3
 
 # Nondimensionalize problem using wave propagation parameters.
 normalizer = spatialdata.units.NondimElasticDynamic
@@ -100,10 +99,10 @@ db_properties.label = Elastic properties
 db_properties.iohandler.filename = matprops.spatialdb
 
 # Set the basis functions and quadrature:
-# 2-D simplex cell with 2nd order quadrature
+# 2-D simplex cell with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 # ----------------------------------------------------------------------
 # boundary conditions
@@ -119,10 +118,10 @@ db.label = Absorbing BC +x
 db.iohandler.filename = matprops.spatialdb
 
 # Specify the basis functions and quadrature:
-# 1-D simplex cell in 2-D space with 2nd order quadrature
+# 1-D simplex cell in 2-D space with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.x_neg]
 # Absorbing boundary condition on -x face of bar
@@ -135,10 +134,10 @@ db.label = Absorbing BC -x
 db.iohandler.filename = matprops.spatialdb
 
 # Specify the basis functions and quadrature:
-# 1-D simplex cell in 2-D space with 2nd order quadrature
+# 1-D simplex cell in 2-D space with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 [pylithapp.timedependent.bc.all_nofault]
 # Dirichlet boundary condition on all vertices except fault vertices
@@ -172,10 +171,10 @@ id = 100
 label = fault
 
 # Specify the basis functions and quadrature:
-# 1-D simplex cell in 2-D space with 2nd order quadrature
+# 1-D simplex cell in 2-D space with 1st order quadrature
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 # Switch to Brune slip time function
 eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn