vp.linalg.schur

The real Schur decomposition A = Z T Z^T: T quasi-triangular, Z orthogonal. As scipy.linalg.schur(A, output="real"). T is not unique, so compare through the residual and the eigenvalues, not entry by entry.

Part of the distribution vipyr; it runs from a checkout of the repository. Planned on PyPI as pip install vipyr.