Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866909220720345088 |
|---|---|
| author | Palitta, Davide Schweitzer, Marcel Simoncini, Valeria |
| author_facet | Palitta, Davide Schweitzer, Marcel Simoncini, Valeria |
| contents | Thanks to its great potential in reducing both computational cost and memory requirements, combining sketching and Krylov subspace techniques has attracted a lot of attention in the recent literature on projection methods for linear systems, matrix function approximations, and eigenvalue problems. Applying this appealing strategy in the context of linear matrix equations turns out to be far more involved than a straightforward generalization. These difficulties include analyzing well-posedness of the projected problem and deriving possible error estimates depending on the sketching properties. Further computational complications include the lack of a natural residual norm estimate and of an explicit basis for the generated subspace. In this paper we propose a new sketched-and-truncated polynomial Krylov subspace method for Sylvester equations that aims to address all these issues. The potential of our novel approach, in terms of both computational time and storage demand, is illustrated with numerical experiments. Comparisons with a state-of-the-art projection scheme based on rational Krylov subspaces are also included. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2311_16019 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations Palitta, Davide Schweitzer, Marcel Simoncini, Valeria Numerical Analysis 65F45, 68W20, 65F25, 65F50 Thanks to its great potential in reducing both computational cost and memory requirements, combining sketching and Krylov subspace techniques has attracted a lot of attention in the recent literature on projection methods for linear systems, matrix function approximations, and eigenvalue problems. Applying this appealing strategy in the context of linear matrix equations turns out to be far more involved than a straightforward generalization. These difficulties include analyzing well-posedness of the projected problem and deriving possible error estimates depending on the sketching properties. Further computational complications include the lack of a natural residual norm estimate and of an explicit basis for the generated subspace. In this paper we propose a new sketched-and-truncated polynomial Krylov subspace method for Sylvester equations that aims to address all these issues. The potential of our novel approach, in terms of both computational time and storage demand, is illustrated with numerical experiments. Comparisons with a state-of-the-art projection scheme based on rational Krylov subspaces are also included. |
| title | Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations |
| topic | Numerical Analysis 65F45, 68W20, 65F25, 65F50 |
| url | https://arxiv.org/abs/2311.16019 |