Abstract
In this paper, circulant preconditioners are studied for discretized matrices arising from finite difference schemes for a kind of spatial fractional diffusion equations. The fractional differential operator is comprised of left-sided and right-sided derivatives with order in (12,1). The resulting discretized matrices preserve Toeplitz-like structure and hence their matrix-vector multiplications can be computed efficiently by the fast Fourier transform. Theoretically, the spectra of the circulant preconditioned matrices are shown to be clustered around 1 under some conditions. Numerical experiments are presented to demonstrate that the preconditioning technique is very efficient.
Original language | English |
---|---|
Pages (from-to) | 729-747 |
Number of pages | 19 |
Journal | Numerical Algorithms |
Volume | 82 |
Issue number | 2 |
DOIs | |
Publication status | Published - 1 Oct 2019 |
Scopus Subject Areas
- Applied Mathematics
User-Defined Keywords
- Circulant preconditioner
- Fast Fourier transform
- Fractional diffusion equation
- Krylov subspace methods
- Toeplitz matrix