$\textrm {PECE}$ algorithms for the solution of stiff systems of ordinary differential equations
HTML articles powered by AMS MathViewer
- by R. W. Klopfenstein and C. B. Davis PDF
- Math. Comp. 25 (1971), 457-473 Request permission
Abstract:
This paper presents a study of a class of PECE algorithms consisting of an application of a predictor followed by application of one iteration of a pseudo Newton-Raphson method to a corrector. Such algorithms require precisely two evaluations of the derivative function for each forward step. Theorems 1 and 4 show that the stability properties of such algorithms compare favorably with those obtained with application of the Newton-Raphson method to the corrector iterated to convergence. A subclass of these algorithms have local truncation error of second order and some have local truncation error of third order. Theorems 2 and 3 exhibit members of this subclass wherein an estimate of the local truncation error is explicit in the algorithm at each step. Initially (in Theorem 1) these algorithms are characterized in terms of their stability properties in the limit as the interval of integration becomes indefinitely large. In Section 5, their properties for other intervals of integration are discussed through the study of some enclosure properties.References
-
F. H. Branin, Jr., βComputer methods of network analysis,β Proc. IEEE, v. 55, 1967, pp. 1787-1801.
D. A. Callahan, βA stable, accurate method of numerical integration for nonlinear systems,β Proc. IEEE, v. 56, 1968, p. 744.
- Germund G. Dahlquist, A special stability problem for linear multistep methods, Nordisk Tidskr. Informationsbehandling (BIT) 3 (1963), 27β43. MR 170477, DOI 10.1007/bf01963532
- C. W. Gear, The automatic integration of stiff ordinary differential equations. , Information Processing 68 (Proc. IFIP Congress, Edinburgh, 1968) North-Holland, Amsterdam, 1969, pp.Β 187β193. MR 0260180
- C. F. Haines, Implicit integration processes with error estimate for the numerical solution of differential equations, Comput. J. 12 (1969/70), 183β187. MR 246513, DOI 10.1093/comjnl/12.2.183
- W. Liniger, A criterion for $A$-stability of linear multistep integration formulae, Computing (Arch. Elektron. Rechnen) 3 (1968), 280β285 (English, with German summary). MR 239763, DOI 10.1007/bf02235394
- Werner Liniger and Ralph A. Willoughby, Efficient integration methods for stiff systems of ordinary differential equations, SIAM J. Numer. Anal. 7 (1970), 47β66. MR 260181, DOI 10.1137/0707002
- Robert D. Richtmyer, Difference methods for initial-value problems, Interscience Tracts in Pure and Applied Mathematics, Tract 4, Interscience Publishers, Inc., New York, 1957. MR 0093918
- H. H. Rosenbrock, Some general implicit processes for the numerical solution of differential equations, Comput. J. 5 (1962/63), 329β330. MR 155434, DOI 10.1093/comjnl/5.4.329
- I. W. Sandberg and H. Shichman, Numerical integration of systems of stiff nonlinear differential equations, Bell System Tech. J. 47 (1968), 511β527. MR 233515, DOI 10.1002/j.1538-7305.1968.tb02483.x S. R. Sedore, SCEPTRE: A Second Generation Transient Analysis Program, Proc. Computer-Aided Circuit Design Seminar, (NASA), Cambridge, Mass., 1967, pp. 55-61.
- Olof B. Widlund, A note on unconditionally stable linear multistep methods, Nordisk Tidskr. Informationsbehandling (BIT) 7 (1967), 65β70. MR 215533, DOI 10.1007/bf01934126
- J. H. Wilkinson, Rounding errors in algebraic processes, Prentice-Hall, Inc., Englewood Cliffs, N.J., 1963. MR 0161456
Additional Information
- © Copyright 1971 American Mathematical Society
- Journal: Math. Comp. 25 (1971), 457-473
- MSC: Primary 65L99
- DOI: https://doi.org/10.1090/S0025-5718-1971-0298956-3
- MathSciNet review: 0298956