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
11 using namespace std;
12 using namespace Eigen;
13
14 #ifndef REPEAT
15 #define REPEAT 40000000
16 #endif
17
18 #ifndef SCALAR
19 #define SCALAR double
20 #endif
21
main(int argc,char * argv[])22 int 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