GPU implementation of Krylov solvers for block-tridiagonal eigenvalue problems
In an eigenvalue problem defined by one or two matrices with block-tridiagonal structure, if only a few eigenpairs are required it is interesting to consider iterative methods based on Krylov subspaces, even if matrix blocks are dense. In this context, using the GPU for the associated dense linear a...
| Authors: | , |
|---|---|
| Format: | book part |
| Publication Date: | 2016 |
| Country: | España |
| Institution: | Universitat Politècnica de València (UPV) |
| Repository: | RiuNet. Repositorio Institucional de la Universitat Politécnica de Valéncia |
| Language: | English |
| OAI Identifier: | oai:riunet.upv.es:10251/66516 |
| Online Access: | https://riunet.upv.es/handle/10251/66516 |
| Access Level: | Open access |
| Keyword: | GPU computing Eigenvalue computation Krylov methods Block-tridiagonal linear solvers CIENCIAS DE LA COMPUTACION E INTELIGENCIA ARTIFICIAL |
| Summary: | In an eigenvalue problem defined by one or two matrices with block-tridiagonal structure, if only a few eigenpairs are required it is interesting to consider iterative methods based on Krylov subspaces, even if matrix blocks are dense. In this context, using the GPU for the associated dense linear algebra may provide high performance. We analyze this in an implementation done in the context of SLEPc, the Scalable Library for Eigenvalue Problem Computations. In the case of a generalized eigenproblem or when interior eigenvalues are computed with shift-and-invert, the main computational kernel is the solution of linear systems with a block-tridiagonal matrix. We explore possible implementations of this operation on the GPU, including a block cyclic reduction algorithm. |
|---|