1*bf2c3715SXin Li // g++ -O3 -DNDEBUG benchmarkX.cpp -o benchmarkX && time ./benchmarkX
2*bf2c3715SXin Li
3*bf2c3715SXin Li #include <iostream>
4*bf2c3715SXin Li
5*bf2c3715SXin Li #include <Eigen/Core>
6*bf2c3715SXin Li
7*bf2c3715SXin Li using namespace std;
8*bf2c3715SXin Li using namespace Eigen;
9*bf2c3715SXin Li
10*bf2c3715SXin Li #ifndef REPEAT
11*bf2c3715SXin Li #define REPEAT 10000
12*bf2c3715SXin Li #endif
13*bf2c3715SXin Li
14*bf2c3715SXin Li #ifndef SCALAR
15*bf2c3715SXin Li #define SCALAR float
16*bf2c3715SXin Li #endif
17*bf2c3715SXin Li
main(int argc,char * argv[])18*bf2c3715SXin Li int main(int argc, char *argv[])
19*bf2c3715SXin Li {
20*bf2c3715SXin Li typedef Matrix<SCALAR, Eigen::Dynamic, Eigen::Dynamic> Mat;
21*bf2c3715SXin Li Mat m(100, 100);
22*bf2c3715SXin Li m.setRandom();
23*bf2c3715SXin Li
24*bf2c3715SXin Li for(int a = 0; a < REPEAT; a++)
25*bf2c3715SXin Li {
26*bf2c3715SXin Li int r, c, nr, nc;
27*bf2c3715SXin Li r = Eigen::internal::random<int>(0,10);
28*bf2c3715SXin Li c = Eigen::internal::random<int>(0,10);
29*bf2c3715SXin Li nr = Eigen::internal::random<int>(50,80);
30*bf2c3715SXin Li nc = Eigen::internal::random<int>(50,80);
31*bf2c3715SXin Li m.block(r,c,nr,nc) += Mat::Ones(nr,nc);
32*bf2c3715SXin Li m.block(r,c,nr,nc) *= SCALAR(10);
33*bf2c3715SXin Li m.block(r,c,nr,nc) -= Mat::constant(nr,nc,10);
34*bf2c3715SXin Li m.block(r,c,nr,nc) /= SCALAR(10);
35*bf2c3715SXin Li }
36*bf2c3715SXin Li cout << m[0] << endl;
37*bf2c3715SXin Li return 0;
38*bf2c3715SXin Li }
39