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

#define SIZE 2048

double a[SIZE][SIZE], b[SIZE][SIZE], c[SIZE][SIZE];

int main(int argc, char* argv[])
{
	int i, j, k;
	int size;
	if (argc == 1) {
		size = 100;
	} else {
		size = atoi(argv[1]);
	}

	if ((size <1) || (size > SIZE)) {
		printf("Size (%d) must be between 1 and %d\n", size, SIZE);
		exit(0);
	}
	
	srand((unsigned)time(NULL));

	for (i=0; i<size; i++) {
		for (j=0; j<size; j++) {
			a[i][j] = 0.0;
			b[i][j] = drand48();
			c[i][j] = drand48();
		}
	}
	
	for (i=0; i<size; i++) {
		for (j=0; j<size; j++) {
			for (k=0; k<size; k++) { 
				a[i][j] = a[i][j] + b[i][k] * c[k][j];
			}
		}
	}
}