R equivalent to Matlab reordering of Schur factorization?

Question

Is there an equivalent MATLAB ordschur function (documentation here ) in R?

The function reorders the Schur factorization X = U * T * U 'generated by the schur function and returns the reordered Schur matrix TS and the cumulative orthogonal transformation US such that X = US * TS * US'. I am particularly interested in the "lhp" method, also described in the MATLAB documentation documentation.

Note that in the packet matrix there is a Schur function (see CRAN documentation here ) R, which performs Schur decomposition and eigenvalues ​​of the square matrix. Update: this function also returns the unitary orthogonal matrix U.

Answer 1

As far as I know, MATLAB uses the ?TRSEN function from LAPACK to perform reordering. You can look at the limited implementation here . To bring this functionality to R, you can implement this procedure yourself.