#include <stdio.h>
#include <mkl.h>
#include <mkl_cblas.h>
void matrix_generator(double *,double *, double*, int);

int main()
{
int i, j,n,n2;
n=15000;
n2=200;
double* m = (double*)malloc( n*n*sizeof(double));
double* x = (double*)malloc( n*sizeof(double));
double* y = (double*)malloc( n*sizeof(double));
mkl_set_dynamic(0);
mkl_set_num_threads(4);

matrix_generator(m,x,y,n);


for(i=0;i<n2;i++)
  cblas_dgemv(CblasRowMajor, CblasNoTrans, n, n, 1.0, m, n,x, 1, 0.0, y, 1);

  return 0;
}

void matrix_generator(double *m,double *x, double *y,int n)
{
int i,j;
for(i=0;i<n;i++)
{
x[i] = (i+1)*(i+1)*1.00;
y[i] = 0;
for(j=0;j<n;j++)
m[i*n+j] = (i+j+2)*(j+1)*1.00;
}
}