A note on grid transfer operators for multigrid methods
Theoretical unification for researchers working on multigrid methods for Toeplitz matrices and PDEs, but the result is incremental as it extends known equivalences.
This paper shows that for elliptic PDEs with constant coefficients, the optimality condition for multigrid methods derived from Local Fourier analysis is equivalent to that from Toeplitz matrix analysis, enabling the use of grid transfer operators with different orders. Numerical experiments confirm the analysis.
The Local Fourier analysis (LFA) is a classic tool to prove convergence theorems for multigrid methods (MGMs). In particular, we are interested in optimality that is a convergence speed independent of the size of the involved matrices. For elliptic partial differential equations (PDEs), a well known optimality result requires that the sum of the orders of the grid transfer operators is not lower than the order of the PDE to solve. Analogously, when dealing with MGMs for Toeplitz matrices in the literature an optimality condition on the position and on the order of the zeros of the symbols of the grid transfer operators has been found. In this work we show that in the case of elliptic PDEs with constant coefficients, the two different approaches lead to an equivalent condition. We argue that the analysis for Toeplitz matrices is an algebraic generalization of the LFA, which allows to deal not only with differential problems but also for instance with integral problems. The equivalence of the two approaches gives the possibility of using grid transfer operators with different orders also for MGMs for Toeplitz matrices. We give also a class of grid transfer operators related to the B-spline's refinement equation and we study their geometric properties. This analysis suggests further links between wavelets and multigrid methods. A numerical experimentation confirms the correctness of the proposed analysis.