From 38ed60385aa0c9a2c1b54150e6efcf0fd5350620 Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Tue, 3 Jun 2014 15:42:32 -0700 Subject: [PATCH] linalg had lowers missing for windows libraries to work correctly --- GPy/util/linalg.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/GPy/util/linalg.py b/GPy/util/linalg.py index 661a2b47..6a38a7c3 100644 --- a/GPy/util/linalg.py +++ b/GPy/util/linalg.py @@ -30,6 +30,7 @@ if config.getboolean('anaconda', 'installed') and config.getboolean('anaconda', dsyrk = mkl_rt.dsyrk dsyr = mkl_rt.dsyr _blas_available = True + print 'anaconda installed and mkl is loaded' except: _blas_available = False else: @@ -150,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=0) + R, info = lapack.dpotri(A, lower=1) symmetrify(R) return R, info @@ -217,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) + Ai, _ = lapack.dpotri(L, lower=1) # Ai = np.tril(Ai) + np.tril(Ai,-1).T symmetrify(Ai)