From 13438f1c10dbabca6c4218a23364e13d29a4685b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?David=20Men=C3=A9ndez=20Hurtado?=
 <david.menendez.hurtado@scilifelab.se>
Date: Mon, 24 Aug 2015 13:46:20 +0200
Subject: [PATCH] Fixing the blas arguments for DSYRK

---
 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 4c5f1b79..acceeb16 100644
--- a/GPy/util/linalg.py
+++ b/GPy/util/linalg.py
@@ -311,7 +311,7 @@ def tdot_blas(mat, out=None):
     # # Call to DSYRK from BLAS
     mat = np.asfortranarray(mat)
     out = blas.dsyrk(alpha=1.0, a=mat, beta=0.0, c=out, overwrite_c=1,
-                     trans=0, lower=1, trans=0)
+                     trans=0, lower=0)
 
     symmetrify(out, upper=True)
     return np.ascontiguousarray(out)
@@ -329,7 +329,7 @@ def DSYR_blas(A, x, alpha=1.):
     :param alpha: scalar
 
     """
-    blas.dsyr(lower=1, x=x, a=A, alpha=alpha)
+    A = blas.dsyr(lower=0, x=x, a=A, alpha=alpha, overwrite_a=True)
     symmetrify(A, upper=True)
 
 def DSYR_numpy(A, x, alpha=1.):