A Backward Stable Algorithm for Computing the CS Decomposition via the Polar Decomposition

Evan S. Gawlik, Yuji Nakatsukasa, Brian D. Sutton

We introduce a backward stable algorithm for computing the CS decomposition of a partitioned $2n\times n$ matrix with orthonormal columns, or a rank-deficient partial isometry. Th…