2010-04-13 01:14:32 +08:00
|
|
|
MatrixXf A = MatrixXf::Random(4,4);
|
|
|
|
|
RealSchur<MatrixXf> schur(4);
|
2010-06-01 01:17:47 +08:00
|
|
|
schur.compute(A, /* computeU = */ false);
|
2010-04-13 01:14:32 +08:00
|
|
|
cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl;
|
2010-06-01 01:17:47 +08:00
|
|
|
schur.compute(A.inverse(), /* computeU = */ false);
|
2010-04-13 01:14:32 +08:00
|
|
|
cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl;
|
