Instability of the Sherman-Morrison formula and stabilization by iterative refinement
Fuente:
arXiv
Salvato in:
| Autori principali: | , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2025
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866914071573430272 |
|---|---|
| author | Hashemi, Behnam Nakatsukasa, Yuji |
| author_facet | Hashemi, Behnam Nakatsukasa, Yuji |
| contents | Owing to its simplicity and efficiency, the Sherman-Morrison (SM) formula has seen widespread use across various scientific and engineering applications for solving rank-one perturbed linear systems of the form $(A+uv^T)x = b$. Although the formula dates back at least to 1944, its numerical stability properties have remained an open question and continue to be a topic of current research. We analyze the backward stability of the SM, demonstrate its instability in a scenario increasingly common in scientific computing and address an open question posed by Nick Higham on the proportionality of the backward error bound to the condition number of $A$. We then incorporate fixed-precision iterative refinement into the SM framework reusing the previously computed decompositions and prove that, under reasonable assumptions, it achieves backward stability without sacrificing the efficiency of the SM formula. While our theory does not prove the SM formula with iterative refinement always outputs a backward stable solution, empirically it is observed to eventually produce a backward stable solution in all our numerical experiments. We conjecture that with iterative refinement, the SM formula yields a backward stable solution provided that $κ_2(A), κ_2(A+uv^T)$ are both bounded safely away from $ε_M^{-1}$, where $ε_M$ is the unit roundoff. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2510_01696 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Instability of the Sherman-Morrison formula and stabilization by iterative refinement Hashemi, Behnam Nakatsukasa, Yuji Numerical Analysis 65Gxx Owing to its simplicity and efficiency, the Sherman-Morrison (SM) formula has seen widespread use across various scientific and engineering applications for solving rank-one perturbed linear systems of the form $(A+uv^T)x = b$. Although the formula dates back at least to 1944, its numerical stability properties have remained an open question and continue to be a topic of current research. We analyze the backward stability of the SM, demonstrate its instability in a scenario increasingly common in scientific computing and address an open question posed by Nick Higham on the proportionality of the backward error bound to the condition number of $A$. We then incorporate fixed-precision iterative refinement into the SM framework reusing the previously computed decompositions and prove that, under reasonable assumptions, it achieves backward stability without sacrificing the efficiency of the SM formula. While our theory does not prove the SM formula with iterative refinement always outputs a backward stable solution, empirically it is observed to eventually produce a backward stable solution in all our numerical experiments. We conjecture that with iterative refinement, the SM formula yields a backward stable solution provided that $κ_2(A), κ_2(A+uv^T)$ are both bounded safely away from $ε_M^{-1}$, where $ε_M$ is the unit roundoff. |
| title | Instability of the Sherman-Morrison formula and stabilization by iterative refinement |
| topic | Numerical Analysis 65Gxx |
| url | https://arxiv.org/abs/2510.01696 |