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.

Block Splitting Preconditioner for Time-Space Fractional Diffusion Equations

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

On Structure Preserving and Circulant Preconditioners for the Space Fractional Coupled Nonlinear Schrödinger Equations

A Class of Preconditioners Based on Positive-Definite Operator Splitting Iteration Methods for Variable-Coefficient Space-Fractional Diffusion Equations

Preconditioning Techniques for Diagonal-times-Toeplitz Matrices in Fractional Diffusion Equations

On Preconditioners Based on HSS for the Space Fractional CNLS Equations

A -Preconditioner for Space Fractional Diffusion Equation with Non-separable Variable Coefficients

A Preconditioned MINRES Method for Block Lower Triangular Toeplitz Systems

Quasi‐Toeplitz Splitting Iteration Methods for Unsteady Space‐fractional Diffusion Equations

An Efficient Preconditioner for Linear Systems Arising from High-Order Accurate Schemes of Time Fractional Diffusion Equations

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

Circulant Preconditioned Iterations for Fractional Diffusion Equations Based on Hermitian and Skew-Hermitian Splittings

Preconditioned Modified Hermitian and Skew-Hermitian Splitting Iteration Methods for Fractional Nonlinear Schrödinger Equations

Approximate Inverse Preconditioners for Linear Systems Arising from Spatial Balanced Fractional Diffusion Equations

A $τ$-preconditioner for space fractional diffusion equation with non-separable variable coefficients

Diagonal and Toeplitz Splitting Iteration Methods for Diagonal-Plus-toeplitz Linear Systems from Spatial Fractional Diffusion Equations.

Fast Iterative Solvers for Linear Systems Arising from Time-Dependent Space-Fractional Diffusion Equations.

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

A Twisted Block Tangential Filtering Decomposition Preconditioner

A τ Matrix Based Approximate Inverse Preconditioning for Tempered Fractional Diffusion Equations