92384bb371d2ee8a3ffbacbf603acd5328855464,pymc3/distributions/multivariate.py,LKJCholeskyCov,logp,#LKJCholeskyCov#Any#,566

Before Change

        corr_diag = x[diag_idxs] / rowlengths
        corr_logdet = np.log(corr_diag).sum()

        det_invjac = np.log(0.5 * rowlengths).sum()

        return (self.n - 1) * corr_logdet + logp_sd + det_invjac

After Change

        corr_diag = x[diag_idxs] / sd_vals
        corr_logdet = np.log(corr_diag).sum()

        count = np.arange(self.n - 1)
        det_invjac = - (count * tt.log(sd_vals[1:])).sum()
        det_invjac += - tt.log(x[diag_idxs]).sum() + tt.log(x[0])

        return (self.n - 1) * corr_logdet + logp_sd + det_invjac

Instances