xref: /aosp_15_r20/external/eigen/doc/examples/TutorialLinAlgSetThreshold.cpp (revision bf2c37156dfe67e5dfebd6d394bad8b2ab5804d4)
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 Li int 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