#include <iostream> | |
#include <Eigen/Dense> | |
using namespace std; | |
using namespace Eigen; | |
int main() | |
{ | |
Matrix2d A; | |
A << 2, 1, | |
2, 0.9999999999; | |
FullPivLU<Matrix2d> lu(A); | |
cout << "By default, the rank of A is found to be " << lu.rank() << endl; | |
lu.setThreshold(1e-5); | |
cout << "With threshold 1e-5, the rank of A is found to be " << lu.rank() << endl; | |
} |