On CSCS-based iteration method for tempered fractional diffusion equations (Q2396901)

The authors use a CSCS (circulant and skew-circulant splitting) iteration method to solve linear systems arising from the discretization by finite difference methods for time fractional differential equations in one space dimension. The method is shown to be unconditionally convergent and the convergence rate is fast in numerical tests. In each iteration, a circulant system and a skew-circulant system are required to be solved which cost only \(O(N\log N)\) operations by the fast Fourier transform, where \(N\) is the number of interior mesh points in space. Moreover, the induced preconditioner possesses circulant-times-skew-circulant structure so that it can be inverted in \(O(N\log N)\) operation. Some numerical experiments are presented in which the preconditioner performs well with a simple choice of the method parameter.