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 |
User-Defined Keywords
- Circulant preconditioner
- Fast Fourier transform
- Fractional diffusion equation
- Krylov subspace methods
- Toeplitz matrix