#include "header.h" int main (void){ int m, n, p, runs; int i, j, k; m=n=p=1000; runs = 1; double my_t, dgemm_t; dgemm_t = gsl_dgemm_bench(m,n,p,runs); my_t = myMult_bench(m,n,p,runs); printf("time %e (Seconds) in GSL\n", dgemm_t); printf("time %e (Seconds) in mine\n", my_t); printf("GSL speed up %e\n", my_t/dgemm_t); /* for(i=100;i<=m;i+=50){ for(j=100;j<=n;j+=50){ for(k=100;k<=p;k+=50){ dgemm_t = dgemm_bench(i,j,k,runs); my_t = myMult_bench(i,j,k,runs); printf("%d\t%d\t%d\t%e\t%e\t%e\n",i,j,k,my_t,dgemm_t,my_t/dgemm_t); } } printf("\n"); } */ return 0; }