A Jacobi–Davidson type method for the generalized singular value problem

We discuss a new method for the iterative computation of some of the generalized singular values and vectors of a large sparse matrix. Our starting point is the augmented matrix formulation of the GSVD. The subspace expansion is performed by (approximately) solving a Jacobi–Davidson type correction equation, while we give several alternatives for the subspace extraction. Numerical experiments illustrate the performance of the method.