From 43b720c848ad9a8f76a7334c349d88ff7d2f11aa Mon Sep 17 00:00:00 2001 From: James Hensman Date: Thu, 25 Apr 2013 16:01:36 +0100 Subject: [PATCH] more minor simplifications --- GPy/models/sparse_GP.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/GPy/models/sparse_GP.py b/GPy/models/sparse_GP.py index 56a764af..e158e026 100644 --- a/GPy/models/sparse_GP.py +++ b/GPy/models/sparse_GP.py @@ -109,8 +109,10 @@ class sparse_GP(GP): self.psi1V = np.dot(self.psi1, self.V) #tmp = np.dot(self.Lmi.T, self.LBi.T) - tmp = linalg.lapack.clapack.dtrtrs(self.Lm.T,np.asarray(self.LBi.T,order='C'),lower=0)[0] - self.C = np.dot(tmp,tmp.T) #TODO: tmp is triangular. replace with dtrmm (blas) when available + #tmp = linalg.lapack.clapack.dtrtrs(self.Lm.T,np.asarray(self.LBi.T,order='C'),lower=0)[0] + #self.C = np.dot(tmp,tmp.T) #TODO: tmp is triangular. replace with dtrmm (blas) when available + tmp = linalg.lapack.flapack.dtrtrs(self.Lm,np.asfortranarray(self.Bi),lower=1,trans=1)[0] + self.C = linalg.lapack.flapack.dtrtrs(self.Lm,np.asfortranarray(tmp.T),lower=1,trans=1)[0] self.Cpsi1V = np.dot(self.C,self.psi1V) self.Cpsi1VVpsi1 = np.dot(self.Cpsi1V,self.psi1V.T) #self.E = np.dot(self.Cpsi1VVpsi1,self.C)/sf2