

def loglkl(self, cosmo, data):
    difference = []
    for (z, Dv) in zip(self.z_eff, self.Dv):
        da = cosmo.angular_distance(z)
        dr = (rs / cosmo.Hubble(z))
        dv = pow((((da ** 2) * ((1.0 + z) ** 2)) * dr), (1.0 / 3))
        rs = cosmo.rs_drag()
        difference.append((Dv - ((dv / rs) * self.rs_fiducial)))
    difference = np.array(difference)
    likelihood = np.exp(((- 0.5) * np.dot(difference, np.dot(self.inverse_covmat, difference.T))))
    return log(likelihood)

dataset/ETHPy150Open baudren/montepython_public/montepython/likelihoods/WiggleZ_bao/__init__.py WiggleZ_bao.loglkl/VarMisuse@51/63 `z`->`rs`