1*bf2c3715SXin Li #include <iostream> 2*bf2c3715SXin Li #include <Eigen/Core> 3*bf2c3715SXin Li using namespace Eigen; 4*bf2c3715SXin Li using namespace std; 5*bf2c3715SXin Li main()6*bf2c3715SXin Liint main() 7*bf2c3715SXin Li { 8*bf2c3715SXin Li typedef Matrix<SCALAR,Eigen::Dynamic,1> Vec; 9*bf2c3715SXin Li Vec v(SIZE); 10*bf2c3715SXin Li v.setZero(); 11*bf2c3715SXin Li v[0] = 1; 12*bf2c3715SXin Li v[1] = 2; 13*bf2c3715SXin Li for(int i = 0; i < 1000000; i++) 14*bf2c3715SXin Li { 15*bf2c3715SXin Li v.coeffRef(0) += v.sum() * SCALAR(1e-20); 16*bf2c3715SXin Li } 17*bf2c3715SXin Li cout << v.sum() << endl; 18*bf2c3715SXin Li } 19