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