#include "matrixmul.h"

void matrixmul(
      mat_a_t a[MAT_A_ROWS][MAT_A_COLS],
      mat_b_t b[MAT_B_ROWS][MAT_B_COLS],
      result_t res[MAT_A_ROWS][MAT_B_COLS])
{

result_t temp;

// Generate the expected result
// Iterate over the rows of the A matrix
I_LOOP:  for(int i = 0; i < MAT_A_ROWS; i++) {
J_LOOP:     for(int j = 0; j < MAT_B_COLS; j++) {
     // Iterate over the columns of the B matrix
     // Do the inner product of a row of A and col of B
   K_LOOP:      for(int k = 0; k < MAT_B_ROWS; k++) {
#pragma AP PIPELINE
       if(k==0)
    temp = 0;
       temp += a[i][k] * b[k][j];
       if(k == (MAT_B_ROWS - 1))
    res[i][j] = temp;
     }
   }
}
}