Brian Silverman | 72890c2 | 2015-09-19 14:37:37 -0400 | [diff] [blame] | 1 | Matrix3f m = Matrix3f::Random(); |
2 | Matrix3f y = Matrix3f::Random(); | ||||
3 | cout << "Here is the matrix m:" << endl << m << endl; | ||||
4 | cout << "Here is the matrix y:" << endl << y << endl; | ||||
5 | Matrix3f x; | ||||
6 | x = m.colPivHouseholderQr().solve(y); | ||||
7 | assert(y.isApprox(m*x)); | ||||
8 | cout << "Here is a solution x to the equation mx=y:" << endl << x << endl; |