[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 full-scale test for embedded fault + refinement. (1c6555b)

[cig_noreply at geodynamics.org]

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 1c6555be39a79a031be6ee0949038af870fc03cb
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Tue Mar 11 13:31:17 2014 -0700

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


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

1c6555be39a79a031be6ee0949038af870fc03cb
 ...earDispNoSlip.py => TestShearDispNoSlipRefine.py} | 20 +++++++++++++-------
 .../3d/{tet4 => hex8}/sheardispnosliprefine.cfg      | 10 +++++-----
 2 files changed, 18 insertions(+), 12 deletions(-)

diff --git a/tests_auto/3d/hex8/TestShearDispNoSlip.py b/tests_auto/3d/hex8/TestShearDispNoSlipRefine.py
similarity index 87%
copy from tests_auto/3d/hex8/TestShearDispNoSlip.py
copy to tests_auto/3d/hex8/TestShearDispNoSlipRefine.py
index d36e4c1..c291965 100644
--- a/tests_auto/3d/hex8/TestShearDispNoSlip.py
+++ b/tests_auto/3d/hex8/TestShearDispNoSlipRefine.py
@@ -16,10 +16,10 @@
 # ----------------------------------------------------------------------
 #
 
-## @file tests/3d/hex8/TestShearDispNoSlip.py
+## @file tests/3d/hex8/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 TestHex8 import TestHex8
@@ -31,7 +31,7 @@ from pylith.tests.Fault import check_vertex_fields
 from pylith.apps.PyLithApp import PyLithApp
 class ShearApp(PyLithApp):
   def __init__(self):
-    PyLithApp.__init__(self, name="sheardispnoslip")
+    PyLithApp.__init__(self, name="sheardispnosliprefine")
     return
 
 
@@ -53,7 +53,7 @@ def run_pylith():
   return
 
 
-class TestShearDispNoSlip(TestHex8):
+class TestShearDispNoSlipRefine(TestHex8):
   """
   Test suite for testing pylith with 2-D shear extension.
   """
@@ -63,15 +63,21 @@ class TestShearDispNoSlip(TestHex8):
     Setup for test.
     """
     TestHex8.setUp(self)
+    self.mesh = {'ncells-elastic': 768*8,
+                 'ncells-viscoelastic': 1280*8,
+                 'ncorners': 8,
+                 'nvertices': 2601,
+                 'spaceDim': 3,
+                 'tensorSize': 6}
     self.nverticesO = self.mesh['nvertices']
     self.mesh['nvertices'] += 44
     self.faultMesh = {'nvertices': 44,
                       'spaceDim': 3,
-                      'ncells': 30,
+                      'ncells': 30*4,
                       'ncorners': 4}
 
     run_pylith()
-    self.outputRoot = "sheardispnoslip"
+    self.outputRoot = "sheardispnosliprefine"
 
     self.soln = AnalyticalSoln()
     return
@@ -169,7 +175,7 @@ class TestShearDispNoSlip(TestHex8):
 # ----------------------------------------------------------------------
 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/sheardispnosliprefine.cfg b/tests_auto/3d/hex8/sheardispnosliprefine.cfg
similarity index 96%
copy from tests_auto/3d/tet4/sheardispnosliprefine.cfg
copy to tests_auto/3d/hex8/sheardispnosliprefine.cfg
index a3084cd..a43897b 100644
--- a/tests_auto/3d/tet4/sheardispnosliprefine.cfg
+++ b/tests_auto/3d/hex8/sheardispnosliprefine.cfg
@@ -16,7 +16,7 @@ command = mpirun -np ${nodes}
 #implicitelasticity = 1
 #quadrature3d = 1
 #faultcohesivekin = 1
-#fiatsimplex = 1
+#fiatlagrange = 1
 
 # ----------------------------------------------------------------------
 # mesh_generator
@@ -37,7 +37,7 @@ coordsys.space_dim = 3
 dimension = 3
 bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
 
-normalizer.length_scale = 2.5*km
+normalizer.length_scale = 5.0*km
 
 [sheardispnosliprefine.timedependent.formulation.time_step]
 total_time = 0.0*s
@@ -55,7 +55,7 @@ label = Elastic material
 id = 1
 db_properties.label = Elastic properties
 db_properties.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
 
 [sheardispnosliprefine.timedependent.materials.viscoelastic]
@@ -63,7 +63,7 @@ label = Elastic material
 id = 2
 db_properties.label = Elastic properties
 db_properties.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
 
 # ----------------------------------------------------------------------
@@ -113,7 +113,7 @@ interfaces = [fault]
 [sheardispnosliprefine.timedependent.interfaces.fault]
 id = 100
 label = fault_x_trim
-quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
 [sheardispnosliprefine.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]