A Ritz method for solution of parametric generalized EVPs

This work deals with approximate solution of generalized eigenvalue problem with coefficient matrix that is an affine function of d-parameters. The coefficient matrix is assumed to be symmetric positive definite and spectrally equivalent to an average matrix for all parameters in a given set. We develop a Ritz method for rapidly approximating the eigenvalues on the spectral interval of interest (0, ) for given parameter value. The Ritz subspace is the same for all parameters and it is designed based on the observation that any eigenvector can be split into two components. The first component belongs to a subspace spanned by some eigenvectors of the average matrix. The second component is defined by a correction operator that is a d + 1 dimensional analytic function. We use this structure and build our Ritz subspace from eigenvectors of the average matrix and samples of the correction operator. The samples are evaluated at interpolation points related to a sparse polynomial interpolation method. We show that the resulting Ritz subspace can approximate eigenvectors of the original problem related to the spectral interval of interest with the same accuracy as the sparse polynomial interpolation approximates the correction operator. Bound for Ritz eigenvalue error follows from this and known results. Theoretical results are illustrated by numerical examples. The advantage of our approach is that the analysis treats multiple eigenvalues and eigenvalue crossings that typically have posed technical challenges in similar works.