1*bf2c3715SXin Li #include <iostream> 2*bf2c3715SXin Li #include <Eigen/Dense> 3*bf2c3715SXin Li 4*bf2c3715SXin Li using namespace std; 5*bf2c3715SXin Li using namespace Eigen; 6*bf2c3715SXin Li main()7*bf2c3715SXin Liint main() 8*bf2c3715SXin Li { 9*bf2c3715SXin Li Matrix2d A; 10*bf2c3715SXin Li A << 2, 1, 11*bf2c3715SXin Li 2, 0.9999999999; 12*bf2c3715SXin Li FullPivLU<Matrix2d> lu(A); 13*bf2c3715SXin Li cout << "By default, the rank of A is found to be " << lu.rank() << endl; 14*bf2c3715SXin Li lu.setThreshold(1e-5); 15*bf2c3715SXin Li cout << "With threshold 1e-5, the rank of A is found to be " << lu.rank() << endl; 16*bf2c3715SXin Li } 17