Convergence rates for adaptive approximation of ordinary differential equations
2003 (English)In: Numerische Mathematik, ISSN 0029-599X, E-ISSN 0945-3245, Vol. 96, no 1, 99-129 p.Article in journal (Refereed) Published
This paper constructs an adaptive algorithm for ordinary differential equations and analyzes its asymptotic behavior as the error tolerance parameter tends to zero. An adaptive algorithm, based on the error indicators and successive subdivision of time steps, is proven to stop with the optimal number, N, of steps up to a problem independent factor defined in the algorithm. A version of the algorithm with decreasing tolerance also stops with the total number of steps, including all refinement levels, bounded by O(N). The alternative version with constant tolerance stops with O(N log N) total steps. The global error is bounded by the tolerance parameter asymptotically as the tolerance tends to zero. For a p-th order accurate method the optimal number of adaptive steps is proportional to the p-th root of the L 1/p+1 quasi-norm of the error density, while the number of uniform steps, with the same error, is proportional to the p-th root of the larger L-1-norm of the error density.
Place, publisher, year, edition, pages
2003. Vol. 96, no 1, 99-129 p.
finite-element solution, time-step control, error estimation
IdentifiersURN: urn:nbn:se:kth:diva-23000DOI: 10.1007/s00211-003-0466-9ISI: 000186956500005OAI: oai:DiVA.org:kth-23000DiVA: diva2:341698
QC 201005252010-08-102010-08-10Bibliographically approved