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

#define DIMENSION 1000


int main(int argc, char** argv) {
    double A[DIMENSION][DIMENSION], B[DIMENSION][DIMENSION], C[DIMENSION][DIMENSION];
    int i, k, j;
    clock_t start, end;

    for (i=0; i<DIMENSION; i++) 
	for (k=0; k<DIMENSION; k++) {
	    A[i][k] = ((double) rand()) / ((double) RAND_MAX);
	    B[i][k] = ((double) rand()) / ((double) RAND_MAX);
	    C[i][k] = 0;
	}

    start = clock();
    for (i=0; i<DIMENSION; i++) 
	for (k=0; k<DIMENSION; k++) 
	    for (j=0; j<DIMENSION; j++) 
		C[i][j] += A[i][k] * B[k][j];
    end = clock();

    printf("Matrix-matrix multiplication on %dx%d elements: %d clock ticks\n",DIMENSION,DIMENSION,end-start);

   

    return 0;
}