mt/cs_matmul.c

   1 #include "stdlib.h"
   2
   3 #include "util.h"
   4
   5 #include "dataset.h"
   6 void __attribute__((noinline)) matmul(const int coreid, const int ncores, const int lda,  const data_t A[], const data_t B[], data_t C[] )
   7 {
   8    int i, j, k, n, m, c1, c2;
   9    for ( j = coreid; j < lda; j += 2*ncores ) {
  10       for ( i = 0; i < lda; i += 1 ){
  11          c1 = 0;     //global vars c1, c2
  12          c2 = 0;
  13          for ( k = 0; k < lda; k += 1 ) {
  14             c1 += A[j * lda + k] * B[k*lda + i];
  15             c2 += A[(j+ncores) * lda + k] * B[k*lda + i];
  16          }
  17
  18          C[i + j * lda] = c1;
  19          C[i + (j+ncores) * lda] = c2;
  20          barrier(ncores);
  21       }
  22    }
  23
  24 }