trying optimize 3mm

polybench-c-4.2.1-beta/linear-algebra/kernels/3mm/3mm_optimized.c

@@ -91,19 +91,23 @@ void kernel_3mm(int ni, int nj, int nk, int nl, int nm,
       for (j = 0; j < _PB_NL; j++)
         F[i][j] += C[i][k] * D[k][j];
 
-  /* E := A*B */  
-  /* G := E*F */
-  for (i = 0; i < _PB_NI; i++){
+  /* E := A*B */
+  for (i = 0; i < _PB_NI; i++)
     for (j = 0; j < _PB_NJ; j++)
     	E[i][j] = SCALAR_VAL(0.0);
+
+  for (i = 0; i < _PB_NI; i++)
+    for (k = 0; k < _PB_NK; ++k)
+      for (j = 0; j < _PB_NJ; j++)
+        E[i][j] += A[i][k] * B[k][j];
+  
+  /* G := E*F */
+  for (i = 0; i < _PB_NI; i++){
     for (j = 0; j < _PB_NL; j++)
 	    G[i][j] = SCALAR_VAL(0.0);
   }
 
   for (i = 0; i < _PB_NI; i++){
-    for (k = 0; k < _PB_NK; ++k)
-      for (j = 0; j < _PB_NJ; j++)
-        E[i][j] += A[i][k] * B[k][j];
     for (k = 0; k < _PB_NJ; ++k)
       for (j = 0; j < _PB_NL; j++)
         G[i][j] += E[i][k] * F[k][j];