1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
| #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;
| }
|
|
