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