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.

Reduced basis approximation of parametric eigenvalue problems in presence of clusters and intersections

Reduced basis approximation of parametric eigenvalue problems in presence of clusters and intersections

On the effect of different samplings to the solution of parametric PDE Eigenvalue Problems

Theory and numerics of subspace approximation of eigenvalue problems

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

Affine approximation of parametrized kernels and model order reduction for nonlocal and fractional Laplace models

Nonlinear compressive reduced basis approximation for multi-parameter elliptic problem

On reduced basis methods for eigenvalue problems, and on its coupling with perturbation theory

Reduced Basis Approximation for Maxwell's Eigenvalue Problem and Parameter-Dependent Domains

Match-based solution of general parametric eigenvalue problems

Reduced basis approximation and a posteriori error estimation: applications to elasticity problems in several parametric settings

On the Approximation of Kernel functions

Reduced basis method for non-symmetric eigenvalue problems: application to the multigroup neutron diffusion equations

Reduced Krylov Basis Methods for Parametric Partial Differential Equations

Reduced basis methods- an application to variational discretization of parametrized elliptic optimal control problems

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

High Order Numerical Methods To Approximate The Singular Value Decomposition

DISCRETE APPROXIMATION OF NON-COMPACT OPERATORS DESCRIBING CONTINUUM-OF-ALLELES MODELS

Semilinear elliptic eigenvalue problem: Parametric analyticity and the uncertainty quantification

Parametric PDEs: Sparse or Low-Rank Approximations?