Towards a more robust algorithm for computing the restricted singular value decomposition

arXiv2002.04828MaRDI QIDQ6334643FDOQ6334643

Authors: Ian Zwaan 

Publication date: 12 February 2020

Abstract: A new algorithm to compute the restricted singular value decomposition of dense matrices is presented. Like Zha's method cite{Zha92}, the new algorithm uses an implicit Kogbetliantz iteration, but with four major innovations. The first innovation is a useful quasi-upper triangular generalized Schur form that just requires orthonormal transformations to compute. Depending on the application, this Schur form can be used instead of the full decomposition. The second innovation is a new preprocessing phase that requires fewer rank determinations than previous methods. The third innovation is a numerically stable RSVD algorithm for $2imes2$ upper-triangular matrices, which forms a key component of the implicit Kogbetliantz iteration. The fourth innovation is an alternative scaling for the restricted singular triplets that results in elegant formulas for their computation. Beyond these four innovations, the qualitative (numerical) characteristics of the algorithm are discussed extensively. Some numerical challenges in the (optional) postprocessing phase are considered too; though, their solutions require further research. Numerical tests and examples confirm the effectiveness of the method.

Mathematics Subject Classification ID

Computational methods for sparse matrices (65F50) Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Ill-posedness and regularization problems in numerical linear algebra (65F22) Numerical methods for inverse problems for integral equations (65R32) Numerical methods for ill-posed problems for integral equations (65R30)

This page was built for publication: Towards a more robust algorithm for computing the restricted singular value decomposition

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6334643)