A Cross-Product Free Jacobi–Davidson Type Method for Computing a Partial Generalized Singular Value Decomposition of a Large Matrix Pair

Huang Jinzhi,Jia Zhongxiao
DOI: https://doi.org/10.1007/s10915-022-02053-w
2022-01-01
Journal of Scientific Computing
Abstract:A Cross-Product Free (CPF) Jacobi-Davidson (JD) type method is proposed tocompute a partial generalized singular value decomposition (GSVD) of a largeregular matrix pair (A,B). It implicitly solves the mathematically equivalentgeneralized eigenvalue problem of (A^TA,B^TB) but does not explicitly formthe cross-product matrices and thus avoids the possible accuracy loss of thecomputed generalized singular values and generalized singular vectors. Themethod is an inner-outer iteration method, where the expansion of the rightsearching subspace forms the inner iterations that approximately solve thecorrection equations involved and the outer iterations extract approximate GSVDcomponents with respect to the subspaces. Some convergence results areestablished for the inner and outer iterations, based on some of whichpractical stopping criteria are designed for the inner iterations. Athick-restart CPF-JDGSVD algorithm with deflation is developed to computeseveral GSVD components. Numerical experiments illustrate the efficiency of thealgorithm.
What problem does this paper attempt to address?