def METHOD_NAME(self):
    vals = list(zip(self.ra, self.dec, self.pol, self.time))
    for ifo in self.d:
        fp = []
        fc = []
        for ra1, dec1, pol1, time1 in vals:
            gmst = lal.GreenwichMeanSiderealTime(time1)
            fp1, fc1 = tuple(lal.ComputeDetAMResponse(ifo.response, ra1, dec1, pol1, gmst))
            fp.append(fp1)
            fc.append(fc1)
        fp2, fc2 = ifo.antenna_pattern(self.ra, self.dec, self.pol, self.time)
        fp = numpy.array(fp)
        fc = numpy.array(fc)
        diff1 = fp - fp2
        diff2 = fc - fc2
        diff = abs(numpy.concatenate([diff1, diff2]))
        tolerance = 2e-4
        print("Max antenna diff:", ifo.name, diff.max())
        self.assertLess(diff.max(), tolerance)