def setUp(self):
"""
Setup for test.
"""
TestTet4.setUp(self)
self.mesh = {
'ncells-elastic': 1266 * 8,
'ncells-viscoelastic': 1276 * 8,
'ncorners': 4,
'nvertices': 4016,
'spaceDim': 3,
'tensorSize': 6
}
self.nverticesO = self.mesh['nvertices']
self.mesh['nvertices'] += 36
self.faultMesh = {
'nvertices': 57,
'spaceDim': 3,
'ncells': 80,
'ncorners': 3
}
run_pylith(ShearApp, GenerateDB, nprocs=4)
self.outputRoot = "sheardispnosliprefine"
self.soln = AnalyticalSoln()
return
def setUp(self):
"""
Setup for test.
"""
TestHex8.setUp(self)
self.mesh = {
'ncells-elastic': 180 * 8,
'ncells-viscoelastic': 180 * 8,
'ncorners': 8,
'nvertices': 3591,
'spaceDim': 3,
'tensorSize': 6
}
self.nverticesO = self.mesh['nvertices']
self.mesh['nvertices'] += 44
self.faultMesh = {
'nvertices': 65,
'spaceDim': 3,
'ncells': 48,
'ncorners': 4
}
run_pylith(ShearApp, GenerateDB, nprocs=4)
self.outputRoot = "sheardispnosliprefine"
self.soln = AnalyticalSoln()
return
def setUp(self):
"""
Setup for test.
"""
TestTet4.setUp(self)
self.nverticesO = self.mesh['nvertices']
# Fault x
self.mesh['nvertices'] += 8
self.faultMeshX = {
'nvertices': 19,
'spaceDim': 3,
'ncells': 20,
'ncorners': 3
}
# Fault y
self.mesh['nvertices'] += 2
self.faultMeshY = {
'nvertices': 9,
'spaceDim': 3,
'ncells': 8,
'ncorners': 3
}
run_pylith(FaultsIntersectNoSlipApp, GenerateDB)
self.outputRoot = "faultsintersectnoslip"
self.soln = AnalyticalSoln()
return
def setUp(self):
"""
Setup for test.
"""
TestTri3.setUp(self)
self.mesh = {
'ncells': 142 * 4,
'ncorners': 3,
'nvertices': 317,
'spaceDim': 2,
'tensorSize': 3
}
self.nverticesO = self.mesh['nvertices']
self.mesh['nvertices'] += 3
self.faultMesh = {
'nvertices': 5,
'spaceDim': 2,
'ncells': 2 * 2,
'ncorners': 2
}
run_pylith(ShearApp, GenerateDB, nprocs=3)
self.outputRoot = "sheardispnosliprefine"
self.soln = AnalyticalSoln()
return
def setUp(self):
"""
Setup for test.
"""
TestTri3.setUp(self)
run_pylith(ShearApp, GenerateDB)
self.outputRoot = "sheardisp"
self.soln = AnalyticalSoln()
return
def setUp(self):
"""
Setup for test.
"""
TestTet4.setUp(self)
self.nverticesO = self.mesh['nvertices']
self.mesh['nvertices'] += 8
self.faultMesh = {'nvertices': 19,
'spaceDim': 3,
'ncells': 20,
'ncorners': 3}
run_pylith(ShearApp, GenerateDB, nprocs=3)
self.outputRoot = "sheardispnoslip"
self.soln = AnalyticalSoln()
return
def run(self):
"""
Generate the database.
"""
# Domain
x = numpy.arange(-4000.0, 4000.1, 1000.0)
y = numpy.arange(-4000.0, 4000.1, 1000.0)
npts = x.shape[0]
xx = x * numpy.ones((npts, 1), dtype=numpy.float64)
yy = y * numpy.ones((npts, 1), dtype=numpy.float64)
xy = numpy.zeros((npts**2, 2), dtype=numpy.float64)
xy[:, 0] = numpy.ravel(xx)
xy[:, 1] = numpy.ravel(numpy.transpose(yy))
from sheardisp_soln import AnalyticalSoln
soln = AnalyticalSoln()
disp = soln.displacement(xy)
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
data = {
'points':
xy,
'coordsys':
cs,
'data_dim':
2,
'values': [{
'name': "displacement-x",
'units': "m",
'data': numpy.ravel(disp[0, :, 0])
}, {
'name': "displacement-y",
'units': "m",
'data': numpy.ravel(disp[0, :, 1])
}]
}
from spatialdata.spatialdb.SimpleIOAscii import SimpleIOAscii
io = SimpleIOAscii()
io.inventory.filename = "shear_disp.spatialdb"
io._configure()
io.write(data)
return
def setUp(self):
"""
Setup for test.
"""
TestHex8.setUp(self)
self.nverticesO = self.mesh['nvertices']
self.mesh['nvertices'] += 10
self.faultMesh = {'nvertices': 21,
'spaceDim': 3,
'ncells': 12,
'ncorners': 4}
run_pylith(ShearApp, GenerateDB, nprocs=3)
self.outputRoot = "frictionnoslip"
self.soln = AnalyticalSoln()
return
