source: trunk/platforms/caba-ring-ccxcachev4_memcachev4-mips32el/soft/matrix_multiplication/matrix_multiplication.c @ 134

#include "matrix_multiplication.h"
#include "stdio.h"
#include "system.h"

void matrix_multiplication_st (int ** a,
                               int ** b,
                               int ** d,
                               unsigned int n)
{
  int x, y;

  for (x=0; x<n; ++x)
    {
      printf("%d, ",x);

      for (y=0; y<n; ++y)
        {
          int i;

          int tmp = 0;
          for (i=0; i<n; ++i)
            tmp += a[x][i] * b[i][y];
          d[x][y] = tmp;
        }
    }
  
  printf("finish\n");
}

static int lock;
static int next;

void matrix_multiplication_mt (int ** a,
                               int ** b,
                               int ** d,
                               unsigned int n,
                               int lock_by_line)
{
  int x, y;

  while (1)
    {
      lock_lock(&lock);
      
      x=next++;
      
      lock_unlock(&lock);

      if (x >= n)
        break;
      
      printf("%d, ",x);
      for (y=0; y<n; ++y)
        {
          int line  =(lock_by_line)?x:y;
          int column=(lock_by_line)?y:x;
          
          int i;
          int tmp = 0;
          for (i=0; i<n; ++i)
            tmp += a[column][i] * b[i][line];
          d[column][line] = tmp;
        }
    }
  
  printf("finish\n");
}