From 61a79c504196e7fc1a9ae5dadbccbe42a1200935 Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Thu, 16 May 2013 13:47:33 +0100 Subject: [PATCH] clipping bounds released to 1e6 --- GPy/models/sparse_GP.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/GPy/models/sparse_GP.py b/GPy/models/sparse_GP.py index ed5a5004..2aafab16 100644 --- a/GPy/models/sparse_GP.py +++ b/GPy/models/sparse_GP.py @@ -76,7 +76,7 @@ class sparse_GP(GP): # psi2_beta_scaled = (self.psi2 * (self.likelihood.precision.flatten().reshape(self.N, 1, 1) / sf2)).sum(0) psi2_beta_scaled = (self.psi2 * (self.likelihood.precision.flatten().reshape(self.N, 1, 1))).sum(0) evals, evecs = linalg.eigh(psi2_beta_scaled) - clipped_evals = np.clip(evals, 0., 1e15) # TODO: make clipping configurable + clipped_evals = np.clip(evals, 0., 1e6) # TODO: make clipping configurable if not np.allclose(evals, clipped_evals): print "Warning: clipping posterior eigenvalues" tmp = evecs * np.sqrt(clipped_evals)