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
;
