Search

In the present study, we consider the numerical method for Toeplitz-like linear systems arising from the $d$-dimensional Riesz space fractional diffusion equations (RSFDEs). We apply the Crank-Nicolson (CN) technique to discretize the temporal derivative and apply a quasi-compact finite difference method to discretize the Riesz space fractional derivatives. For the $d$-dimensional problem, the corresponding coefficient matrix is the sum of a product of a (block) tridiagonal matrix multiplying a diagonal matrix and a $d$-level Toeplitz matrix. We develop a sine transform based preconditioner to accelerate the convergence of the GMRES method. Theoretical analyses show that the upper bound of relative residual norm of the preconditioned GMRES method with the proposed preconditioner is mesh-independent, which leads to a linear convergence rate. Numerical results are presented to confirm the theoretical results regarding the preconditioned matrix and to illustrate the efficiency of the proposed preconditioner., Comment: 17 pages, 3 figures