This work was partially supported by AFOSR grant FA9550-17-1-0090, Croatian Science Foundation grant HRZZ-9345, bilateral Croatian-USA grant (administered jointly by Croatian-MZO and NSF) and NSF grant DMS-1414365. The numerical studies were facilitated by the equipment acquired using NSF's Major Research Instrumentation grant DMS-1624776.
Numerical analysis, Galerkin methods, Discontinuous functions
We consider the problem of computing a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. A rational function of the operator is constructed such that the eigenspace of interest is its dominant eigenspace, and a subspace iteration procedure is used to approximate this eigenspace. The computed space is then used to obtain approximations of the eigenvalues of interest. An eigenvalue and eigenspace convergence analysis that considers both iteration error and dis- cretization error is provided. A realization of the proposed approach for a model second-order elliptic operator is based on a discontinuous Petrov-Galerkin discretization of the resolvent, and a variety of numerical experiments illustrate its performance.
Gopalakrishnan, Jay; Grubišić, Luka; and Ovall, Jeffrey S., "Filtered Subspace Iteration for Selfadjoint Operators" (2017). Portland Institute for Computational Science Publications. 3.