2013/05/28 by Xian‐Ming Gu, Gu, Xian-Ming, Ting‐Zhu Huang +7
Computer Science · Mathematics · Physics and Astronomy · #15A18 #26A33 #65F12 #65L05 #65N22 #Electromagnetic Scattering and Analysis #FOS: Mathematics #Fractional Differential Equations Solutions #Matrix Theory and Algorithms #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.1305.6653
openalex publication_date 2013/05/28 · openalex created_date 2019/06/27 · openalex updated_date 2026/07/28
The finite difference scheme with the shifted Grünwarld formula is employed to semi-discrete the fractional diffusion equations. This spatial discretization can reduce to the large system of ordinary differential equations (ODEs) with initial values. Recently, boundary value method (BVM) was developed as a popular algorithm for solving large systems of ODEs. This method requires the solutions of one or more nonsymmetric, large and sparse linear systems. In this paper, the GMRES method with the block circulant preconditioner is proposed for solving these linear systems. One of the main results is that if an Aν1,ν2-stable boundary value method is used for an m-by-m system of ODEs, then the preconditioner is invertible and the preconditioned matrix can be decomposed as I+L, where I is the identity matrix and the rank of L is at most 2m(ν1+ν2). It means that when the GMRES method is applied to solve the preconditioned linear systems, the method will converge in at most 2m(ν1+ν2)+1 iterations.Finally, extensive numerical experiments are reported to illustrate the effectiveness of our methods for solving the fractional diffusion equations.