def test_akirichards_alt(self):
vp1 = 12250.
vp2 = 11600.
vs1 = 6200.
vs2 = 6650.
rho1 = 2.66
rho2 = 2.34
theta = 40.
reflect = avo.akirichards_alt(vp1, vs1, rho1, vp2, vs2, rho2, theta)
self.assertAlmostEquals(reflect, -0.112236, places=3)
def test_akirichards_alt(self):
vp1 = 12250.
vp2 = 11600.
vs1 = 6200.
vs2 = 6650.
rho1 = 2.66
rho2 = 2.34
theta = 40.
reflect = avo.akirichards_alt( vp1,vs1,rho1,vp2,
vs2,rho2,theta )
self.assertAlmostEquals( reflect, -0.112236, places = 3 )
def test_akirichards_alt(self):
vp1 = 12250.
vp2 = 11600.
vs1 = 6200.
vs2 = 6650.
rho1 = 2.66
rho2 = 2.34
theta = np.arange(45)
reflect = avo.akirichards_alt(vp1, vs1, rho1, vp2, vs2, rho2, theta)
reflect_zoep = avo.zoeppritz(vp1, vs1, rho1, vp2, vs2, rho2, theta)
# See if it is within .1 of zoep for < 45 deg
test = np.allclose(reflect, reflect_zoep, rtol=self.tolerance)
self.assertTrue(test)
def akirichards_alt(Rp0, Rp1, theta1):
reflection.akirichards_alt(Rp0.vp, Rp0.vs, Rp0.rho, Rp1.vp, Rp1.vs, Rp1.rho, theta1)