A PRECONDITIONED FAST FINITE DIFFERENCE METHOD FOR SPACE-TIME FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS

We develop a fast space-time finite difference method for space-time fractional diffusion equations by fully utilizing the mathematical structure of the scheme. A circulant block preconditioner is proposed to further reduce the computational costs. The method has optimal-order memory requirement and approximately linear computational complexity. The method is not lossy, as no compression of the underlying numerical scheme has been employed. Consequently, the method retains the stability, accuracy, and, in particular, the conservation property of the underlying numerical scheme. Numerical experiments are presented to show the efficiency and capacity of long time modelling of the new method.