#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <mkl.h>



int main(int argc, char* argv[])
{
    int ret_code;
    MM_typecode matcode;
    FILE *f;
    int M, N, nz;
    int i, *I, *J,*ia,*ja;
    double *val,*acsr;
    clock_t start, stop;

    //Open the file 
    if ((f = fopen(argv[1], "r")) == NULL){
	    printf("Cant open");
            exit(1);

     }

     M=400;
     N=400;
     nz=2860;
    /* reseve memory for matrices */

    I = (int *) malloc(nz * sizeof(int));
    J = (int *) malloc(nz * sizeof(int));
    val = (double *) malloc(nz * sizeof(double));

    for (i=0; i<nz; i++)
    {
        fscanf(f, "%d %d %lf", &I[i],&J[i], &val[i]);
    }
    if (f !=stdin) 
        fclose(f);
    for (i=0; i<10; i++)
       printf("%d %d %lf\n", I[i], J[i], val[i]);
    
    /* Parameters for conversion */

    double* x = (double *)malloc(N*sizeof(double)*N);
    double* y = (double *)malloc(N*sizeof(double)*N);
    
    for(i=0;i<N;i++)
        x[i]=1;

    //Do the SpMV
    printf("Here");
    start = clock();
    mkl_dcoogemv("N", &M, val, I, J,&nz, x, y);
    stop = clock();
    printf("Here");

    //Print the results
    printf("\n---------------RESULT----------------\n");
    for(i=0;i<10;i++)
	printf("%lf\n",y[i]);
    printf("Dgemm_multiply(). Elapsed time = %g seconds\n",((double)(stop - start)) / CLOCKS_PER_SEC);

    return 0;
}