| Benoit Jacob | e093b43 | 2009-07-06 17:12:10 +0200 | [diff] [blame] | 1 | typedef Matrix<float,3,3> Matrix3x3; |
| 2 | Matrix3x3 m = Matrix3x3::Random(); |
| 3 | Matrix3f y = Matrix3f::Random(); |
| 4 | cout << "Here is the matrix m:" << endl << m << endl; |
| 5 | cout << "Here is the matrix y:" << endl << y << endl; |
| 6 | Matrix3f x; |
| Gael Guennebaud | 7579360 | 2009-07-27 18:50:39 +0200 | [diff] [blame] | 7 | m.householderQr().solve(y, &x); |
| Benoit Jacob | e093b43 | 2009-07-06 17:12:10 +0200 | [diff] [blame] | 8 | assert(y.isApprox(m*x)); |
| 9 | cout << "Here is a solution x to the equation mx=y:" << endl << x << endl; |