Code for computing the complete (2-by-2) CS decomposition of a partitioned unitary matrix.
My CSD code is now included in LAPACK! To build from source, download the latest LAPACK release from Netlib. (See below.) The driver routines are SORCSD, DORCSD, CUNCSD, and ZUNCSD. The major support routines are SORBDB, DORBDB, CUNBDB, ZUNBDB, SBBCSD, DBBCSD, CBBCSD, and ZBBCSD.
Native Matlab code is also available. It is numerically stable but slow, and it places some artificial restrictions on the input. If you would like a copy, please send me an email (firstname.lastname@example.org).
This material is based upon work supported by the National Science Foundation under Grant No. 0914559. Any opinions, findings and conclusions or recomendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation (NSF).home
Copyright 2012, Brian Sutton.