00001 #include <iostream> 00002 #include <Eigen/Dense> 00003 00004 using namespace std; 00005 using namespace Eigen; 00006 00007 int main() 00008 { 00009 Matrix2d A; 00010 A << 2, 1, 00011 2, 0.9999999999; 00012 FullPivLU<Matrix2d> lu(A); 00013 cout << "By default, the rank of A is found to be " << lu.rank() << endl; 00014 lu.setThreshold(1e-5); 00015 cout << "With threshold 1e-5, the rank of A is found to be " << lu.rank() << endl; 00016 }