Sharp estimates for the convergence rate of Orthomin(k) for a class of linear systems
Abstract
In this work we show that the convergence rate of Orthomin($k$) applied to systems of the form $(I+\rho U) x = b$, where $U$ is a unitary operator and $0<\rho<1$, is less than or equal to $\rho$. Moreover, we give examples of operators $U$ and $\rho>0$ for which the asymptotic convergence rate of Orthomin($k$) is exactly $\rho$, thus showing that the estimate is sharp. While the systems under scrutiny may not be of great interest in themselves, their existence shows that, in general, Orthomin($k$) does not converge faster than Orthomin(1). Furthermore, we give examples of systems for which Orthomin($k$) has the same asymptotic convergence rate as Orthomin(2) for $k\ge 2$, but smaller than that of Orthomin(1). The latter systems are related to the numerical solution of certain partial differential equations.
- Publication:
-
arXiv e-prints
- Pub Date:
- September 2011
- DOI:
- 10.48550/arXiv.1109.2669
- arXiv:
- arXiv:1109.2669
- Bibcode:
- 2011arXiv1109.2669D
- Keywords:
-
- Mathematics - Numerical Analysis
- E-Print:
- 20 pages, 6 figures