Abstract:Various iterative eigenvalue solvers have been developed to compute parts of the spectrum for a large sparse matrix, including the power method, Krylov subspace methods, contour integral methods, and preconditioned solvers such as the so called LOBPCG method. All of these solvers rely on random matrices to determine, e.g., starting vectors that have, with high probability, a non-negligible overlap with the eigenvectors of interest. For this purpose, a safe and common choice are unstructured Gaussian random matrices. In this work, we investigate the use of random Khatri-Rao products in eigenvalue solvers. On the one hand, we establish a novel subspace embedding property that provides theoretical justification for the use of such structured random matrices. On the other hand, we highlight the potential algorithmic benefits when solving eigenvalue problems with Kronecker product structure, as they arise frequently from the discretization of eigenvalue problems for differential operators on tensor product domains. In particular, we consider the use of random Khatri-Rao products within a contour integral method and LOBPCG. Numerical experiments indicate that the gains for the contour integral method strongly depend on the ability to efficiently and accurately solve (shifted) matrix equations with low-rank right-hand side. The flexibility of LOBPCG to directly employ preconditioners makes it easier to benefit from Khatri-Rao product structure, at the expense of having less theoretical justification.

On the application of Laguerre's method to the polynomial eigenvalue problem

Numerical methods for eigenvalues of singular polynomial eigenvalue problems

A THICK-RESTART LANCZOS ALGORITHM WITH POLYNOMIAL FILTERING FOR HERMITIAN EIGENVALUE PROBLEMS

Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

A Ulm-like Method for Inverse Eigenvalue Problems

Computing the Action of the Generating Function of Bernoulli Polynomials on a Matrix with An Application to Non-local Boundary Value Problems

Fast computation of spectral densities for generalized eigenvalue problems

To solving the eigenvalue problem for polynomial matrices of general form

Automatic rational approximation and linearization of nonlinear eigenvalue problems

An integral method for solving nonlinear eigenvalue problems

Root-Finding with Eigen-Solving

Randomized methods for computing joint eigenvalues, with applications to multiparameter eigenvalue problems and root finding

The infinite Lanczos method for symmetric nonlinear eigenvalue problems

A Note on the Backward Errors for Hermite Eigenvalue Problems

Subspace embedding with random Khatri-Rao products and its application to eigensolvers

Spectral Method with the Tensor-Product Nodal Basis for the Steklov Eigenvalue Problem

Restarted Generalized Krylov Subspace Methods for Solving Large-Scale Polynomial Eigenvalue Problems

Successive linear Newton interpolation methods for solving the large-scale nonlinear eigenvalue problems

Advances in Solution of Classical Generalized Eigenvalue Problem

On matrix balancing and eigenvector computation