def test_xcorr_3c(self):
st = read()
st2 = read()
for tr in st2:
tr.data = -5 * np.roll(tr.data, 50)
shift, value, x = xcorr_3c(st, st2, 200, full_xcorr=True)
self.assertEqual(shift, -50)
self.assertAlmostEqual(value, -0.998, 3)
def test_xcorr_3c(self):
st = read()
st2 = read()
for tr in st2:
tr.data = -5 * np.roll(tr.data, 50)
shift, value, x = xcorr_3c(st, st2, 200, full_xcorr=True)
self.assertEqual(shift, -50)
self.assertAlmostEqual(value, -0.998, 3)
def test_xcorr_3c(self):
st = read()
st2 = read()
for tr in st2:
tr.data = -5 * np.roll(tr.data, 50)
shift, value, x = xcorr_3c(st, st2, 200, full_xcorr=True)
assert shift == -50
assert round(abs(value - -0.998), 3) == 0