An image algorithm for the bidiagonal SVD Original Research Article

The relatively robust representations (RRR) algorithm computes the eigendecomposition of a symmetric tridiagonal matrix T with an image complexity. This article discusses how this method can be extended to the bidiagonal SVD B=UΣVT. It turns out that using the RRR algorithm as a black box to compute BTB=VΣ2VT and BBT=UΣ2UT separately may give poor results for short parallelBV−UΣshort parallel. The use of the standard Jordan–Wielandt representation can fail as well if clusters of tiny singular values are present. A solution is to work on BTB and to keep factorizations of BBT implicitly. We introduce a set of coupling transformations which allow us to replace the representation image by a more stable representation image , where image is a diagonal matrix. Numerical results of our implementation are compared with the LAPACK routines DSTEGR, DBDSQR and DBDSDC.