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)