Polynomial Preconditioning for the Action of the Matrix Square Root and Inverse Square Root
Fuente:
arXiv
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2024
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866914639806201856 |
|---|---|
| author | Frommer, Andreas Ramirez-Hidalgo, Gustavo Schweitzer, Marcel Tsolakis, Manuel |
| author_facet | Frommer, Andreas Ramirez-Hidalgo, Gustavo Schweitzer, Marcel Tsolakis, Manuel |
| contents | While preconditioning is a long-standing concept to accelerate iterative methods for linear systems, generalizations to matrix functions are still in their infancy. We go a further step in this direction, introducing polynomial preconditioning for Krylov subspace methods which approximate the action of the matrix square root and inverse square root on a vector. Preconditioning reduces the subspace size and therefore avoids the storage problem together with -- for non-Hermitian matrices -- the increased computational cost per iteration that arises in the unpreconditioned case. Polynomial preconditioning is an attractive alternative to current restarting or sketching approaches since it is simpler and computationally more efficient. We demonstrate this for several numerical examples. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2401_06684 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Polynomial Preconditioning for the Action of the Matrix Square Root and Inverse Square Root Frommer, Andreas Ramirez-Hidalgo, Gustavo Schweitzer, Marcel Tsolakis, Manuel Numerical Analysis 65F60, 65F08, 65F50, 15A16 While preconditioning is a long-standing concept to accelerate iterative methods for linear systems, generalizations to matrix functions are still in their infancy. We go a further step in this direction, introducing polynomial preconditioning for Krylov subspace methods which approximate the action of the matrix square root and inverse square root on a vector. Preconditioning reduces the subspace size and therefore avoids the storage problem together with -- for non-Hermitian matrices -- the increased computational cost per iteration that arises in the unpreconditioned case. Polynomial preconditioning is an attractive alternative to current restarting or sketching approaches since it is simpler and computationally more efficient. We demonstrate this for several numerical examples. |
| title | Polynomial Preconditioning for the Action of the Matrix Square Root and Inverse Square Root |
| topic | Numerical Analysis 65F60, 65F08, 65F50, 15A16 |
| url | https://arxiv.org/abs/2401.06684 |