- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- BIRS Workshop Lecture Videos /
- Filtered subspace iteration for selfadjoint operator...
Open Collections
BIRS Workshop Lecture Videos
BIRS Workshop Lecture Videos
Filtered subspace iteration for selfadjoint operator eigenvalue problems Ovall, Jeff
Description
Subspace iteration for computing eigenvalues and eigenvectors, a natural generalization of the power method, is among the most straight-forward methods to analyze and implement. A variant of subpace iteration, called FEAST, that uses rational filters to accelerate convergence toward a targeted invariant subspace (typically selected by enclosing the eigenvalues of interest by a simple contour), has gained significant interest in recent years, having been adopted as part of Intel's Math Kernel Library for matrix eigenvalue problems. Such rational filters are often derived from quadrature approximations of contour integrals---think of approximating Cauchy's integral formula for the indicator function of the region enclosed by the contour. The rational filter governs the iterative convergence of the method, so it is perhaps surprising that quadrature error plays little role in the analysis of the algorithm, and we will explain why this is so. Inspired by this approach, we consider filtered subspace iteration for (possibly) unbounded selfadjoint operators, having in mind differential operators as motivational examples. In this broader context, it is natural to consider errors in approximating the invariant subspace with respect to different norms, and we provide a fairly general framework for analyzing both iteration and discretization errors for eigenvalues and invariant subspaces. The bulk of the computational effort in the algorithm involves approximating the action of the resolvent at a few points along a contour enclosing the eigenvalues of interest. In our numerical examples, such finite-rank approximations of the resolvent are obtained by finite element discretizations.
Item Metadata
| Title |
Filtered subspace iteration for selfadjoint operator eigenvalue problems
|
| Creator | |
| Publisher |
Banff International Research Station for Mathematical Innovation and Discovery
|
| Date Issued |
2018-07-03T09:02
|
| Description |
Subspace iteration for computing eigenvalues and eigenvectors, a natural generalization of the power method, is among the most straight-forward methods to analyze and implement. A variant of subpace iteration, called FEAST, that uses rational filters to accelerate convergence toward a targeted invariant subspace (typically selected by enclosing the eigenvalues of interest by a simple contour), has gained significant interest in recent years, having been adopted as part of Intel's Math Kernel Library for matrix eigenvalue problems. Such rational filters are often derived from quadrature approximations of contour integrals---think of approximating Cauchy's integral formula for the indicator function of the region enclosed by the contour. The rational filter governs the iterative convergence of the method, so it is perhaps surprising that quadrature error plays little role in the analysis of the algorithm, and we will explain why this is so. Inspired by this approach, we consider filtered subspace iteration for (possibly) unbounded selfadjoint operators, having in mind differential operators as motivational examples. In this broader context, it is natural to consider errors in approximating the invariant subspace with respect to different norms, and we provide a fairly general framework for analyzing both iteration and discretization errors for eigenvalues and invariant subspaces. The bulk of the computational effort in the algorithm involves approximating the action of the resolvent at a few points along a contour enclosing the eigenvalues of interest. In our numerical examples, such finite-rank approximations of the resolvent are obtained by finite element discretizations.
|
| Extent |
48.0
|
| Subject | |
| Type | |
| File Format |
video/mp4
|
| Language |
eng
|
| Notes |
Author affiliation: Portland State University
|
| Series | |
| Date Available |
2019-03-20
|
| Provider |
Vancouver : University of British Columbia Library
|
| Rights |
Attribution-NonCommercial-NoDerivatives 4.0 International
|
| DOI |
10.14288/1.0377210
|
| URI | |
| Affiliation | |
| Peer Review Status |
Unreviewed
|
| Scholarly Level |
Faculty
|
| Rights URI | |
| Aggregated Source Repository |
DSpace
|
Item Media
Item Citations and Data
Rights
Attribution-NonCommercial-NoDerivatives 4.0 International