Splitting preconditioning based on sine transform for time-dependent Riesz space fractional diffusion equations

Abstract

We study the sine-transform-based splitting preconditioning technique for the linear systems arising in the numerical discretization of time-dependent one dimensional and two dimensional Riesz space fractional diffusion equations. Those linear systems are Toeplitz-like. By making use of diagonal-plus-Toeplitz splitting iteration technique, a sine-transform-based splitting preconditioner is proposed to accelerate the convergence rate efficiently when the Krylov subspace method is implemented. Theoretically, we prove that the spectrum of the preconditioned matrix of the proposed method is clustering around 1. In practical computations, by the fast sine transform the computational complexity at each time level can be done in $${{\mathcal {O}}}(n\log n)$$ operations where n is the matrix size. Numerical examples are presented to illustrate the effectiveness of the proposed algorithm.