#include <stdio.h>
#include <stdlib.h>
#include <math.h>


void mat_vec(double* a, double* x, double* y, int n)
{

  int i,j;

  for(j=0;j<n;j++){
    y[i]=0.0;
    for(i=0;i<n;i+=8){
      y[i]+=a[i+n*j]*x[j];
      y[i+1]+=a[(i+1)+n*j]*x[j];
      y[i+2]+=a[(i+2)+n*j]*x[j];
      y[i+3]+=a[(i+3)+n*j]*x[j];
      y[i+4]+=a[(i+4)+n*j]*x[j];
      y[i+5]+=a[(i+5)+n*j]*x[j];
      y[i+6]+=a[(i+6)+n*j]*x[j];
      y[i+7]+=a[(i+7)+n*j]*x[j];
    }
  }

}