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