期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 293, 期 -, 页码 359-369出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2014.08.021
关键词
Anomalous diffusion; Circulant matrix; Conjugate gradient squared method; Fast Fourier transform; Neumann boundary condition; Space-fractional diffusion equation; Toeplitz matrix
资金
- National Science Foundation [EAR-0934747, DMS-1216923]
- National Natural Science Foundation of China [91130010, 11471194]
- State Scholarship Fund from China Scholarship Council
- Direct For Mathematical & Physical Scien
- Division Of Mathematical Sciences [1216923] Funding Source: National Science Foundation
Numerical methods for space-fractional diffusion equations often generate dense or even full stiffness matrices. Traditionally, these methods were solved via Gaussian type direct solvers, which requires O(N-3) of computational work per time step and O(N-2) of memory to store where N is the number of spatial grid points in the discretization. In this paper we develop a preconditioned fast Krylov subspace iterative method for the efficient and faithful solution of finite difference methods (both steady-state and time-dependent) space-fractional diffusion equations with fractional derivative boundary conditions in one space dimension. The method requires O(N) of memory and O(NlogN) of operations per iteration. Due to the application of effective preconditioners, significantly reduced numbers of iterations were achieved that further reduces the computational cost of the fast method. Numerical results are presented to show the utility of the method. (C) 2014 Elsevier Inc. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据