some simplification of the gradient expressions

This commit is contained in:
James Hensman 2012-12-16 17:31:51 +00:00
parent 7a34e1c446
commit 43af1c174d

View file

@ -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)) self.dL_dpsi2 = 0.5 * self.beta * self.D * (self.Kmmi - mdot(self.Kmmi,self.q_u_expectation[1],self.Kmmi))
# Compute dL_dKmm # 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) tmp = 0.5*self.beta*mdot(self.psi2,self.Kmmi,self.q_u_expectation[1])
tmptmp = - np.dot(self.q_u_expectation[0],self.psi1V.T)*0.5 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 tmp += tmptmp + tmptmp.T
self.dL_dKmm = mdot(self.Kmmi,tmp,self.Kmmi) self.dL_dKmm = mdot(self.Kmmi,tmp,self.Kmmi)