From d20222a90bac8181f39fac50f851d048e91c5bac Mon Sep 17 00:00:00 2001
From: James Hensman <james.hensman@gmail.com>
Date: Wed, 4 Jun 2014 15:43:03 +0100
Subject: [PATCH] reverting Maxs linalg changes

---
 GPy/util/linalg.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/GPy/util/linalg.py b/GPy/util/linalg.py
index 6a38a7c3..b6b4a204 100644
--- a/GPy/util/linalg.py
+++ b/GPy/util/linalg.py
@@ -151,7 +151,7 @@ def dpotri(A, lower=1):
     assert lower==1, "scipy linalg behaviour is very weird. please use lower, fortran ordered arrays"
 
     A = force_F_ordered(A)
-    R, info = lapack.dpotri(A, lower=1)
+    R, info = lapack.dpotri(A, lower=0) #needs to be zero here, seems to be a scipy bug
     symmetrify(R)
     return R, info
 
@@ -218,7 +218,7 @@ def pdinv(A, *args):
     L = jitchol(A, *args)
     logdet = 2.*np.sum(np.log(np.diag(L)))
     Li = dtrtri(L)
-    Ai, _ = lapack.dpotri(L, lower=1)
+    Ai, _ = dpotri(L, lower=1)
     # Ai = np.tril(Ai) + np.tril(Ai,-1).T
     symmetrify(Ai)