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

//gsl blas libs
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_cblas.h>

double 	 gsl_dgemm_bench( int m, int n, int p, int runs);
double 	 mkl_dgemm_bench( int m, int n, int p, int runs);
double       myMult_bench(    int m, int n, int p, int runs);
gsl_matrix * Create_Rand_Matrix(int m, int n);
void 		 my_dgemm(gsl_matrix * A, gsl_matrix * B, gsl_matrix * C);
void 		 gsl_dgemm(gsl_matrix * A, gsl_matrix * B, gsl_matrix * C);
void 		 mkl_dgemm(gsl_matrix * A, gsl_matrix * B, gsl_matrix * C);