Polynomial eigenvalue problems are often found in scientific computing applications. When the coefficient matrices of the polynomial are large and sparse, usually only a few eigenpairs are required and projection methods are the best choice. We focus on Krylov methods that operate on the companion linearization of the polynomial but exploit the block structure with the aim of being memory-efficient in the representation of the Krylov subspace basis. The problem may appear in the form of a low-degree polynomial (quartic or quintic, say) expressed in the monomial basis, or a high-degree polynomial (coming from interpolation of a nonlinear eigenproblem) expressed in a nonmonomial basis. We have implemented a parallel solver in SLEPc covering both cases that is able to compute exterior as well as interior eigenvalues via spectral transformation. We discuss important issues such as scaling and restart and illustrate the robustness and performance of the solver with some numerical experiments.
%0 Journal Article
%1 campos2016parallel
%A Campos, Carmen
%A Roman, Jose E.
%D 2016
%I Society for Industrial & Applied Mathematics (SIAM)
%J SIAM Journal on Scientific Computing
%K 65-04-numerical-analysis-software-source-code 65f15-numerical-eigenvalues-eigenvectors 65y05-parallel-computation 65y15-packaged-numerical-methods
%N 5
%P S385--S411
%R 10.1137/15m1022458
%T Parallel Krylov Solvers for the Polynomial Eigenvalue Problem in SLEPc
%U https://epubs.siam.org/doi/10.1137/15M1022458
%V 38
%X Polynomial eigenvalue problems are often found in scientific computing applications. When the coefficient matrices of the polynomial are large and sparse, usually only a few eigenpairs are required and projection methods are the best choice. We focus on Krylov methods that operate on the companion linearization of the polynomial but exploit the block structure with the aim of being memory-efficient in the representation of the Krylov subspace basis. The problem may appear in the form of a low-degree polynomial (quartic or quintic, say) expressed in the monomial basis, or a high-degree polynomial (coming from interpolation of a nonlinear eigenproblem) expressed in a nonmonomial basis. We have implemented a parallel solver in SLEPc covering both cases that is able to compute exterior as well as interior eigenvalues via spectral transformation. We discuss important issues such as scaling and restart and illustrate the robustness and performance of the solver with some numerical experiments.
@article{campos2016parallel,
abstract = {Polynomial eigenvalue problems are often found in scientific computing applications. When the coefficient matrices of the polynomial are large and sparse, usually only a few eigenpairs are required and projection methods are the best choice. We focus on Krylov methods that operate on the companion linearization of the polynomial but exploit the block structure with the aim of being memory-efficient in the representation of the Krylov subspace basis. The problem may appear in the form of a low-degree polynomial (quartic or quintic, say) expressed in the monomial basis, or a high-degree polynomial (coming from interpolation of a nonlinear eigenproblem) expressed in a nonmonomial basis. We have implemented a parallel solver in SLEPc covering both cases that is able to compute exterior as well as interior eigenvalues via spectral transformation. We discuss important issues such as scaling and restart and illustrate the robustness and performance of the solver with some numerical experiments.},
added-at = {2021-07-09T06:05:04.000+0200},
author = {Campos, Carmen and Roman, Jose E.},
biburl = {https://www.bibsonomy.org/bibtex/238a08ff148b7fbb8ac651659458f0216/gdmcbain},
doi = {10.1137/15m1022458},
interhash = {cbbd2bca73708f690c53b0c1509ff57f},
intrahash = {38a08ff148b7fbb8ac651659458f0216},
journal = {{SIAM} Journal on Scientific Computing},
keywords = {65-04-numerical-analysis-software-source-code 65f15-numerical-eigenvalues-eigenvectors 65y05-parallel-computation 65y15-packaged-numerical-methods},
month = jan,
number = 5,
pages = {S385--S411},
publisher = {Society for Industrial {\&} Applied Mathematics ({SIAM})},
timestamp = {2021-07-09T06:05:04.000+0200},
title = {Parallel Krylov Solvers for the Polynomial Eigenvalue Problem in {SLEPc}},
url = {https://epubs.siam.org/doi/10.1137/15M1022458},
volume = 38,
year = 2016
}