COMMENTS ON HIGH-ORDER INTEGRATORS EMBEDDED WITHIN INTEGRAL DEFERRED CORRECTION METHODS

Spectral deferred correction (SDC) methods for solving ordinary differential equations (ODEs) were introduced by Dutt, Greengard and Rokhlin, [3]. In this paper, we study the properties of these integral deferred correction methods, constructed using high order integrators in the prediction and correction loops, and various distributions of quadrature nodes. The smoothness of the error vector associated with a method, is a key indicator of the order of convergence that can be expected from a scheme, [1, 7, 19]. We will demonstrate that using an rth order method in the correction loop, doesn’t always result in an r-order increase in accuracy. Examples include integral deferred correction methods constructed using non self-starting multi-step integrators, and methods constructed using a non-uniform 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. In brief, known RK methods are more efficient than integral correction methods for low order (< 6th order ) schemes. For high order schemes (≥ 8th order), comparable efficiency is observed numerically. However, it should be noted that integral deferred correction methods offer a much larger region of absolute stability, which may be beneficial in many problems. Additionally, as the order of the embedded integrator is increased, the stability region of these integral deferred correction method increases.

[1]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[2]  Rolf Jeltsch,et al.  Largest disk of stability of explicit Runge-Kutta methods , 1978 .

[3]  L. Greengard,et al.  Spectral Deferred Correction Methods for Ordinary Differential Equations , 2000 .

[4]  Robert D. Skeel,et al.  A Theoretical Framework for Proving Accuracy Results for Deferred Corrections , 1982 .

[5]  Anita T. Layton,et al.  On the choice of correctors for semi-implicit Picard deferred correction methods , 2008 .

[6]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[7]  Yinhua Xia,et al.  Efficient time discretization for local discontinuous Galerkin methods , 2007 .

[8]  A. R. Curtis,et al.  An eighth order Runge-Kutta process with eleven function evaluations per step , 1970 .

[9]  H. A. Luther.  An explicit sixth-order Runge-Kutta formula , 1968 .

[10]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[11]  C.-T. Pan.  Some stability results for explicit Runge-Kutta methods , 1990 .

[12]  M. Minion.  Semi-implicit spectral deferred correction methods for ordinary differential equations , 2003 .

[13]  Michael L. Minion,et al.  Implications of the choice of predictors for semi-implicit Picard integral deferred correction methods , 2007 .

[14]  Thomas Hagstrom,et al.  On the spectral deferred correction of splitting methods for initial value problems , 2006 .

[15]  Michael L. Minion,et al.  Implications of the Choice of Quadrature Nodes for Picard Integral Deferred Corrections Methods for Ordinary Differential Equations , 2005 .

[16]  Jingfang Huang,et al.  Arbitrary order Krylov deferred correction methods for differential algebraic equations , 2007, J. Comput. Phys..

[17]  Jingfang Huang,et al.  Accelerating the convergence of spectral deferred correction methods , 2006, J. Comput. Phys..

[18]  M. E. Hosea.  A new recurrence for computing Runge-Kutta truncation error coefficients , 1995 .

[19]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .