#include <stdlib.h>
#include <iostream>
#include <vector>
#include "nbstopwatch.hh"

/*
 * Funktion zur Matrixmultiplikation 
 * C=A*B 
 * fuer Matrizen der Groesse NxN 
 */
void matmult(int N,
    std::vector<double> & A,
    std::vector<double> & B,
    std::vector<double> & C)
{
    // TODO
}

/*
 * Hauptprogramm
 */
int main()
{
    // Variiere die Groesse der Matrix
    for (int N = 2; N <= 512; N*=2)
    {
        // mehrfach multiplizieren
        int mehr = 512*512*512 / N/N/N;
        
        std::vector<double> A(N*N);
        std::vector<double> B(N*N);
        std::vector<double> C(N*N);
        
        // Zufallszahlengenerator initialisieren
        srand(666);
        
        // Matrizen initialisieren
        for (int i=0; i<N; i++)
            for (int j=0; j<N; j++)
            {
                A[i*N+j] = 1.0 / (1.0 + rand() );
                B[i*N+j] = 1.0 / (1.0 + rand() );
                C[i*N+j] = 0.0;
            }

        double start = get_time();
        // Matrix Multiplikation
        for (int i=0; i<mehr; i++)
        {
            matmult(N, A, B, C);
        }
        double dauer = get_time() - start;
        // Ergebnis im gnuplot format ausgeben
        std::cout << N << "    " << 2.0*N*N*N*mehr/dauer << std::endl;
    }
    
    return 0;
}