An approximate inverse preconditioner for spatial fractional diffusion equations with piecewise continuous coefficients

• Published in: International journal of computer mathematics Vol. 97; no. 3; pp. 523 - 545
• Main Authors: Fang, Zhi-Wei, Sun, Hai-Wei, Wei, Hui-Qin
• Format: Journal Article
• Language: English
• Published: Abingdon Taylor & Francis 03.03.2020; Taylor & Francis Ltd
• Subjects: Algorithms; approximate inverse; circulant matrix; Coefficients; Diffusion; Discontinuity; Discretization; Eigenvalues; fast Fourier transform; Finite difference method; Fractional diffusion equation; Interpolation; Krylov subspace methods; Linear systems; piecewise continuous coefficients; Polynomials; Toeplitz matrix
• ISSN: 0020-7160; 1029-0265
• DOI: 10.1080/00207160.2019.1579313

In this paper, we study the discretized linear systems arising from the space-fractional diffusion equations with piecewise continuous coefficients. Using the implicit finite difference scheme with the shifted Grünwald discretization, the resulting linear systems are Toeplitz-like which can be written as the sum of a scaled identity matrix and two diagonal-times-Toeplitz matrices. Standard circulant preconditioners and the existing approximate circulant-inverse preconditioner do not work for such Toeplitz-like linear systems since the discontinuous diffusion coefficients cannot be well approximated by interpolation polynomials. The main aim of this paper is to propose a new approximate circulant-inverse preconditioner to handle the fractional diffusion equations when the diffusion coefficients are piecewise continuous with finite jump discontinuities. Our idea is to approximate the eigenvalues of circulant matrices by the interpolation formula instead of approximating the diffusion coefficients as done by the existing algorithms. Therefore, the discontinuity of the diffusion coefficients does not influence the efficiency of the preconditioner. Theoretically, the spectra of the resulting preconditioned matrices are shown to be clustered around one, which can guarantee the fast convergence rate of the proposed preconditioner. Numerical examples are provided to demonstrate the effectiveness of our method.