Brian Silverman | 72890c2 | 2015-09-19 14:37:37 -0400 | [diff] [blame] | 1 | MatrixXf A = MatrixXf::Random(4,4); |
2 | RealSchur<MatrixXf> schur(4); | ||||
3 | schur.compute(A, /* computeU = */ false); | ||||
4 | cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; | ||||
5 | schur.compute(A.inverse(), /* computeU = */ false); | ||||
6 | cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl; |