ground/gcs/src/libs/eigen/bench/benchmark.cpp

   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
  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 }