#include <iostream>
#include <cstdlib>  // for random number generator
#include <ctime>
#include <fstream> // for writing to the file
using namespace std;

int main()
{

    int array_sizes[6] = {32,64,128,256,512,1024};

    int m;
    ofstream results; //output the results
    results.open ("results.txt");

    for(int i = 0; i < 6; i++){
        for (int k = 0; k < 3; k ++){

            clock_t start = clock();
            m = array_sizes[i];

            int **A = new int*[m];
            int **B = new int*[m];
            int **product = new int*[m];

            for ( int i = 0 ; i < m ; i++)
            {
                A[i] = new int[m];
                B[i] = new int[m];
                product[i] = new int[m];
            }

            // Initialize A[I][K]
            // I x K == row * inner
            for(int row = 0; row < m; row++) {
                for(int inner = 0; inner < m; inner++){
                    A[row][inner] = 5;
                }
            }


            for(int row = 0; row < m; row++) {
                for(int inner = 0; inner < m; inner++){
                    B[row][inner] = 3;
                    //cout << B[row][inner] << " ";
                }
                //cout << "\n";
            }

            // main routine for multiplication
            for(int row = 0; row < m; row++) {
                for(int col = 0; col < m; col++) {
                    product[row][col] = 0;
                    for(int inner = 0; inner < m; inner++) {
                        product[row][col] += A[row][inner] * B[inner][col];
                    }
                    //cout << product[row][col] << " ";
                }
                //cout << "\n";
            }

            //Deallocate
            for ( int i = 0 ; i < m ; i++)
            {
                delete[] A[i];
                delete[] B[i];
                delete[] product[i];
            }
            delete[] A;
            delete[] B;
            delete[] product;


            clock_t end = clock();
            double cpu_time = static_cast<double>(end - start)/CLOCKS_PER_SEC;

            results << cpu_time << endl;
        }
        results << m << endl;
    }

    results.close();
    return 0;
}