Abstract:In this work, we develop a novel multilevel Tau matrix-based preconditioned method for a class of non-symmetric multilevel Toeplitz systems. This method not only accounts for but also improves upon an ideal preconditioner pioneered by [J. Pestana. Preconditioners for symmetrized Toeplitz and multilevel Toeplitz matrices. SIAM J. Matrix Anal. Appl., 40(3):870-887, 2019]. The ideal preconditioning approach was primarily examined numerically in that study, and an effective implementation was not included. To address these issues, we first rigorously show in this study that this ideal preconditioner can indeed achieve optimal convergence when employing the MINRES method, with a convergence rate is that independent of the mesh size. Then, building on this preconditioner, we develop a practical and optimal preconditioned MINRES method. To further illustrate its applicability and develop a fast implementation strategy, we consider solving Riemann-Liouville fractional diffusion equations as an application. Specifically, following standard discretization on the equation, the resultant linear system is a non-symmetric multilevel Toeplitz system, affirming the applicability of our preconditioning method. Through a simple symmetrization strategy, we transform the original linear system into a symmetric multilevel Hankel system. Subsequently, we propose a symmetric positive definite multilevel Tau preconditioner for the symmetrized system, which can be efficiently implemented using discrete sine transforms. Theoretically, we demonstrate that mesh-independent convergence can be achieved. In particular, we prove that the eigenvalues of the preconditioned matrix are bounded within disjoint intervals containing $\pm 1$, without any outliers. Numerical examples are provided to critically discuss the results, showcase the spectral distribution, and support the efficacy of our strategy.

The Solution of Ill-Conditioned Symmetric Toeplitz Systems Via Two-Grid and Wavelet Methods

The Modified Solutoin for Large Sparse Symmetric Linear Systems in Electromagnetic Field Analysis

An Aggregation-Based Two-Grid Method for Multilevel Block Toeplitz Linear Systems

Multigrid method for symmetric Toeplitz tridiagonal matrix: Convergence analysis and application to fractional Feynman-Kac equation

Analysis on aggregation and block smoothers in multigrid methods for block Toeplitz linear systems

On HSS and AHSS Iteration Methods for Nonsymmetric Positive Definite Toeplitz Systems.

Numerical algorithms for the fast and reliable solution of periodic tridiagonal Toeplitz linear systems

Unconditional superconvergence analysis of a two-grid finite element method for nonlinear wave equations

On The Hss Iteration Methods For Positive Definite Toeplitz Linear Systems (Vol 224, Pg 709, 2009)

Dct- and Dst-Based Splitting Methods for Toeplitz Systems

Two-Grid Deflated Krylov Methods for Linear Equations

On Τ Matrix-Based Approximate Inverse Preconditioning Technique for Diagonal-Plus-toeplitz Linear Systems from Spatial Fractional Diffusion Equations

Two-grid algorithms for some linear and nonlinear elliptic systems

On K-Step CSCS-based Polynomial Preconditioners for Toeplitz Linear Systems with Application to Fractional Diffusion Equations

A Limited-Memory Block Bi-Diagonal Toeplitz Preconditioner for Block Lower Triangular Toeplitz System from Time-Space Fractional Diffusion Equation

Diagonal and normal with Toeplitz-block splitting iteration method for space fractional coupled nonlinear Schrödinger equations with repulsive nonlinearities

An approximate inverse preconditioner for Toeplitz systems with multiple right-hand sides.

A note on the convergence of multigrid methods for the Riesz-space equation and an application to image deblurring

An Efficient Iteration Method for Toeplitz-Plus-Band Triangular Systems Generated from Fractional Ordinary Differential Equation

Multilevel Tau preconditioners for symmetrized multilevel Toeplitz systems with applications to solving space fractional diffusion equations

New splitting iteration methods for toeplitz-plus-hankel systems