Mixed precision iterative refinement for least squares with linear equality constraints and generalized least squares problems
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2024
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866912575037374464 |
|---|---|
| author | Gao, Bowen Ma, Yuxin Shao, Meiyue |
| author_facet | Gao, Bowen Ma, Yuxin Shao, Meiyue |
| contents | Recent development on mixed precision techniques has largely enhanced the performance of various linear algebra solvers, one of which being the solver for the least squares problem $\min_{x}\lVert b-Ax\rVert_{2}$. By transforming least squares problems into augmented linear systems, mixed precision techniques are capable of refining the lower precision solution to the working precision. In this paper, we propose mixed precision iterative refinement algorithms for two variants of least squares problems -- the least squares problem with linear equality constraints (LSE) and the generalized least squares problem (GLS). Both classical and GMRES-based iterative refinement can be applied to augmented systems of these two problems to improve the accuracy of the solution. For reasonably well-conditioned problems, our algorithms reduce the execution time by a factor of 40% on average compared to the fixed precision ones from LAPACK on the x86-64 architecture. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2406_16499 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Mixed precision iterative refinement for least squares with linear equality constraints and generalized least squares problems Gao, Bowen Ma, Yuxin Shao, Meiyue Numerical Analysis 65F05, 65F08, 65F10 Recent development on mixed precision techniques has largely enhanced the performance of various linear algebra solvers, one of which being the solver for the least squares problem $\min_{x}\lVert b-Ax\rVert_{2}$. By transforming least squares problems into augmented linear systems, mixed precision techniques are capable of refining the lower precision solution to the working precision. In this paper, we propose mixed precision iterative refinement algorithms for two variants of least squares problems -- the least squares problem with linear equality constraints (LSE) and the generalized least squares problem (GLS). Both classical and GMRES-based iterative refinement can be applied to augmented systems of these two problems to improve the accuracy of the solution. For reasonably well-conditioned problems, our algorithms reduce the execution time by a factor of 40% on average compared to the fixed precision ones from LAPACK on the x86-64 architecture. |
| title | Mixed precision iterative refinement for least squares with linear equality constraints and generalized least squares problems |
| topic | Numerical Analysis 65F05, 65F08, 65F10 |
| url | https://arxiv.org/abs/2406.16499 |