tests/c_files/matmul.c

   1 /* matmul.c  128*128 matrix multiply */
   2 #include <stdio.h>
   3 #define N 128
   4
   5 int main()
   6 {
   7   int   i,j,k;
   8   float a[N][N], b[N][N], c[N][N];
   9
  10   /* We read matrix a */
  11   for(i=0; i<N; i++)
  12     for(j=0; j<N; j++)
  13       scanf(" %f",&a[i][j]);
  14
  15   /* We read matrix b */
  16   for(i=0; i<N; i++)
  17     for(j=0; j<N; j++)
  18       scanf(" %f",&b[i][j]);
  19
  20   /* c = a * b */
  21 #pragma scop
  22   for(i=0; i<N; i++)
  23     for(j=0; j<N; j++)
  24     {
  25       c[i][j] = 0.0;
  26       for(k=0; k<N; k++)
  27         c[i][j] = c[i][j] + a[i][k]*b[k][j];
  28     }
  29 #pragma endscop
  30
  31   /* We print matrix c */
  32   for(i=0; i<N; i++)
  33   {
  34     for(j=0; j<N; j++)
  35       printf("%6.2f ",c[i][j]);
  36     printf("\n");
  37   }
  38
  39   return 0;
  40 }