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