期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 214, 期 2, 页码 633-656出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2005.10.004
关键词
spectral deferred correction methods; stiff equations; Krylov subspace methods; GMRES
In the recent paper by Dutt, Greengard and Rokhlin, a variant of deferred or defect correction methods is presented which couples Gaussian quadrature with the Picard integral equation formulation of the initial value ordinary differential equation. The resulting spectral deferred correction (SDC) methods have been shown to possess favorable accuracy and stability properties even for versions with very high order of accuracy. In this paper, we show that for linear problems, the iterations in the SDC algorithm are equivalent to constructing a preconditioned Neumann series expansion for the solution of the standard collocation discretization of the ODE. This observation is used to accelerate the convergence of SDC using the GMRES Krylov subspace method. For nonlinear problems, the GMRES acceleration is coupled with a linear implicit approach. Stability and accuracy analyses show the accelerated scheme provides an improvement in the accuracy, efficiency, and stability of the original SDC approach. Furthermore, preliminary numerical experiments show that accelerating the convergence of SDC methods can effectively eliminate the order reduction previously observed for stiff ODE systems. (c) 2005 Elsevier Inc. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据