| Jitse Niesen | 73d3a27 | 2010-04-12 18:14:32 +0100 | [diff] [blame^] | 1 | MatrixXf A = MatrixXf::Random(4,4); |
| 2 | RealSchur<MatrixXf> schur(4); | ||||
| 3 | schur.compute(A); | ||||
| 4 | cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; | ||||
| 5 | schur.compute(A.inverse()); | ||||
| 6 | cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl; | ||||