[cig-commits] [commit] baagaard/fix-cubit-cellsize, baagaard/fix-fullscale-3d-tests, baagaard/fix-packaging, knepley/fix-parallel-mult-faults, knepley/upgrade-petsc-interface, master, next: Added embedded fault + refinement full-scale tet4 test. (096d4cb)

Wed Apr 9 03:02:32 PDT 2014

Repository : ssh://geoshell/pylith

On branches: baagaard/fix-cubit-cellsize,baagaard/fix-fullscale-3d-tests,baagaard/fix-packaging,knepley/fix-parallel-mult-faults,knepley/upgrade-petsc-interface,master,next
Link       : https://github.com/geodynamics/pylith/compare/a213c3005450d915f40c7137ff7d8dbbb439d334...1b3d6d3bc246edc4235d0051142d675d91e9be41

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

commit 096d4cbeaf5370d43b58aa6c5bd55ed51e51d3ec
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Tue Mar 11 13:23:25 2014 -0700

    Added embedded fault + refinement full-scale tet4 test.
    
    Exposes error.


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

096d4cbeaf5370d43b58aa6c5bd55ed51e51d3ec
 ...rDispNoSlip.py => TestShearDispNoSlipRefine.py} | 20 +++++---
 ...eardispnoslip.cfg => sheardispnosliprefine.cfg} | 57 +++++++++++-----------
 2 files changed, 42 insertions(+), 35 deletions(-)

diff --git a/tests_auto/3d/tet4/TestShearDispNoSlip.py b/tests_auto/3d/tet4/TestShearDispNoSlipRefine.py
similarity index 85%
copy from tests_auto/3d/tet4/TestShearDispNoSlip.py
copy to tests_auto/3d/tet4/TestShearDispNoSlipRefine.py
index 795312f..f72fc5b 100644
--- a/tests_auto/3d/tet4/TestShearDispNoSlip.py
+++ b/tests_auto/3d/tet4/TestShearDispNoSlipRefine.py
@@ -16,10 +16,10 @@
 # ----------------------------------------------------------------------
 #
 
-## @file tests/3d/tet4/TestShearDispNoSlip.py
+## @file tests/3d/tet4/TestShearDispNoSlipRefine.py
 ##
 ## @brief Test suite for testing pylith with 3-D shear motion with no
-## fault slip.
+## fault slip and mesh refinement.
 
 import numpy
 from TestTet4 import TestTet4
@@ -29,7 +29,7 @@ from sheardisp_soln import AnalyticalSoln
 from pylith.apps.PyLithApp import PyLithApp
 class ShearApp(PyLithApp):
   def __init__(self):
-    PyLithApp.__init__(self, name="sheardispnoslip")
+    PyLithApp.__init__(self, name="sheardispnosliprefine")
     return
 
 
@@ -51,7 +51,7 @@ def run_pylith():
   return
 
 
-class TestShearDispNoSlip(TestTet4):
+class TestShearDispNoSlipRefine(TestTet4):
   """
   Test suite for testing pylith with 2-D shear extension.
   """
