From cb1ab89d8a5eebca414284ef35d635ccafffe516 Mon Sep 17 00:00:00 2001 From: Neil Lawrence Date: Sun, 31 Dec 2017 22:54:43 +0100 Subject: [PATCH] Update gp.py Sample return seemed to have been based on number of training data, not number of posterior samples requested. --- GPy/core/gp.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/GPy/core/gp.py b/GPy/core/gp.py index 7bad7648..4481d3bc 100644 --- a/GPy/core/gp.py +++ b/GPy/core/gp.py @@ -564,7 +564,7 @@ class GP(Model): if self.output_dim == 1: return sim_one_dim(m, v) else: - fsim = np.empty((self.output_dim, self.num_data, size)) + fsim = np.empty((self.output_dim, X.shape[1], size)) for d in range(self.output_dim): if full_cov and v.ndim == 3: fsim[d] = sim_one_dim(m[:, d], v[:, :, d])