From 3853721355f792481890115d00c1b2c5c9a1cf6d Mon Sep 17 00:00:00 2001 From: Andreas Date: Tue, 28 May 2013 15:50:18 +0100 Subject: [PATCH] Fixed symmetrify for when C/F compiler doesn't work --- GPy/util/linalg.py | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/GPy/util/linalg.py b/GPy/util/linalg.py index abfb1900..da23a0a8 100644 --- a/GPy/util/linalg.py +++ b/GPy/util/linalg.py @@ -313,7 +313,10 @@ def symmetrify(A,upper=False): elif A.flags['F_CONTIGUOUS'] and not upper: weave.inline(f_contig_code,['A','N'], extra_compile_args=['-O3']) else: - tmp = np.tril(A) + if upper: + tmp = np.tril(A.T) + else: + tmp = np.tril(A) A[:] = 0.0 A += tmp A += np.tril(tmp,-1).T