#include <iostream.h>
#include <iomanip.h>
#include <time.h>
#include <stdlib.h>

#define BLOCCO_A_FINE

#define DIM 610

#define TITOLO(x1,x2,x3) "Mat. "<<x1<<"x"<<x2<<" per mat. "<<x2<<"x"<<x3<<"\n"

int operand1 [DIM][DIM];
int operand2 [DIM][DIM];
int result   [DIM][DIM];

void randomFill(int m[DIM][DIM]);

void mulSqrMatr(int op1[DIM][DIM], int op2[DIM][DIM], int res[DIM][DIM]);

void sqrMatrPrint(int m[DIM][DIM]);

void main() {
   // generazione del seme per rand()
   unsigned int seed = time(NULL);
   long t1, t2;
   double durationMs;

   srand(seed);

   // inserimento valori in operand1
   randomFill(operand1);
   //sqrMatrPrint(operand1);
   cout << endl;
   // inserimento valori in operand2
   randomFill(operand2);
   //sqrMatrPrint(operand2);
   cout << endl;
   // calcolo di result = operand1 * operand2
   t1= clock();
   mulSqrMatr(operand1, operand2, result);
   t2= clock();
   durationMs=1000.0*(t2-t1)/CLOCKS_PER_SEC;
   cout << "t1=" << t1 << endl;
   cout << "t2=" << t2 << endl;
   cout << "differenza = " << durationMs << endl;
   cout << "coeff: " << 1000*durationMs/(DIM*DIM*DIM) << endl;
   // visualizzazione
   //sqrMatrPrint(result);
   cout << "fatto!" <<endl;

#ifdef BLOCCO_A_FINE
char c;
cin >> c;
#endif
}

void randomFill(int m[DIM][DIM]) {
   int i,j;
   for(i=0; i<DIM; i++)
      for(j=0; j<DIM; m[i][j++]=(rand()%200)-100);
}

void mulSqrMatr(int op1[DIM][DIM], int op2[DIM][DIM], int res[DIM][DIM]){
   int i, j, k;
   for(i=0; i<DIM; i++) {
       for(j=0; j<DIM; j++) {
           res[i][j]=0;
           for(k=0; k<DIM; k++) {
               res[i][j]+=op1[i][k]*op2[k][j];
           }
       }
   }
}

void sqrMatrPrint(int m[DIM][DIM]) {
   int i,j;
   for(i=0; i<DIM; i++) {
       for(j=0; j<DIM; j++)
           cout << setw(5) << m[i][j] << " ";
       cout <<endl;
   }
}