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