From 542759353b63a5ca3a04fe9f0c96c2b7972c1f55 Mon Sep 17 00:00:00 2001
From: James Hensman <james.hensman@gmail.com>
Date: Mon, 13 May 2013 17:23:49 +0100
Subject: [PATCH] minor changes to the symmetrify function

---
 GPy/kern/linear.py | 2 +-
 GPy/util/linalg.py | 8 ++++++--
 2 files changed, 7 insertions(+), 3 deletions(-)

diff --git a/GPy/kern/linear.py b/GPy/kern/linear.py
index 16ef2499..af3e60ea 100644
--- a/GPy/kern/linear.py
+++ b/GPy/kern/linear.py
@@ -203,7 +203,7 @@ class linear(kernpart):
                 target_mu(n,q) += factor*tmp;
                 target_S(n,q) += factor*AZZA_2(q,m,mm,q);
               }
-            }  
+            }
           }
         }
         """
diff --git a/GPy/util/linalg.py b/GPy/util/linalg.py
index 64155dca..555345d9 100644
--- a/GPy/util/linalg.py
+++ b/GPy/util/linalg.py
@@ -257,16 +257,20 @@ def symmetrify(A,upper=False):
     N,M = A.shape
     assert N==M
     c_contig_code = """
+    int iN;
     for (int i=1; i<N; i++){
+      iN = i*N;
       for (int j=0; j<i; j++){
-        A[i+j*N] = A[i*N+j];
+        A[i+j*N] = A[iN+j];
       }
     }
     """
     f_contig_code = """
+    int iN;
     for (int i=1; i<N; i++){
+      iN = i*N;
       for (int j=0; j<i; j++){
-        A[i*N+j] = A[i+j*N];
+        A[iN+j] = A[i+j*N];
       }
     }
     """