@@ -61,14 +61,20 @@ class TestShearDispNoSlip(TestTet4):
     Setup for test.
     """
     TestTet4.setUp(self)
+    self.mesh = {'ncells-elastic': 7068*8,
+                 'ncells-viscoelastic': 11131*8,
+                 'ncorners': 4,
+                 'nvertices': 3652,
+                 'spaceDim': 3,
+                 'tensorSize': 6}
     self.nverticesO = self.mesh['nvertices']
     self.mesh['nvertices'] += 51
     self.faultMesh = {'nvertices': 51,
                       'spaceDim': 3,
-                      'ncells': 30,
+                      'ncells': 30*4,
                       'ncorners': 3}
     run_pylith()
-    self.outputRoot = "sheardispnoslip"
+    self.outputRoot = "sheardispnosliprefine"
 
     self.soln = AnalyticalSoln()
     return
@@ -138,7 +144,7 @@ class TestShearDispNoSlip(TestTet4):
 # ----------------------------------------------------------------------
 if __name__ == '__main__':
   import unittest
-  from TestShearDispNoSlip import TestShearDispNoSlip as Tester
+  from TestShearDispNoSlipRefine import TestShearDispNoSlipRefine as Tester
 
   suite = unittest.TestSuite()
   suite.addTest(unittest.makeSuite(Tester))
diff --git a/tests_auto/3d/tet4/sheardispnoslip.cfg b/tests_auto/3d/tet4/sheardispnosliprefine.cfg
similarity index 73%
copy from tests_auto/3d/tet4/sheardispnoslip.cfg
copy to tests_auto/3d/tet4/sheardispnosliprefine.cfg
index bc7a2c7..a3084cd 100644
--- a/tests_auto/3d/tet4/sheardispnoslip.cfg
+++ b/tests_auto/3d/tet4/sheardispnosliprefine.cfg
@@ -1,12 +1,12 @@
-[sheardispnoslip]
+[sheardispnosliprefine]
 
-[sheardispnoslip.launcher] # WARNING: THIS IS NOT PORTABLE
+[sheardispnosliprefine.launcher] # WARNING: THIS IS NOT PORTABLE
 command = mpirun -np ${nodes}
 
 # ----------------------------------------------------------------------
 # journal
 # ----------------------------------------------------------------------
-[sheardispnoslip.journal.info]
+[sheardispnosliprefine.journal.info]
 #timedependent = 1
 #implicit = 1
 #petsc = 1
@@ -21,35 +21,36 @@ command = mpirun -np ${nodes}
 # ----------------------------------------------------------------------
 # mesh_generator
 # ----------------------------------------------------------------------
-[sheardispnoslip.mesh_generator]
+[sheardispnosliprefine.mesh_generator]
 reader = pylith.meshio.MeshIOCubit
 reorder_mesh = True
+refiner = pylith.topology.RefineUniform
 
-[sheardispnoslip.mesh_generator.reader]
+[sheardispnosliprefine.mesh_generator.reader]
 filename = mesh.exo
 coordsys.space_dim = 3
 
 # ----------------------------------------------------------------------
 # problem
 # ----------------------------------------------------------------------
-[sheardispnoslip.timedependent]
+[sheardispnosliprefine.timedependent]
 dimension = 3
 bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
 
-normalizer.length_scale = 5.0*km
+normalizer.length_scale = 2.5*km
 
-[sheardispnoslip.timedependent.formulation.time_step]
+[sheardispnosliprefine.timedependent.formulation.time_step]
 total_time = 0.0*s
 
 # ----------------------------------------------------------------------
 # materials
 # ----------------------------------------------------------------------
-[sheardispnoslip.timedependent]
+[sheardispnosliprefine.timedependent]
 materials = [elastic,viscoelastic]
 materials.elastic = pylith.materials.ElasticIsotropic3D
 materials.viscoelastic = pylith.materials.ElasticIsotropic3D
 
-[sheardispnoslip.timedependent.materials.elastic]
+[sheardispnosliprefine.timedependent.materials.elastic]
 label = Elastic material
 id = 1
 db_properties.label = Elastic properties
@@ -57,7 +58,7 @@ db_properties.iohandler.filename = matprops.spatialdb
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 3
 
-[sheardispnoslip.timedependent.materials.viscoelastic]
+[sheardispnosliprefine.timedependent.materials.viscoelastic]
 label = Elastic material
 id = 2
 db_properties.label = Elastic properties
@@ -68,35 +69,35 @@ quadrature.cell.dimension = 3
 # ----------------------------------------------------------------------
 # boundary conditions
 # ----------------------------------------------------------------------
-[sheardispnoslip.timedependent.bc.x_pos]
+[sheardispnosliprefine.timedependent.bc.x_pos]
 bc_dof = [1]
 label = face_xpos
 db_initial = spatialdata.spatialdb.SimpleDB
 db_initial.label = Dirichlet BC +x edge
 db_initial.iohandler.filename = shear_dispy.spatialdb
 
-[sheardispnoslip.timedependent.bc.x_neg]
+[sheardispnosliprefine.timedependent.bc.x_neg]
 bc_dof = [1]
 label = face_xneg
 db_initial = spatialdata.spatialdb.SimpleDB
 db_initial.label = Dirichlet BC -x edge
 db_initial.iohandler.filename = shear_dispy.spatialdb
 
-[sheardispnoslip.timedependent.bc.y_pos]
+[sheardispnosliprefine.timedependent.bc.y_pos]
 bc_dof = [0]
 label = face_ypos
 db_initial = spatialdata.spatialdb.SimpleDB
 db_initial.label = Dirichlet BC +y edge
 db_initial.iohandler.filename = shear_dispx.spatialdb
 
-[sheardispnoslip.timedependent.bc.y_neg]
+[sheardispnosliprefine.timedependent.bc.y_neg]
 bc_dof = [0]
 label = face_yneg
 db_initial = spatialdata.spatialdb.SimpleDB
 db_initial.label = Dirichlet BC -y edge
 db_initial.iohandler.filename = shear_dispx.spatialdb
 
-[sheardispnoslip.timedependent.bc.z_neg]
+[sheardispnosliprefine.timedependent.bc.z_neg]
 bc_dof = [2]
 label = face_zneg
 db_initial = spatialdata.spatialdb.SimpleDB
@@ -106,16 +107,16 @@ db_initial.iohandler.filename = shear_dispz.spatialdb
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------
-[sheardispnoslip.timedependent]
+[sheardispnosliprefine.timedependent]
 interfaces = [fault]
 
-[sheardispnoslip.timedependent.interfaces.fault]
+[sheardispnosliprefine.timedependent.interfaces.fault]
 id = 100
 label = fault_x_trim
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.dimension = 2
 
-[sheardispnoslip.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+[sheardispnosliprefine.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
 slip = spatialdata.spatialdb.UniformDB
 slip.label = Final slip
 slip.values = [left-lateral-slip,reverse-slip,fault-opening]
@@ -129,7 +130,7 @@ slip_time.data = [0.0*s]
 # ----------------------------------------------------------------------
 # PETSc
 # ----------------------------------------------------------------------
-[sheardispnoslip.petsc]
+[sheardispnosliprefine.petsc]
 pc_type = asm
 
 # Change the preconditioner settings.
@@ -150,20 +151,20 @@ ksp_gmres_restart = 50
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-[sheardispnoslip.problem.formulation.output.output]
+[sheardispnosliprefine.problem.formulation.output.output]
 writer = pylith.meshio.DataWriterHDF5
-writer.filename = sheardispnoslip.h5
+writer.filename = sheardispnosliprefine.h5
 
-[sheardispnoslip.timedependent.interfaces.fault.output]
+[sheardispnosliprefine.timedependent.interfaces.fault.output]
 writer = pylith.meshio.DataWriterHDF5
-writer.filename = sheardispnoslip-fault.h5
+writer.filename = sheardispnosliprefine-fault.h5
 
-[sheardispnoslip.timedependent.materials.elastic.output]
+[sheardispnosliprefine.timedependent.materials.elastic.output]
 cell_filter = pylith.meshio.CellFilterAvg
 writer = pylith.meshio.DataWriterHDF5
-writer.filename = sheardispnoslip-elastic.h5
+writer.filename = sheardispnosliprefine-elastic.h5
 
-[sheardispnoslip.timedependent.materials.viscoelastic.output]
+[sheardispnosliprefine.timedependent.materials.viscoelastic.output]
 cell_filter = pylith.meshio.CellFilterAvg
 writer = pylith.meshio.DataWriterHDF5
-writer.filename = sheardispnoslip-viscoelastic.h5
+writer.filename = sheardispnosliprefine-viscoelastic.h5

