Communications in Applied Mathematics and Computational Science
- Commun. Appl. Math. Comput. Sci.
- Volume 4, Number 1 (2009), 27-56.
Comments on high-order integrators embedded within integral deferred correction methods
A class of novel deferred correction methods, integral deferred correction (IDC) methods, is studied. This class of methods is an extension of ideas introduced by Dutt, Greengard and Rokhlin on spectral deferred correction (SDC) methods for solving ordinary differential equations (ODEs). The novel nature of this class of defect correction methods is that the correction of the defect is carried out using an accurate integral form of the residual instead of the more familiar differential form. As a family of methods, these schemes are capable of matching the efficiency of popular high-order RK methods.
The smoothness of the error vector associated with an IDC method is an important indicator of the order of convergence that can be expected from a scheme (Christlieb, Ong, and Qiu; Hansen and Strain; Skeel). It is demonstrated that embedding an -th order integrator in the correction loop of an IDC method does not always result in an -th order increase in accuracy. Examples include IDC methods constructed using non-self-starting multistep integrators, and IDC methods constructed using a nonuniform distribution of quadrature nodes.
Additionally, the integral deferred correction concept is reposed as a framework to generate high-order Runge–Kutta (RK) methods; specifically, we explain how the prediction and correction loops can be incorporated as stages of a high-order RK method. This alternate point of view allows us to utilize standard methods for quantifying the performance (efficiency, accuracy and stability) of integral deferred correction schemes. It is found that IDC schemes constructed using uniformly distributed nodes and high-order integrators are competitive in terms of efficiency with IDC schemes constructed using Gauss–Lobatto nodes and forward Euler integrators. With respect to regions of absolute stability, however, IDC methods constructed with uniformly distributed nodes and high-order integrators are far superior. It is observed that as the order of the embedded integrator increases, the stability region of the IDC method increases.
Commun. Appl. Math. Comput. Sci., Volume 4, Number 1 (2009), 27-56.
Received: 13 November 2008
Revised: 20 January 2009
Accepted: 15 March 2009
First available in Project Euclid: 20 December 2017
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Primary: 65L05: Initial value problems 65L06: Multistep, Runge-Kutta and extrapolation methods 65L20: Stability and convergence of numerical methods 65L70: Error bounds
Secondary: 65B05: Extrapolation to the limit, deferred corrections
Christlieb, Andrew; Ong, Benjamin; Qiu, Jing-Mei. Comments on high-order integrators embedded within integral deferred correction methods. Commun. Appl. Math. Comput. Sci. 4 (2009), no. 1, 27--56. doi:10.2140/camcos.2009.4.27. https://projecteuclid.org/euclid.camcos/1513798576