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.

Preconditioning Technique for Symmetric M-matrices

A Note on Preconditioning for M-Matrix

A Generalization of T. Chan’s Preconditioner

A CONSTRAINT PRECONDITIONER FOR SOLVING SYMMETRIC POSITIVE DEFINITE SYSTEMS AND APPLICATION TO THE HELMHOLTZ EQUATIONS AND POISSON EQUATIONS

A Modified SSOR-like Preconditioner for Non-Hermitian Positive Definite Matrices

A preconditioned version of the MBP iteration method for a class of complex symmetric linear systems

Circulant Preconditioners for Solving Differential Equations with Multidelays

Compact quasi‐Newton preconditioners for symmetric positive definite linear systems

CIMGS: an Incomplete Orthogonal FactorizationPreconditioner

Optimal preconditioners for nonsymmetric multilevel Toeplitz systems with application to solving non-local evolutionary partial differential equations

Circulant Preconditioners for Solving Singular Perturbation Delay Differential Equations

Preconditioning techniques for generalized Sylvester matrix equations

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

Practical Restrictively Preconditioned Conjugate Gradient Methods for a Class of Block Two-by-Two Linear Systems

Preconditioned iterative method for nonsymmetric saddle point linear systems

A Preconditioned MINRES Method for Block Lower Triangular Toeplitz Systems

Efficient Preconditioner Updates for Unsymmetric Shifted Linear Systems

Preconditioned AOR Iterative Method for H-matrices

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

Structure Preprocessing Method for the System of Unclosed Linear Algebraic Equations