Technical Report CS-2015-02

TR#:CS-2015-02
Class:CS
Title: Minimal Polynomial and Reduced Rank Extrapolation Methods Are Related
Authors: Avram Sidi
PDFCurrently accessibly only within the Technion network
Abstract: Minimal Polynomial Extrapolation (MPE) and Reduced Rank Extrapolation (RRE) are two polynomial methods used for accelerating the convergence of sequences of vectors $\{x_m\}$. They are applied successfully in conjunction with fixed-point iterative schemes in the solution of large and sparse systems of linear and nonlinear equations in different disciplines of science and engineering.

Both methods produce approximations $s_k$ to the limit or antilimit of $\{x_m\}$ that are of the form $s_k=\sum^k_{i=0}\gamma_i x_i$ with $\sum^k_{i=0}\gamma_i=1$, for some scalars $\gamma_i$. The way the two methods are derived suggests that they might, somehow, be related to each other; this has not been explored so far, however.

In this work, we tackle this issue and show that the vectors $s_{k}^\text{MPE}$ and $s_k^\text{RRE}$ produced by the two methods are related in more than one way, and independently of the way the $x_m$ are generated. One of our results states that RRE stagnates, in the sense that $s_k^\text{RRE}=s_{k-1}^\text{RRE}$, if and only if $s_{k}^\text{MPE}$ does not exist. Another result states that, when $s_{k}^\text{MPE}$ exists, there holds $$\mu_k s_k^\text{RRE}=\mu_{k-1} s_{k-1}^\text{RRE}+\nu_k s_{k}^\text{MPE}\quad \text{with}\quad \mu_k=\mu_{k-1}+\nu_k,$$ for some positive scalars $\mu_k$, $\mu_{k-1}$, and $\nu_k$ that depend only on $s_k^\text{RRE}$, $s_{k-1}^\text{RRE}$, and $s_{k}^\text{MPE}$, respectively. Our results are valid when MPE and RRE are defined in any weighted inner product and the norm induced by it. They also contain as special cases the known results pertaining to the connection between the method of Arnoldi and the method of generalized minimal residuals, two important Krylov subspace methods for solving nonsingular linear systems.

CopyrightThe above paper is copyright by the Technion, Author(s), or others. Please contact the author(s) for more information

Remark: Any link to this technical report should be to this page (http://www.cs.technion.ac.il/users/wwwb/cgi-bin/tr-info.cgi/2015/CS/CS-2015-02), rather than to the URL of the PDF files directly. The latter URLs may change without notice.

To the list of the CS technical reports of 2015
To the main CS technical reports page

Computer science department, Technion
admin