| Matrix3f A(3,3); | |
| A << 1,2,3, 4,5,6, 7,8,10; | |
| Matrix<float,3,2> B; | |
| B << 3,1, 3,1, 4,1; | |
| Matrix<float,3,2> X; | |
| A.partialLu().solve(B, &X); | |
| cout << "The solution with right-hand side (3,3,4) is:" << endl; | |
| cout << X.col(0) << endl; | |
| cout << "The solution with right-hand side (1,1,1) is:" << endl; | |
| cout << X.col(1) << endl; |