class TestFaultXYP(TestTri3): """ Test suite for testing sense of slip. """ def setUp(self): """ Setup for test. """ TestTri3.setUp(self) run_pylith(FaultXYPApp, GenDBFaultXYP, nprocs=3) self.mesh['nvertices'] += 23 self.outputRoot = "faultxyp" self.soln = SolnFaultXYP() return def calcStateVar(self, name, vertices, cells): """ Calculate state variable. """ ncells = self.mesh['ncells'] pts = numpy.zeros( (ncells, 3), dtype=numpy.float64) if name == "total_strain": stateVar = self.soln.strain(pts) elif name == "stress": stateVar = self.soln.stress(pts) else: raise ValueError("Unknown state variable '%s'." % name) return stateVar
class TestFaultXYP(TestTri3): """ Test suite for testing sense of slip. """ def setUp(self): """ Setup for test. """ TestTri3.setUp(self) run_pylith(FaultXYPApp, GenDBFaultXYP, nprocs=3) self.mesh['nvertices'] += 23 self.outputRoot = "faultxyp" self.soln = SolnFaultXYP() return def calcStateVar(self, name, vertices, cells): """ Calculate state variable. """ ncells = self.mesh['ncells'] pts = numpy.zeros((ncells, 3), dtype=numpy.float64) if name == "total_strain": stateVar = self.soln.strain(pts) elif name == "stress": stateVar = self.soln.stress(pts) else: raise ValueError("Unknown state variable '%s'." % name) return stateVar