Computes the Schur decomposition of a square matrix.
(u, s) = SCHUR(a)
A square matrix.
(u, s) = schur(matrix) returns the separate unit Schur and Schur components such that
For matrix A, SCHUR and USCHUR produce matrixes such that:
If the matrix is real, SCHUR returns the real Schur form which has the real eigenvalues on the diagonal and the complex eigenvalues in 2-by-2 blocks on the diagonal.
If the matrix is complex, SCHUR returns the complex Schur form which is upper triangular with the eigenvalues of the matrix on the diagonal.
See DADiSP/MatrixXL to significantly optimize SCHUR.