def _ppcf(self, samples): nrvs = norm.ppf(samples) vals = norm.cdf(nrvs[:, 0] * np.sqrt(1 - self.theta**2) + self.theta * nrvs[:, 1]) return vals