From 43af1c174d43ab9d93826aa56c38118c704f570b Mon Sep 17 00:00:00 2001 From: James Hensman Date: Sun, 16 Dec 2012 17:31:51 +0000 Subject: [PATCH] some simplification of the gradient expressions --- GPy/models/uncollapsed_sparse_GP.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/GPy/models/uncollapsed_sparse_GP.py b/GPy/models/uncollapsed_sparse_GP.py index f0193a31..fb728e7b 100644 --- a/GPy/models/uncollapsed_sparse_GP.py +++ b/GPy/models/uncollapsed_sparse_GP.py @@ -61,8 +61,10 @@ class uncollapsed_sparse_GP(sparse_GP_regression): self.dL_dpsi2 = 0.5 * self.beta * self.D * (self.Kmmi - mdot(self.Kmmi,self.q_u_expectation[1],self.Kmmi)) # Compute dL_dKmm - tmp = -0.5*self.beta*self.D*self.psi2 -0.5*self.D*self.Kmm +0.5*self.D*(self.q_u_expectation[1]) + 0.5*self.beta*mdot(self.psi2,self.Kmmi,self.q_u_expectation[1]) + 0.5*self.beta*mdot(self.q_u_expectation[1],self.Kmmi, self.psi2) - tmptmp = - np.dot(self.q_u_expectation[0],self.psi1V.T)*0.5 + tmp = 0.5*self.beta*mdot(self.psi2,self.Kmmi,self.q_u_expectation[1]) + tmp += tmp.T + tmp += 0.5*self.D*(-self.beta*self.psi2 - self.Kmm + self.q_u_expectation[1]) + tmptmp = - 0.5*np.dot(self.q_u_expectation[0],self.psi1V.T) tmp += tmptmp + tmptmp.T self.dL_dKmm = mdot(self.Kmmi,tmp,self.Kmmi)