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];
       }
     }
     """