A. H. Bentbib,K. Jbilou,R. Tahiri
Abstract:In this paper, we present a new framework for the recent multidimensional extrapolation methods: Tensor Global Minimal Polynomial (TG-MPE) and Tensor Global Reduced Rank Extrapolation (TG-RRE) methods. We develop a new approach to the one presented in \cite{17}. The proposed framework highlights, in addition their polynomial feature, the connection of TG-MPE and TG-RRE with nonlinear Krylov subspace methods. A unified algorithm is proposed for their implemention. Theoretical results are given and some numerical experiments on linear and nonlinear problems are considered to confirm the performance of the proposed algorithms.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: propose a new framework to reveal the polynomial characteristics of recently developed multi - dimensional extrapolation methods (such as Tensor Global Minimal Polynomial Extrapolation (TG - MPE) and Tensor Global Reduced - Rank Extrapolation (TG - RRE)), and prove their connection with nonlinear Krylov subspace methods. Specifically, these methods aim to accelerate the convergence rate of multi - dimensional sequences, especially for those sequences generated by linear and nonlinear iterative processes.
### Problem Background
In numerical analysis, extrapolation methods are widely used to accelerate the convergence of sequences. However, existing extrapolation methods are usually only applicable to specific types of sequences, and each method has its limitations. Therefore, researchers are constantly exploring new extrapolation methods to adapt to different types of sequences and improve their convergence rates.
### Core Contributions of the Paper
1. **New Framework**: The paper proposes a new framework to understand and interpret the TG - MPE and TG - RRE methods, emphasizing their polynomial characteristics.
2. **Connection with Krylov Subspace Methods**: The paper shows how these methods can be regarded as generalizations of nonlinear Krylov subspace methods when applied to linear and nonlinear problems.
3. **Unified Algorithm**: Provides a unified algorithm to implement these extrapolation methods, making them more conveniently applicable to practical problems.
4. **Theoretical and Experimental Verification**: Through theoretical analysis and numerical experiments, the effectiveness and feasibility of the proposed algorithm are verified.
### Mathematical Formula Representation
To better understand these methods, the following are some key formulas involved in the paper:
- **Einstein Product**:
\[
(A *_{M} B)_{i_{1}, \ldots, i_{N}, k_{1}, \ldots, k_{L}}=\sum_{j_{1}, \ldots, j_{M}} A_{i_{1}, \ldots, i_{N}, j_{1}, \ldots, j_{M}} B_{j_{1}, \ldots, j_{M}, k_{1}, \ldots, k_{L}}
\]
- **Minimal Polynomial**:
\[
\bar{P}_{[E_{n}]}^{M}(M) *_{N} E_{n} = 0
\]
\[
\bar{P}_{[D_{n}]}^{M}(M) *_{N} D_{n} = 0
\]
- **Limit Calculation**:
\[
\bar{X}=\frac{\sum_{i = 0}^{d}\theta_{i}^{(n)} X_{k + i}}{\sum_{i = 0}^{d}\theta_{i}^{(n)}}
\]
Through these formulas, the paper describes in detail how to use the minimal polynomial of tensors to determine the limit of a sequence and shows the potential of these methods in accelerating convergence.