def sample(self, k): #samples = self.distribution.rvs(k).reshape(k,p) samples = self.rng.multivariate_normal(self.mean, self.cov, k) return samples