Abstract:We consider the approximation of the smallest eigenvalue of a large parameter-dependent Hermitian matrix over a continuum compact domain. Our approach is based on approximating the smallest eigenvalue by the one obtained by projecting the large matrix onto a suitable small subspace, a practice widely employed in the literature. The projection subspaces are constructed iteratively (to reduce the error of the approximation where it is large) with the addition of the eigenvectors of the parameter-dependent matrix at the parameter values where a surrogate error is maximal. The surrogate error is the gap between the approximation and a lower bound for the smallest eigenvalue proposed in [Sirkovic and Kressner, SIAM J. Matrix Anal. Appl., 37(2), 2016]. Unlike the classical approaches, such as the successive constraint method, that maximize such surrogate errors over a discrete and finite set, we maximize the surrogate error over the continuum of all permissible parameter values globally. We put particular attention to the lower bound, which enables us to formally prove the global convergence of our framework both in finite-dimensional and infinite-dimensional settings. In the second part, we focus on the approximation of the smallest singular value of a large parameter-dependent matrix, in case it is non-Hermitian, and propose another subspace framework to construct a small parameter-dependent non-Hermitian matrix whose smallest singular value approximates the original large-scale smallest singular value. We perform numerical experiments on synthetic examples, as well as on real examples arising from parametric PDEs. The numerical experiments show that the proposed techniques are able to drastically reduce the size of the large parameter-dependent matrix, while ensuring an approximation error for the smallest eigenvalue/singular value below the prescribed tolerance.

Theory and numerics of subspace approximation of eigenvalue problems

A subspace approximation method for the quadratic eigenvalue problem

Subspace method based on neural networks for eigenvalue problems

Fast Calculation of Laplace-Beltrami Eigenproblems Via Subdivision Linear Subspace

Approximating matrix eigenvalues by subspace iteration with repeated random sparsification

Advances in Classical Subspace Iteration Method for Eigenvalue Problem

A Parallel Augmented Subspace Method for Eigenvalue Problems.

Uniform Approximation of Eigenproblems of a Large-Scale Parameter-Dependent Hermitian Matrix

Numerical Approximations of a Nonlinear Eigenvalue Problem and Applications to a Density Functional Model

Subspace iterative methods for eigenvalue problems

Numerical Analysis of the Parallel Orbital-Updating Approach for Eigenvalue Problems

Advances in Solution of Classical Generalized Eigenvalue Problem

Subspace Methods with Local Refinements for Eigenvalue Computation Using Low-Rank Tensor-Train Format.

A numerical projection technique for large-scale eigenvalue problems

Sublinear Time Eigenvalue Approximation via Random Sampling

An Enhanced Formulation Of Error Bound In Subspace Iteration Method

On an Eigenvector-Dependent Nonlinear Eigenvalue Problem.

A theory of quantum subspace diagonalization

A numerical approach for investigating multiple eigenpairs of a quasilinear elliptic system

Gradient-type subspace iteration methods for the symmetric eigenvalue problem

Estimating the numerical range with a Krylov subspace