Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
6931200 | Journal of Computational Physics | 2015 | 21 Pages |
Abstract
Numerical methods for fractional differential equations generate full stiffness matrices, which were traditionally solved via Gaussian type direct solvers that require O(N3) of computational work and O(N2) of memory to store where N is the number of spatial grid points in the discretization. We develop a preconditioned fast Krylov subspace iterative method for the efficient and faithful solution of finite volume schemes defined on a locally refined composite mesh for fractional differential equations to resolve boundary layers of the solutions. Numerical results are presented to show the utility of the method.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Jinhong Jia, Hong Wang,