bench/bench_sum.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

#include <iostream>
#include <Eigen/Core>
using namespace Eigen;
using namespace std;

int main() 
{
  typedef Matrix<SCALAR,Eigen::Dynamic,1> Vec;
  Vec v(SIZE);
  v.setZero();
  v[0] = 1;
  v[1] = 2;
  for(int i = 0; i < 1000000; i++)
  {
    v.coeffRef(0) += v.sum() * SCALAR(1e-20);
  }
  cout << v.sum() << endl;
}