1// g++ -O3 -DNDEBUG -DMATSIZE=<x> benchmark.cpp -o benchmark && time ./benchmark
2
3#include <iostream>
4
5#include <Eigen/Core>
6
7#ifndef MATSIZE
8#define MATSIZE 3
9#endif
10
11using namespace std;
12using namespace Eigen;
13
14#ifndef REPEAT
15#define REPEAT 40000000
16#endif
17
18#ifndef SCALAR
19#define SCALAR double
20#endif
21
22int main(int argc, char *argv[])
23{
24    Matrix<SCALAR,MATSIZE,MATSIZE> I = Matrix<SCALAR,MATSIZE,MATSIZE>::Ones();
25    Matrix<SCALAR,MATSIZE,MATSIZE> m;
26    for(int i = 0; i < MATSIZE; i++)
27        for(int j = 0; j < MATSIZE; j++)
28        {
29            m(i,j) = (i+MATSIZE*j);
30        }
31    asm("#begin");
32    for(int a = 0; a < REPEAT; a++)
33    {
34        m = Matrix<SCALAR,MATSIZE,MATSIZE>::Ones() + 0.00005 * (m + (m*m));
35    }
36    asm("#end");
37    cout << m << endl;
38    return 0;
39}
40