| Matrix3f A(3,3); | |
| A << 1,2,3, 4,5,6, 7,8,10; | |
| PartialPivLU<Matrix3f> luOfA(A); // compute LU decomposition of A | |
| Vector3f b; | |
| b << 3,3,4; | |
| Vector3f x; | |
| x = luOfA.solve(b); | |
| cout << "The solution with right-hand side (3,3,4) is:" << endl; | |
| cout << x << endl; | |
| b << 1,1,1; | |
| x = luOfA.solve(b); | |
| cout << "The solution with right-hand side (1,1,1) is:" << endl; | |
| cout << x << endl; |