Band-Toeplitz preconditioners for ill-conditioned Toeplitz systems

Journal article

Preconditioning for Toeplitz systems has been an active research area over the past few decades. Along this line of research, circulant preconditioners have been recently proposed for the Toeplitz-like system arising from discretizing fractional diffusion equations. A common approach is to combine a circulant preconditioner with the preconditioned conjugate gradient normal residual (PCGRN) method for the coefficient system. In this work, instead of using PCGRN for the normal equation system, we propose a simple yet effective preconditioning approach for solving the original system using the preconditioned minimal residual (PMINRES) method that can achieve convergence guarantees depending only on eigenvalues. Namely, for a large class of ill-conditioned Toeplitz systems, we propose a number of preconditioners that attain the overall O(nlog⁡n) complexity. We first symmetrize the given Toeplitz system by using a permutation matrix and construct a band-Toeplitz plus circulant preconditioner for the modified system. Then, under certain assumptions, we show that the eigenvalues of the preconditioned system are clustered around ±1 except a number of outliers and hence superlinear convergence rate of PMINRES can be achieved. Particularly, we indicate that our solver can be applied to solve certain fractional diffusion equations. An extension of this work to the block Toeplitz case is also included. Numerical examples are provided to demonstrate the effectiveness of our proposed method.

Authors: Hon, S; Serra-Capizzano, S; Wathen, A

• Publication status: Published
• Peer review status: Peer reviewed
• Publisher: Springer
• Journal: BIT Numerical Mathematics
• Volume: 62
• Issue: 2
• Pages: 465–491
• Publication date: 2021-08-09
• Acceptance date: 2021-07-31
• DOI: 10.1007/s10543-021-00889-6
• EISSN: 1572-9125
• ISSN: 0006-3835
• Language: English
• Keywords: fractional diffusion equations; singular value; FFR; Band-Toeplitz; Toeplitz; eigenvalue distribution; block matrices; Hankel matrices; Krylov subspace methods; circulant preconditioners