[cig-commits] [commit] knepley/upgrade-petsc-interface: Added full-scale test for embedded fault + refinement. (1c6555b)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue Mar 11 13:31:24 PDT 2014
Repository : ssh://geoshell/pylith
On branch : knepley/upgrade-petsc-interface
Link : https://github.com/geodynamics/pylith/compare/096d4cbeaf5370d43b58aa6c5bd55ed51e51d3ec...1c6555be39a79a031be6ee0949038af870fc03cb
>---------------------------------------------------------------
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]
More information about the CIG-COMMITS
mailing list