main()
{
double * a, *b, *c, *savedC;
int s, i;
long t;
struct timeval start, end;
for(s = 500; s <= 5000; s+= 500) {
a = (double*)malloc(s * s * sizeof(double));
b = (double*)malloc(s * s * sizeof(double));
c = (double*)malloc(s * s * sizeof(double));
seedMatrixRandomly(a, s);
seedMatrixRandomly(b, s);
seedMatrixRandomly(c, s);
printf("Starting multiplication [size %d]\n", s);
gettimeofday(&start, NULL);
// Doing the multiplication //////
char n = 'n';
double scalar = 1.0;
dgemm_(&n, &n, &s, &s, &s, &scalar, a, &s, b, &s, &scalar, c, &s);
//////////////////////////////////
gettimeofday(&end, NULL);
t = ((end.tv_sec - start.tv_sec) * 1000 + (end.tv_usec - start.tv_usec)/1000.0) + 0.5;
printf("Finished multiplication [time %d], calculating difference\n", (int)t);
if(1) {
savedC = c;
c = (double*)malloc(s * s * sizeof(double));
seedMatrixRandomly(c, s);
doMatrixMultiplication(a, b, c, s);
double totalOff = 0;
for(i = 0 ; i < s * s; i++) {
totalOff += c[i] - savedC[i];
}
free(savedC);
printf("Difference calculated [value %.13f]\n\n", totalOff);
}
free(a);
free(b);
free(c);
}
return 0;
}
int main(void){
doMatrixMultiplication(3);
return 0;
}
