Fast Flexible LSQR with a Hybrid Variant for Efficient Large-Scale Regularization
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866908419609329664 |
|---|---|
| author | Mikušová, Eva Hnětynková, Iveta |
| author_facet | Mikušová, Eva Hnětynková, Iveta |
| contents | A wide range of applications necessitates solving large-scale ill-posed problems contaminated by noise. Krylov subspace regularization methods are particularly advantageous in this context, as they rely solely on matrix-vector multiplication. Among the most widely used techniques are LSQR and CGLS, both of which can be extended with flexible preconditioning to enforce solution properties such as nonnegativity or sparsity. Flexible LSQR (FLSQR) can also be further combined with direct methods to create efficient hybrid approaches. The Flexible Golub-Kahan bidiagonalization underlying FLSQR requires two long-term recurrences. In this paper, we introduce a novel Fast Flexible Golub-Kahan bidiagonalization method that employs one long-term and one short-term recurrence. Using this, we develop the Fast Flexible LSQR (FaFLSQR) algorithm, which offers comparable computational cost to FCGLS while also supporting hybrid regularization like FLSQR. We analyze the properties of FaFLSQR and prove its mathematical equivalence to FCGLS. Numerical experiments demonstrate that in floating point arithmetic FaFLSQR outperforms both FCGLS and FLSQR in terms of computational efficiency. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2506_19666 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Fast Flexible LSQR with a Hybrid Variant for Efficient Large-Scale Regularization Mikušová, Eva Hnětynková, Iveta Numerical Analysis 65F22, 65F08, 65F10, 15A29 A wide range of applications necessitates solving large-scale ill-posed problems contaminated by noise. Krylov subspace regularization methods are particularly advantageous in this context, as they rely solely on matrix-vector multiplication. Among the most widely used techniques are LSQR and CGLS, both of which can be extended with flexible preconditioning to enforce solution properties such as nonnegativity or sparsity. Flexible LSQR (FLSQR) can also be further combined with direct methods to create efficient hybrid approaches. The Flexible Golub-Kahan bidiagonalization underlying FLSQR requires two long-term recurrences. In this paper, we introduce a novel Fast Flexible Golub-Kahan bidiagonalization method that employs one long-term and one short-term recurrence. Using this, we develop the Fast Flexible LSQR (FaFLSQR) algorithm, which offers comparable computational cost to FCGLS while also supporting hybrid regularization like FLSQR. We analyze the properties of FaFLSQR and prove its mathematical equivalence to FCGLS. Numerical experiments demonstrate that in floating point arithmetic FaFLSQR outperforms both FCGLS and FLSQR in terms of computational efficiency. |
| title | Fast Flexible LSQR with a Hybrid Variant for Efficient Large-Scale Regularization |
| topic | Numerical Analysis 65F22, 65F08, 65F10, 15A29 |
| url | https://arxiv.org/abs/2506.19666 |