apunkt/GPy
1
0
Fork 0
mirror of https://github.com/SheffieldML/GPy.git synced 2026-05-14 14:32:37 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

oops, that was silly bug. Don't code tired, kids.

Browse source 
problem summing psi2 in sparse_GP
This commit is contained in:
James Hensman 2013-07-30 11:43:07 +01:00
parent b9c9f5ac83
commit 632dda8397
1 changed files with 2 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

7
GPy/core/sparse_gp.py
View file

@ -80,8 +80,6 @@ class SparseGP(GPBase):
self.psi2 = None self.psi2 = None
def _computations(self): def _computations(self):
# factor Kmm # factor Kmm
self.Lm = jitchol(self.Kmm) self.Lm = jitchol(self.Kmm)
@ -90,10 +88,9 @@ class SparseGP(GPBase):
if self.likelihood.is_heteroscedastic: if self.likelihood.is_heteroscedastic:
psi2_beta = (self.psi2 * (self.likelihood.precision.flatten().reshape(self.num_data, 1, 1))).sum(0) psi2_beta = (self.psi2 * (self.likelihood.precision.flatten().reshape(self.num_data, 1, 1))).sum(0)
else: else:
#psi2_beta = self.psi2.sum(0) * self.likelihood.precision psi2_beta = self.psi2.sum(0) * self.likelihood.precision
psi2_beta = self.psi2[0] * (self.num_data * self.likelihood.precision)
evals, evecs = linalg.eigh(psi2_beta) evals, evecs = linalg.eigh(psi2_beta)
clipped_evals = np.clip(evals, 0., 1e16) # TODO: make clipping configurable clipped_evals = np.clip(evals, 0., 1e6) # TODO: make clipping configurable
if not np.array_equal(evals, clipped_evals): if not np.array_equal(evals, clipped_evals):
pass#print evals pass#print evals
tmp = evecs * np.sqrt(clipped_evals) tmp = evecs * np.sqrt(clipped_evals)