def test_get_fault_vertices_3d(self):
p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
Line([Point(10., 45.2),
Point(10.5, 45.3),
Point(10., 45.487783)]),
upper_seismogenic_depth=0.,
lower_seismogenic_depth=14.,
dip=30,
index_patch=1.)
p4, p5, p6, p7 = SimpleFaultSurface.get_fault_patch_vertices(
Line([Point(10., 45.2),
Point(10.5, 45.3),
Point(10., 45.487783)]),
upper_seismogenic_depth=0.,
lower_seismogenic_depth=14.,
dip=30,
index_patch=2.)
# Test for the first patch
self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
self.assertAlmostEqual(p0.latitude, 45.2, delta=0.1)
self.assertAlmostEqual(p0.depth, 0., delta=0.1)
self.assertAlmostEqual(p1.longitude, 10.5, delta=0.1)
self.assertAlmostEqual(p1.latitude, 45.3, delta=0.1)
self.assertAlmostEqual(p1.depth, 0., delta=0.1)
self.assertAlmostEqual(p2.longitude, 10.81, delta=0.1)
self.assertAlmostEqual(p2.latitude, 45.2995, delta=0.1)
self.assertAlmostEqual(p2.depth, 14., delta=0.1)
self.assertAlmostEqual(p3.longitude, 10.31, delta=0.1)
self.assertAlmostEqual(p3.latitude, 45.1995, delta=0.1)
self.assertAlmostEqual(p3.depth, 14., delta=0.1)
# Test for the second patch
self.assertAlmostEqual(p4.longitude, 10.5, delta=0.1)
self.assertAlmostEqual(p4.latitude, 45.3, delta=0.1)
self.assertAlmostEqual(p4.depth, 0., delta=0.1)
self.assertAlmostEqual(p5.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p5.latitude, 45.4877, delta=0.1)
self.assertAlmostEqual(p5.depth, 0., delta=0.1)
self.assertAlmostEqual(p6.longitude, 10.311, delta=0.1)
self.assertAlmostEqual(p6.latitude, 45.4873, delta=0.1)
self.assertAlmostEqual(p6.depth, 14., delta=0.1)
self.assertAlmostEqual(p7.longitude, 10.81, delta=0.1)
self.assertAlmostEqual(p7.latitude, 45.2995, delta=0.1)
self.assertAlmostEqual(p7.depth, 14., delta=0.1)
def test_get_fault_vertices_3d(self):
p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
Line([Point(10., 45.2), Point(10.5, 45.3), Point(10., 45.487783)]),
upper_seismogenic_depth=0., lower_seismogenic_depth=14.,
dip=30, index_patch=1.
)
p4, p5, p6, p7 = SimpleFaultSurface.get_fault_patch_vertices(
Line([Point(10., 45.2), Point(10.5, 45.3), Point(10., 45.487783)]),
upper_seismogenic_depth=0., lower_seismogenic_depth=14.,
dip=30, index_patch=2.
)
# Test for the first patch
self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
self.assertAlmostEqual(p0.latitude, 45.2, delta=0.1)
self.assertAlmostEqual(p0.depth, 0., delta=0.1)
self.assertAlmostEqual(p1.longitude, 10.5, delta=0.1)
self.assertAlmostEqual(p1.latitude, 45.3, delta=0.1)
self.assertAlmostEqual(p1.depth, 0., delta=0.1)
self.assertAlmostEqual(p2.longitude, 10.81, delta=0.1)
self.assertAlmostEqual(p2.latitude, 45.2995, delta=0.1)
self.assertAlmostEqual(p2.depth, 14., delta=0.1)
self.assertAlmostEqual(p3.longitude, 10.31, delta=0.1)
self.assertAlmostEqual(p3.latitude, 45.1995, delta=0.1)
self.assertAlmostEqual(p3.depth, 14., delta=0.1)
# Test for the second patch
self.assertAlmostEqual(p4.longitude, 10.5, delta=0.1)
self.assertAlmostEqual(p4.latitude, 45.3, delta=0.1)
self.assertAlmostEqual(p4.depth, 0., delta=0.1)
self.assertAlmostEqual(p5.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p5.latitude, 45.4877, delta=0.1)
self.assertAlmostEqual(p5.depth, 0., delta=0.1)
self.assertAlmostEqual(p6.longitude, 10.311, delta=0.1)
self.assertAlmostEqual(p6.latitude, 45.4873, delta=0.1)
self.assertAlmostEqual(p6.depth, 14., delta=0.1)
self.assertAlmostEqual(p7.longitude, 10.81, delta=0.1)
self.assertAlmostEqual(p7.latitude, 45.2995, delta=0.1)
self.assertAlmostEqual(p7.depth, 14., delta=0.1)
def test_get_fault_vertices_3d_upper_seismogenic_depth_with_depth(self):
p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
Line([Point(10., 45.0), Point(10., 45.5)]),
upper_seismogenic_depth=10., lower_seismogenic_depth=14.,
dip=90, index_patch=1.
)
# The value used for this test is by hand calculation
self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
self.assertAlmostEqual(p0.latitude, 45.0, delta=0.1)
self.assertAlmostEqual(p0.depth, 10., delta=0.1)
self.assertAlmostEqual(p1.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p1.latitude, 45.5, delta=0.1)
self.assertAlmostEqual(p1.depth, 10., delta=0.1)
self.assertAlmostEqual(p2.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p2.latitude, 45.5, delta=0.1)
self.assertAlmostEqual(p2.depth, 14., delta=0.1)
self.assertAlmostEqual(p3.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p3.latitude, 45.0, delta=0.1)
self.assertAlmostEqual(p3.depth, 14., delta=0.1)
def test_get_fault_vertices_3d_upper_seismogenic_depth_with_depth(self):
p0, p1, p2, p3 = SimpleFaultSurface.get_fault_patch_vertices(
Line([Point(10., 45.0), Point(10., 45.5)]),
upper_seismogenic_depth=10., lower_seismogenic_depth=14.,
dip=90, index_patch=1.
)
# The value used for this test is by hand calculation
self.assertAlmostEqual(p0.longitude, 10., delta=0.1)
self.assertAlmostEqual(p0.latitude, 45.0, delta=0.1)
self.assertAlmostEqual(p0.depth, 10., delta=0.1)
self.assertAlmostEqual(p1.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p1.latitude, 45.5, delta=0.1)
self.assertAlmostEqual(p1.depth, 10., delta=0.1)
self.assertAlmostEqual(p2.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p2.latitude, 45.5, delta=0.1)
self.assertAlmostEqual(p2.depth, 14., delta=0.1)
self.assertAlmostEqual(p3.longitude, 10.0, delta=0.1)
self.assertAlmostEqual(p3.latitude, 45.0, delta=0.1)
self.assertAlmostEqual(p3.depth, 14., delta=0.1)