Spectral deferred correction methods for ordinary differential equations

Alok Dutt, Leslie Greengard, Vladimir Rokhlin

Research output: Contribution to journalArticle

Abstract

We introduce a new class of methods for the Cauchy problem for ordinary differential equations (ODEs). We begin by converting the original ODE into the corresponding Picard equation and apply a deferred correction procedure in the integral formulation, driven by either the explicit or the implicit Euler marching scheme. The approach results in algorithms of essentially arbitrary order accuracy for both non-stiff and stiff problems; their performance is illustrated with several numerical examples. For non-stiff problems, the stability behavior of the obtained explicit schemes is very satisfactory and algorithms with orders between 8 and 20 should be competitive with the best existing ones. In our preliminary experiments with stiff problems, a simple adaptive implementation of the method demonstrates performance comparable to that of a state-of-the-art extrapolation code (at least, at moderate to high precision). Deferred correction methods based on the Picard equation appear to be promising candidates for further investigation.

Original languageEnglish (US)
Pages (from-to)241-266
Number of pages26
JournalBIT Numerical Mathematics
Volume40
Issue number2
StatePublished - Jun 2000

Fingerprint

Deferred Correction
Ordinary differential equations
Stiff Problems
Ordinary differential equation
Extrapolation
Explicit Scheme
Euler
Cauchy Problem
Numerical Examples
Formulation
Experiments
Arbitrary
Demonstrate
Experiment

Keywords

  • Deferred correction
  • Initial value problems
  • Spectral methods
  • Stiffness

ASJC Scopus subject areas

  • Computer Graphics and Computer-Aided Design
  • Software
  • Applied Mathematics
  • Computational Mathematics

Cite this

Spectral deferred correction methods for ordinary differential equations. / Dutt, Alok; Greengard, Leslie; Rokhlin, Vladimir.

In: BIT Numerical Mathematics, Vol. 40, No. 2, 06.2000, p. 241-266.

Research output: Contribution to journalArticle

@article{ecdd148d37dd401ca8d415ae18d3ecf2,
title = "Spectral deferred correction methods for ordinary differential equations",
abstract = "We introduce a new class of methods for the Cauchy problem for ordinary differential equations (ODEs). We begin by converting the original ODE into the corresponding Picard equation and apply a deferred correction procedure in the integral formulation, driven by either the explicit or the implicit Euler marching scheme. The approach results in algorithms of essentially arbitrary order accuracy for both non-stiff and stiff problems; their performance is illustrated with several numerical examples. For non-stiff problems, the stability behavior of the obtained explicit schemes is very satisfactory and algorithms with orders between 8 and 20 should be competitive with the best existing ones. In our preliminary experiments with stiff problems, a simple adaptive implementation of the method demonstrates performance comparable to that of a state-of-the-art extrapolation code (at least, at moderate to high precision). Deferred correction methods based on the Picard equation appear to be promising candidates for further investigation.",
keywords = "Deferred correction, Initial value problems, Spectral methods, Stiffness",
author = "Alok Dutt and Leslie Greengard and Vladimir Rokhlin",
year = "2000",
month = "6",
language = "English (US)",
volume = "40",
pages = "241--266",
journal = "BIT Numerical Mathematics",
issn = "0006-3835",
publisher = "Springer Netherlands",
number = "2",

}

TY - JOUR

T1 - Spectral deferred correction methods for ordinary differential equations

AU - Dutt, Alok

AU - Greengard, Leslie

AU - Rokhlin, Vladimir

PY - 2000/6

Y1 - 2000/6

N2 - We introduce a new class of methods for the Cauchy problem for ordinary differential equations (ODEs). We begin by converting the original ODE into the corresponding Picard equation and apply a deferred correction procedure in the integral formulation, driven by either the explicit or the implicit Euler marching scheme. The approach results in algorithms of essentially arbitrary order accuracy for both non-stiff and stiff problems; their performance is illustrated with several numerical examples. For non-stiff problems, the stability behavior of the obtained explicit schemes is very satisfactory and algorithms with orders between 8 and 20 should be competitive with the best existing ones. In our preliminary experiments with stiff problems, a simple adaptive implementation of the method demonstrates performance comparable to that of a state-of-the-art extrapolation code (at least, at moderate to high precision). Deferred correction methods based on the Picard equation appear to be promising candidates for further investigation.

AB - We introduce a new class of methods for the Cauchy problem for ordinary differential equations (ODEs). We begin by converting the original ODE into the corresponding Picard equation and apply a deferred correction procedure in the integral formulation, driven by either the explicit or the implicit Euler marching scheme. The approach results in algorithms of essentially arbitrary order accuracy for both non-stiff and stiff problems; their performance is illustrated with several numerical examples. For non-stiff problems, the stability behavior of the obtained explicit schemes is very satisfactory and algorithms with orders between 8 and 20 should be competitive with the best existing ones. In our preliminary experiments with stiff problems, a simple adaptive implementation of the method demonstrates performance comparable to that of a state-of-the-art extrapolation code (at least, at moderate to high precision). Deferred correction methods based on the Picard equation appear to be promising candidates for further investigation.

KW - Deferred correction

KW - Initial value problems

KW - Spectral methods

KW - Stiffness

UR - http://www.scopus.com/inward/record.url?scp=0042185232&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0042185232&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0042185232

VL - 40

SP - 241

EP - 266

JO - BIT Numerical Mathematics

JF - BIT Numerical Mathematics

SN - 0006-3835

IS - 2

ER -