Picture of person typing on laptop with programming code visible on the laptop screen

World class computing and information science research at Strathclyde...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by University of Strathclyde researchers, including by researchers from the Department of Computer & Information Sciences involved in mathematically structured programming, similarity and metric search, computer security, software systems, combinatronics and digital health.

The Department also includes the iSchool Research Group, which performs leading research into socio-technical phenomena and topics such as information retrieval and information seeking behaviour.

Explore

Error control for initial value problems with discontinuities and delays

Higham, D.J. (1993) Error control for initial value problems with discontinuities and delays. Applied Numerical Mathematics, 12 (4). pp. 315-330. ISSN 0168-9274

Full text not available in this repository. Request a copy from the Strathclyde author

Abstract

When using software for ordinary differential equation (ODE) initial value problems, it is not unreasonable to expect the global error to decrease linearly with the user-supplied error tolerance. For standard ODEs, conditions on an algorithm that guarantee such 'tolerance proportionality' asymptotically (as the error tolerance tends to zero) were derived by Stetter. Here we extend the analysis to cover a certain class of ODEs with low-order derivative discontinuities, and the class of ODEs with constant delays. We show that standard error control techniques will be successful if discontinuities are handled correctly and delay terms are calculated with sufficiently accurate interpolants. It is perhaps surprising that several delay ODE algorithms that have been proposed do not use sufficiently accurate interpolants to guarantee asymptotic proportionality. Our theoretical results are illustrated numerically.