Defect estimation in Adams PECE codes

D.J. Higham

Research output: Contribution to journalArticle

Abstract

Many modern codes for solving the nonstiff initial value problem $y'(x) - f(x,y(x)) = 0,y(a)$ given, $a leqq x leqq b$, produce, in addition to a discretised solution, a function $p(x)$ that approximates $y(x)$ over $[a,b]$. The associated defect $delta (x): = p'(x) - f(x,p(x))$ is a natural measure of the error. In this paper the problem of reliably estimating the defect in Adams PECE methods is considered. Attention is focused on the widely used Shampine-Gordon variable order, variable step code fitted with a continuously differentiable interpolant $p(x)$ due to Watts and Shampine [SIAM .J. Sci. Statist. Comput, 7 (1986), pp. 334-345]. It is shown that over each step an asymptotically correct estimate of the defect can be obtained by sampling at a single, suitably chosen point. It is also shown that a valid "free" estimate can be formed without recourse to sampling. Numerical results are given to support the theory.
Original languageEnglish
Pages (from-to)964-976
Number of pages12
JournalSIAM Journal on Scientific Computing
Volume10
Issue number5
DOIs
Publication statusPublished - 1989

Fingerprint

Defects
Sampling
Initial value problems
Continuously differentiable
Interpolants
Estimate
Initial Value Problem
Valid
Numerical Results

Keywords

  • Adams PECE method
  • interpolant
  • defect
  • numerical mathematics

Cite this

Higham, D.J. / Defect estimation in Adams PECE codes. In: SIAM Journal on Scientific Computing. 1989 ; Vol. 10, No. 5. pp. 964-976.
@article{bf42b031d9ae428982ac30dcd11d0e24,
title = "Defect estimation in Adams PECE codes",
abstract = "Many modern codes for solving the nonstiff initial value problem $y'(x) - f(x,y(x)) = 0,y(a)$ given, $a leqq x leqq b$, produce, in addition to a discretised solution, a function $p(x)$ that approximates $y(x)$ over $[a,b]$. The associated defect $delta (x): = p'(x) - f(x,p(x))$ is a natural measure of the error. In this paper the problem of reliably estimating the defect in Adams PECE methods is considered. Attention is focused on the widely used Shampine-Gordon variable order, variable step code fitted with a continuously differentiable interpolant $p(x)$ due to Watts and Shampine [SIAM .J. Sci. Statist. Comput, 7 (1986), pp. 334-345]. It is shown that over each step an asymptotically correct estimate of the defect can be obtained by sampling at a single, suitably chosen point. It is also shown that a valid {"}free{"} estimate can be formed without recourse to sampling. Numerical results are given to support the theory.",
keywords = "Adams PECE method, interpolant, defect, numerical mathematics",
author = "D.J. Higham",
year = "1989",
doi = "10.1137/0910056",
language = "English",
volume = "10",
pages = "964--976",
journal = "SIAM Journal on Scientific Computing",
issn = "1064-8275",
number = "5",

}

Defect estimation in Adams PECE codes. / Higham, D.J.

In: SIAM Journal on Scientific Computing, Vol. 10, No. 5, 1989, p. 964-976.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Defect estimation in Adams PECE codes

AU - Higham, D.J.

PY - 1989

Y1 - 1989

N2 - Many modern codes for solving the nonstiff initial value problem $y'(x) - f(x,y(x)) = 0,y(a)$ given, $a leqq x leqq b$, produce, in addition to a discretised solution, a function $p(x)$ that approximates $y(x)$ over $[a,b]$. The associated defect $delta (x): = p'(x) - f(x,p(x))$ is a natural measure of the error. In this paper the problem of reliably estimating the defect in Adams PECE methods is considered. Attention is focused on the widely used Shampine-Gordon variable order, variable step code fitted with a continuously differentiable interpolant $p(x)$ due to Watts and Shampine [SIAM .J. Sci. Statist. Comput, 7 (1986), pp. 334-345]. It is shown that over each step an asymptotically correct estimate of the defect can be obtained by sampling at a single, suitably chosen point. It is also shown that a valid "free" estimate can be formed without recourse to sampling. Numerical results are given to support the theory.

AB - Many modern codes for solving the nonstiff initial value problem $y'(x) - f(x,y(x)) = 0,y(a)$ given, $a leqq x leqq b$, produce, in addition to a discretised solution, a function $p(x)$ that approximates $y(x)$ over $[a,b]$. The associated defect $delta (x): = p'(x) - f(x,p(x))$ is a natural measure of the error. In this paper the problem of reliably estimating the defect in Adams PECE methods is considered. Attention is focused on the widely used Shampine-Gordon variable order, variable step code fitted with a continuously differentiable interpolant $p(x)$ due to Watts and Shampine [SIAM .J. Sci. Statist. Comput, 7 (1986), pp. 334-345]. It is shown that over each step an asymptotically correct estimate of the defect can be obtained by sampling at a single, suitably chosen point. It is also shown that a valid "free" estimate can be formed without recourse to sampling. Numerical results are given to support the theory.

KW - Adams PECE method

KW - interpolant

KW - defect

KW - numerical mathematics

U2 - 10.1137/0910056

DO - 10.1137/0910056

M3 - Article

VL - 10

SP - 964

EP - 976

JO - SIAM Journal on Scientific Computing

JF - SIAM Journal on Scientific Computing

SN - 1064-8275

IS - 5

ER -