From 57a751d70df7388696dcc271b079b940a8497fa1 Mon Sep 17 00:00:00 2001 From: James Hensman Date: Fri, 15 Feb 2013 16:30:00 +0000 Subject: [PATCH] fixed bug with the bias kernel. --- GPy/models/sparse_GP.py | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/GPy/models/sparse_GP.py b/GPy/models/sparse_GP.py index 3239d462..48b88bc4 100644 --- a/GPy/models/sparse_GP.py +++ b/GPy/models/sparse_GP.py @@ -104,9 +104,11 @@ class sparse_GP(GP): self.dL_dpsi2 += - 0.5 * self.likelihood.precision[:,None,None]/sf2 * self.D * self.C[None,:,:] # dC self.dL_dpsi2 += - 0.5 * self.likelihood.precision[:,None,None]* self.E[None,:,:] # dD else: - self.dL_dpsi2 = 0.5 * self.likelihood.precision * self.D * self.Kmmi[None,:,:] # dB - self.dL_dpsi2 += - 0.5 * self.likelihood.precision/sf2 * self.D * self.C[None,:,:] # dC - self.dL_dpsi2 += - 0.5 * self.likelihood.precision * self.E[None,:,:] # dD + self.dL_dpsi2 = 0.5 * self.likelihood.precision * self.D * self.Kmmi # dB + self.dL_dpsi2 += - 0.5 * self.likelihood.precision/sf2 * self.D * self.C # dC + self.dL_dpsi2 += - 0.5 * self.likelihood.precision * self.E # dD + #repeat for each of the N psi_2 matrices + self.dL_dpsi2 = np.repeat(self.dL_dpsi2[None,:,:],self.N,axis=0) # Compute dL_dKmm self.dL_dKmm = -0.5 * self.D * mdot(self.Lmi.T, self.A, self.Lmi)*sf2 # dB