apunkt/GPy
1
0
Fork 0
mirror of https://github.com/SheffieldML/GPy.git synced 2026-05-12 05:22:38 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

Fixing the blas arguments for DSYRK

Browse source
This commit is contained in:
David Menéndez Hurtado 2015-08-24 13:46:20 +02:00
parent 3746af7729
commit 13438f1c10
1 changed files with 2 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

4
GPy/util/linalg.py
View file

@ -311,7 +311,7 @@ def tdot_blas(mat, out=None):
# # Call to DSYRK from BLAS # # Call to DSYRK from BLAS
mat = np.asfortranarray(mat) mat = np.asfortranarray(mat)
out = blas.dsyrk(alpha=1.0, a=mat, beta=0.0, c=out, overwrite_c=1, 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) symmetrify(out, upper=True)
return np.ascontiguousarray(out) return np.ascontiguousarray(out)
@ -329,7 +329,7 @@ def DSYR_blas(A, x, alpha=1.):
:param alpha: scalar :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) symmetrify(A, upper=True)
def DSYR_numpy(A, x, alpha=1.): def DSYR_numpy(A, x, alpha=1.):