def METHOD_NAME(self):
    chisq_quick = {}
    chisq_full = {}
    chisq_ref = {}
    for ifo in self.ifos:
        nbins = 26
        dof = nbins * 2 - 2
        bins = power_chisq_bins(self.hp, nbins, self.psd[ifo],
                                low_frequency_cutoff=20.0)
        chisq = power_chisq_at_points_from_precomputed \
            (self.corr[ifo], self.snr_unnorm[ifo][27402:27492].data,
             self.norm[ifo], bins, indices=numpy.arange(27402,27492))
        chisq_quick[ifo] = chisq /  dof
        chisq = power_chisq(self.hp, self.data[ifo], nbins,
                            self.psd[ifo], low_frequency_cutoff=20.0)
        chisq_full[ifo] = chisq[27402:27492] / dof
        max_diff = max(abs(chisq_full[ifo] - chisq_quick[ifo]))
        self.assertTrue(max_diff < 1E-5)