matvec-serial.c.txt

#include <stdio.h>
#define N 10000

int matvec(int n, double a[n][n],
           double x[n], double y[n]) {
  int i, j, sum;

#pragma ivdep
  for (i = 0; i < n; i++) {
    sum = 0.0;
    for (j = 0; j < n; j++) {
      sum += a[i][j]*x[j];
    }
    y[i] = sum;
  }
  return 0;
}

int main() {
  double a[N][N];
  double b[N], c[N];
  int i, j;

#pragma ivdep
  for (i=0; i<N; i++) {
     b[i] = i;
     for (j=0; j<N; j++) {
        a[i][j] = 0.0;
     }
     a[i][i] = 1.0;
  }
  matvec(N,a,b,c);
//  for (i=0; i<N; i++) {
      printf("%f\n", c[N-1]);
//  }
  return 0;
}