ComplexSchur_compute.cpp 301 B

123456
  1. MatrixXcf A = MatrixXcf::Random(4,4);
  2. ComplexSchur<MatrixXcf> schur(4);
  3. schur.compute(A);
  4. cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl;
  5. schur.compute(A.inverse());
  6. cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl;