diff --git a/test/gw/detector/interferometer_test.py b/test/gw/detector/interferometer_test.py index 3eea3c59d4667bd3fad5232bfb866ca5b98b636d..ee575317a0c0f075bff6cbfe04e83711c7c2a038 100644 --- a/test/gw/detector/interferometer_test.py +++ b/test/gw/detector/interferometer_test.py @@ -563,6 +563,25 @@ class TestInterferometerAntennaPatternAgainstLAL(unittest.TestCase): with self.subTest(':'.join((ifo_name, pol))): self.assertAlmostEqual(std[m], 0.0, places=7) + def test_time_delay_vs_lal(self): + delays = np.zeros(self.trial) + + for n, ifo_name in enumerate(self.ifo_names): + ifo = self.ifos[n] + det = lal.cached_detector_by_prefix[self.lal_prefixes[ifo_name]] + for i in range(self.trial): + gpstime = np.random.uniform(1205303144, 1405303144) + ra = 2. * np.pi * np.random.uniform() + dec = np.pi * np.random.uniform() - np.pi / 2. + delays[i] = ( + lal.TimeDelayFromEarthCenter(det.location, ra, dec, gpstime) + - ifo.time_delay_from_geocenter(ra, dec, gpstime) + ) + + std = max(abs(delays)) + with self.subTest(ifo_name): + self.assertAlmostEqual(std, 0.0, places=10) + if __name__ == "__main__": unittest.main()