From 8c7b5f5d3ba9e612cc92b2c33259bfeeeb4cf38c Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Sat, 23 May 2015 09:31:15 +0100
Subject: [PATCH] Cosmetic

---
 lib/math/Grid_math_arith_mul.h | 9 ++++-----
 1 file changed, 4 insertions(+), 5 deletions(-)

diff --git a/lib/math/Grid_math_arith_mul.h b/lib/math/Grid_math_arith_mul.h
index 04b05bb4..fcdd60ee 100644
--- a/lib/math/Grid_math_arith_mul.h
+++ b/lib/math/Grid_math_arith_mul.h
@@ -16,19 +16,18 @@ strong_inline void mult(iScalar<rtype> * __restrict__ ret,const iScalar<mtype> *
 template<class rrtype,class ltype,class rtype,int N>
 strong_inline void mult(iMatrix<rrtype,N> * __restrict__ ret,const iMatrix<ltype,N> * __restrict__ lhs,const iMatrix<rtype,N> * __restrict__ rhs){
   for(int c1=0;c1<N;c1++){
-    int c3=0;
     for(int c2=0;c2<N;c2++){
-      mult(&ret->_internal[c1][c2],&lhs->_internal[c1][c3],&rhs->_internal[c3][c2]);
+      mult(&ret->_internal[c1][c2],&lhs->_internal[c1][0],&rhs->_internal[0][c2]);
     }
   }
-  for(int c3=1;c3<N;c3++){
-    for(int c1=0;c1<N;c1++){
+  for(int c1=0;c1<N;c1++){
+    for(int c3=1;c3<N;c3++){
       for(int c2=0;c2<N;c2++){
 	mac(&ret->_internal[c1][c2],&lhs->_internal[c1][c3],&rhs->_internal[c3][c2]);
       }
     }
   }
-    return;
+  return;
 }
 
 template<class rrtype,class ltype,class rtype,int N>