Abstract
In [16], Rosenbrock methods of order four are investigated using an embedded method of order three for step size control. Here, we study such a method using Richardson extrapolation for step size control and compare the two techniques with each other. The usual belief that extrapolation is inferior to embedding is not true. Numerical results for the 25 examples of STIFF DETEST and for some more difficult problems show the following behaviour: For low tolerances (∼10−2) embedding is superior, for moderate tolerances (∼10−4) both techniques are comparable and for high tolerances (<10−5) extrapolation is superior. Under certain conditions the extrapolated value can be used for step continuation without stability problems.
Zusammenfassung
In [16] werden Rosenbrockmethoden der Ordnung 4 untersucht, die zur Schrittweitensteuerung eine eingebettete Methode der Ordnung 3 verwenden. In dieser Arbeit werden neue solche Methoden hergeleitet, die auf Schrittweitensteuerung durch Richardsonextrapolation zugeschnitten sind, und beide Techniken miteinander verglichen. Es zeigt sich, daß im Gegensatz zur üblichen Meinung Extrapolation nicht schlechter ist als Einbettung. Die numerischen Ergebnisse für die 25 Beispiele aus STIFF DETEST und einige schwierigere Beispiele zeigen folgendes Verhalten: Für niedrige Genauigkeiten (∼10−2) sind eingebettete Methoden vorteilhaft, für mittlere Genauigkeiten (∼10−4) sind beide Techniken gleichwertig und für hohe Genauigkeiten (<10−5) ist Extrapolation überlegen. Unter bestimmten Voraussetzungen kann der extrapolierte Wert als Startwert für den nächsten Schritt verwendet werden.
Similar content being viewed by others
References
Bader, G., Deuflhard, P.: A semi-implicit mid-point rule for stiff systems of ordinary differential equations. University of Heidelberg, SFB 123, Technical Report Nr. 114 (1981); Numerische Mathematik41, 373–398 (1983).
Bui, T. D.: Some A-stable und L-stable methods for the numerical integration of stiff ordinary differential equations. J. Assoc. Comput. Mach.26, 483–493 (1979).
Bui, T. D.: A note on the Rosenbrock procedure. Math. Comp.33, 971–975 (1979).
Bui, T. D., Poon, S. W. H.: On the computational aspects of Rosenbrock procedures with built-in errors estimates for stiff systems. BIT21, 168–174 (1981).
Cash, J. R.: Semi-implicit Runge-Kutta procedures with error estimates for the numerical integration of stiff systems of ODEs. JACM23, 455–460 (1976).
Chan, Y. N. I., Birnbaum, I., Lapidus, L.: Solution of stiff differential equations and the use of imbedding techniques. Industr. and Eng. Chemistry Fundamentals17, 133–148 (1978).
Enright, W. H., Hull, T. E.: Test results on initial value methods for non-stiff ordinary differential equations. SIAM J. Numer. Anal.13, 944–961 (1976).
Enright, W. H., Hull, T. E., Lindberg, B.: Comparing numerical methods for stiff ordinary differential equations. BIT15, 10–48 (1975).
Gottwald, B. A., Wanner, G.: A reliable Rosenbrock integrator for stiff differential equations. Computing26, 355–360 (1981).
Hairer, E.: Order conditions for numerical methods for partitioned ordinary differential equations. Numer. Math.36, 431–445 (1981).
Hairer, E., Bader, G., Lubich, Ch.: On the stability of semi-implicit methods for ordinary differential equations. BIT22, 211–232 (1982).
Hearn, A. C.: REDUCE 2, Users Manual. University of Utah, Salt Lake City, March 1973.
Hindmarsh, A. C.: GEAR: Ordinary differential equations solver. UCID-30001 Rev. 3, Lawrence Livermore Laboratory, 1974.
Kaps, P.: Rosenbrock-type methods. In: Numerical Methods for Solving Stiff Initial Value Problems. Proceeding, Oberwolfach 28.6.–4.7.1981 (Dahlquist, G., Jeltsch, R., eds.). Bericht Nr. 9, Institut für Geometrie und Praktische Mathematik, RWTH Aachen, Templergraben 55, D-5100 Aachen.
Kaps, P., Poon, S., Bui, T. D.: Rosenbrock methods for stiff ODEs — a comparison of Richardson extrapolation and embedding technique. Institutsnotiz Nr. 1 (1984), Institut für Mathematik und Geometrie, Technikerstrasse 13, A-6020 Innsbruck.
Kaps, P., Rentrop, P.: Generalized Runge-Kutta methods of order four with stepsize control for stiff ordinary differential equations. Numerische Mathematik33, 55–68 (1979).
Kaps, P., Wanner, G.: A study of Rosenbrock-type methods of high order. Numerische Mathematik38, 279–298 (1981).
Lambert, J. D.: Nonlinear methods for stiff systems of ordinary differential equations. In: Conference on the Numerical Solution of Differential Equations, Dundee 1973 (Watson, G. A., ed.). Lecture Notes in Mathematics, Vol. 363 Berlin-Heidelberg-New York: Springer 1974.
Prokopakis, G. J., Seider, W. D.: Adaptive semi-implicit Runge-Kutta method for solution of stiff ordinary differential equations. Industr. and Engin. Chemistry Fundamentals20, 255–266 (1981).
Nørsett, S. P.: C-polynomials for rational approximation to the exponential function. Numer. Mathematik25, 39–56 (1975).
Nørsett, S. P., Wolfbrandt, A.: Order conditions for Rosenbrock type methods. Numer. Math.32, 1–15 (1979).
Robertson, H. H.: Part II. The solution of a set of reaction rate equations. In: Numerical Analysis (Walsh, J., ed.), pp. 178–182. London: Academic Press 1966.
Scott, M. R., Watts, H. A.: In: Numerical Methods for Differential Systems (Lapidus, L., Schiesser, W. E., eds.) pp. 197–227. New York: Academic Press 1976.
Shampine, L. F.: Evaluation of a test set for stiff ODE solvers. ACM Transactions on Mathematical Software7, 409–420 (1981).
Shampine, L. F.: Implementation of Rosenbrock methods. ACM Transactions on Mathematical Software8, 93–113 (1982).
Steihaug, T., Wolfrandt, A.: An attempt to avoid exact Jacobian and nonlinear equations in the numerical solution of stiff differential equations. Math. Comp.33, 521–534 (1979).
Strehmel, K., Weiner, G.: Lokale Fehlerschätzung mittels modifizierter Richardson-Extrapolation in linear expliziten Einschrittverfahren. Computing (to appear).
Verwer, J. G.: Internal S-stability for generalized Runge-Kutta methods. Report NW21/75, Mathematisch Centrum, Amsterdam (1975).
Verwer, J. G.: Instructive experiments with some Runge-Kutta Rosenbrock methods. Computers and Mathematics with Applications8, 217–229 (1982).
Verwer, J. G., Scholz, S.: Rosenbrock methods and time-lagged Jacobian matrices. Report NW 82/80, Mathematisch Centrum, Amsterdam (1980).
Verwer, J. G., Scholz, S., Blom, J. G., Louter-Nool, M.: A class of Runge-Kutta-Rosenbrock methods for solving stiff differential equations. ZAMM63, 13–20 (1983).
Wanner, G.: Private Communication, 1976.
Wanner, G.: On the integration of stiff differential equations. In: Numerical Analysis (Descloux, J., Marti, J. eds.), ISNM, Vol. 37, pp. 209–226. Basel-Stuttgart: Birkhäuser 1977.
Wanner, G.: On the choice of γ for singly-implicit RK or Rosenbrock methods. BIT20, 102–106 (1980).
Willoughby, R. A.: Stiff Differential Systems. New York: Plenum Press 1974.
Author information
Authors and Affiliations
Additional information
This work has been supported in part by the Deutsche Forschungsgemeinschaft and by the National Research Council of Canada.
Rights and permissions
About this article
Cite this article
Kaps, P., Poon, S.W.H. & Bui, T.D. Rosenbrock methods for Stiff ODEs: A comparison of Richardson extrapolation and embedding technique. Computing 34, 17–40 (1985). https://doi.org/10.1007/BF02242171
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02242171