Abstract:The generalized singular value decomposition (GSVD) of a matrix pair $\{A, L\}$ with $A\in\mathbb{R}^{m\times n}$ and $L\in\mathbb{R}^{p\times n}$ generalizes the singular value decomposition (SVD) of a single matrix. In this paper, we provide a new understanding of GSVD from the viewpoint of SVD, based on which we propose a new iterative method for computing nontrivial GSVD components of a large-scale matrix pair. By introducing two linear operators $\mathcal{A}$ and $\mathcal{L}$ induced by $\{A, L\}$ between two finite-dimensional Hilbert spaces and applying the theory of singular value expansion (SVE) for linear compact operators, we show that the GSVD of $\{A, L\}$ is nothing but the SVEs of $\mathcal{A}$ and $\mathcal{L}$. This result characterizes completely the structure of GSVD for any matrix pair with the same number of columns. As a direct application of this result, we generalize the standard Golub-Kahan bidiagonalization (GKB) that is a basic routine for large-scale SVD computation such that the resulting generalized GKB (gGKB) process can be used to approximate nontrivial extreme GSVD components of $\{A, L\}$, which is named the gGKB\_GSVD algorithm. We use the GSVD of $\{A, L\}$ to study several basic properties of gGKB and also provide preliminary results about convergence and accuracy of gGKB\_GSVD for GSVD computation. Numerical experiments are presented to demonstrate the effectiveness of this method.

The Generalized Triangular Decomposition.

The Geometric Mean Decomposition

Canonical Polyadic Decomposition via the generalized Schur decomposition

Q-less QR Decomposition in Inner Product Spaces

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

Triangular Decomposition of Matrices in a Domain

Generalized cyclic symmetric decompositions for the matrix multiplication tensor

Characterizing GSVD by singular value expansion of linear operators and its computation

Two Harmonic Jacobi–Davidson Methods for Computing a Partial Generalized Singular Value Decomposition of a Large Matrix Pair

Efficient Matrix Decomposition for High-Dimensional Structured Systems: Theory and Applications

The Hadamard decomposition problem

Adaptive Algorithms for Generalized Eigenvalue Decomposition with A Nonquadratic Criterion

A Parallelizable Method for Computing the Singular Subspaces of a General Matrix

Generalized Tensor Function Via the Tensor Singular Value Decomposition Based on the T-product

An Efficient GIPM Algorithm for Computing the Smallest V-Singular Value of the Partially Symmetric Tensor

A Higher-Order Generalized Singular Value Decomposition for Rank Deficient Matrices

Effective Methods of QR-Decompositions of Square Complex Matrices by Fast Discrete Signal-Induced Heap Transforms

Quaternion matrix decomposition and its theoretical implications

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

On Choices of Formulations of Computing the Generalized Singular Value Decomposition of a Large Matrix Pair.

The Joint Bidiagonalization Method for Large GSVD Computations in Finite Precision.