mt/al_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, x;
   9     data_t temp0, temp1, temp2, temp3, temp4, temp5, temp6, temp7;
  10     data_t temp8, temp9, temp10, temp11, temp12, temp13, temp14, temp15;
  11
  12     //complete Q1
  13    if(coreid > 1) return;
  14    if(coreid == 0)   {
  15       for(j = 0; j < 32; j++) {
  16          temp0  = C[j*lda];
  17          temp1  = C[1  + j*lda];
  18          temp2  = C[2  + j*lda];
  19          temp3  = C[3  + j*lda];
  20          temp4  = C[4  + j*lda];
  21          temp5  = C[5  + j*lda];
  22          temp6  = C[6  + j*lda];
  23          temp7  = C[7  + j*lda];
  24          temp8  = C[8  + j*lda];
  25          temp9  = C[9  + j*lda];
  26          temp10 = C[10 + j*lda];
  27          temp11 = C[11 + j*lda];
  28          temp12 = C[12 + j*lda];
  29          temp13 = C[13 + j*lda];
  30          temp14 = C[14 + j*lda];
  31          temp15 = C[15 + j*lda];
  32          for(k = 0; k < 32; k++) {
  33             temp0  += A[j*lda + k] * B[k*lda];
  34             temp1  += A[j*lda + k] * B[1+k*lda];
  35             temp2  += A[j*lda + k] * B[2+k*lda];
  36             temp3  += A[j*lda + k] * B[3+k*lda];
  37             temp4  += A[j*lda + k] * B[4+k*lda];
  38             temp5  += A[j*lda + k] * B[5+k*lda];
  39             temp6  += A[j*lda + k] * B[6+k*lda];
  40             temp7  += A[j*lda + k] * B[7+k*lda];
  41             temp8  += A[j*lda + k] * B[8+k*lda];
  42             temp9  += A[j*lda + k] * B[9+k*lda];
  43             temp10 += A[j*lda + k] * B[10+k*lda];
  44             temp11 += A[j*lda + k] * B[11+k*lda];
  45             temp12 += A[j*lda + k] * B[12+k*lda];
  46             temp13 += A[j*lda + k] * B[13+k*lda];
  47             temp14 += A[j*lda + k] * B[14+k*lda];
  48             temp15 += A[j*lda + k] * B[15+k*lda];
  49          }
  50          C[j*lda] = temp0;
  51          C[1  + j*lda] = temp1;
  52          C[2  + j*lda] = temp2;
  53          C[3  + j*lda] = temp3;
  54          C[4  + j*lda] = temp4;
  55          C[5  + j*lda] = temp5;
  56          C[6  + j*lda] = temp6;
  57          C[7  + j*lda] = temp7;
  58          C[8  + j*lda] = temp8;
  59          C[9  + j*lda] = temp9;
  60          C[10 + j*lda] = temp10;
  61          C[11 + j*lda] = temp11;
  62          C[12 + j*lda] = temp12;
  63          C[13 + j*lda] = temp13;
  64          C[14 + j*lda] = temp14;
  65          C[15 + j*lda] = temp15;
  66       }
  67    }
  68
  69    if( coreid == 1 || ncores == 1) {
  70       for(j = 0; j < 32; j++) {
  71          temp0  = C[16 + j*lda];
  72          temp1  = C[17 + j*lda];
  73          temp2  = C[18 + j*lda];
  74          temp3  = C[19 + j*lda];
  75          temp4  = C[20 + j*lda];
  76          temp5  = C[21 + j*lda];
  77          temp6  = C[22 + j*lda];
  78          temp7  = C[23 + j*lda];
  79          temp8  = C[24 + j*lda];
  80          temp9  = C[25 + j*lda];
  81          temp10 = C[26 + j*lda];
  82          temp11 = C[27 + j*lda];
  83          temp12 = C[28 + j*lda];
  84          temp13 = C[29 + j*lda];
  85          temp14 = C[30 + j*lda];
  86          temp15 = C[31 + j*lda];
  87          for(k = 0; k < 32; k++) {
  88             temp0  += A[j*lda + k] * B[16 + k*lda];
  89             temp1  += A[j*lda + k] * B[17 + k*lda];
  90             temp2  += A[j*lda + k] * B[18 + k*lda];
  91             temp3  += A[j*lda + k] * B[19 + k*lda];
  92             temp4  += A[j*lda + k] * B[20 + k*lda];
  93             temp5  += A[j*lda + k] * B[21 + k*lda];
  94             temp6  += A[j*lda + k] * B[22 + k*lda];
  95             temp7  += A[j*lda + k] * B[23 + k*lda];
  96             temp8  += A[j*lda + k] * B[24 + k*lda];
  97             temp9  += A[j*lda + k] * B[25 + k*lda];
  98             temp10 += A[j*lda + k] * B[26 + k*lda];
  99             temp11 += A[j*lda + k] * B[27 + k*lda];
 100             temp12 += A[j*lda + k] * B[28 + k*lda];
 101             temp13 += A[j*lda + k] * B[29 + k*lda];
 102             temp14 += A[j*lda + k] * B[30 + k*lda];
 103             temp15 += A[j*lda + k] * B[31 + k*lda];
 104          }
 105          C[16 + j*lda] = temp0;
 106          C[17 + j*lda] = temp1;
 107          C[18 + j*lda] = temp2;
 108          C[19 + j*lda] = temp3;
 109          C[20 + j*lda] = temp4;
 110          C[21 + j*lda] = temp5;
 111          C[22 + j*lda] = temp6;
 112          C[23 + j*lda] = temp7;
 113          C[24 + j*lda] = temp8;
 114          C[25 + j*lda] = temp9;
 115          C[26 + j*lda] = temp10;
 116          C[27 + j*lda] = temp11;
 117          C[28 + j*lda] = temp12;
 118          C[29 + j*lda] = temp13;
 119          C[30 + j*lda] = temp14;
 120          C[31 + j*lda] = temp15;
 121       }
 122    }
 123 }