00001 MatrixXcf A = MatrixXcf::Random(4,4); 00002 ComplexSchur<MatrixXcf> schur(4); 00003 schur.compute(A); 00004 cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; 00005 schur.compute(A.inverse()); 00006 cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl;