Abstract
An iteration scheme, for solving the non-linear equations arising in the implementation of implicit Runge-Kutta methods, is proposed. This scheme is particularly suitable for parallel computation and can be applied to any method which has a coefficient matrixA with all eigenvalues real (and positive). For such methods, the efficiency of a modified Newton scheme may often be improved by the use of a similarity transformation ofA but, even when this is the case, the proposed scheme can have advantages for parallel computation. Numerical results illustrate this. The new scheme converges in a finite number of iterations when applied to linear systems of differential equations, achieving this by using the nilpotency of a strictly lower triangular matrixS −1 AS — Λ, with Λ a diagonal matrix. The scheme reduces to the modified Newton scheme whenS −1 AS is diagonal.A convergence result is obtained which is applicable to nonlinear stiff systems.
Zusammenfassung
Wir schlagen ein Iterationsverfahren für die Lösung der nicht-linearen Gleichungen bei der Implementierung von impliziten Runge-Kutta-Verfahren vor, das besonders geeignet für die Parallelverarbeitung ist und immer angewandt werden kann, wenn die KoeffizientenmatrixA nur reelle (und positive) Eigenwerte hat. Für solche Verfahren kann man zwar die Effizienz eines modifizierten Newton-Verfahrens oft auch mit Hilfe einer Ähnlichkeitstransformation vonA verbessern, aber sogar dann hat unser vorgeschlagenes Verfahren Vorteile für die Parallelverarbeitung, wie numerische Beispiele zeigen. Das neue Vorgehen konvergiert für lineare Differentialgleichungssysteme in endlich vielen Schritten wegen der Nilpotenz einer strengen unteren DreiecksmatrixS −1 AS — Λ und reduziert sich auf das modifizierte Newton-Verfahren wennS −1 eine Diagonalmatrix ist. Unser Konvergenzresultat ist auf nichtlineare steife Systeme anwendbar.
Similar content being viewed by others
References
Alexander, R.: The modified Newton method in the solution of stiff ordinary differential equations. Math. Comp.57, 673–701 (1991).
Bickart, T. A.: An efficient solution process for implicit Runge-Kutta methods. SIAM J. Numer. Anal.14, 1022–1027 (1977).
Burrage, K.: A special family of Runge-Kutta methods for solving stiff differential equatons. BIT18, 22–41 (1978).
Burrage, K.: Efficiently implementable algebraically stable Runge-Kutta methods. SIAM J. Numer. Anal.19, 245–258 (1982).
Burrage, K., Butcher, J. C., Chipman, F. H.: An implementation of singly-implicit Runge-Kutta methods. BIT20, 326–340 (1980).
Butcher, J. C.: On the implementation of implicit Runge-Kutta methods. BIT16, 237–240 (1976).
Butcher, J. C.: Some implementation schemes for implicit Runge-Kutta methods. Proc. Dundee Conf. on Numerical Analysis, Lecture notes in Mathematics pp. 12–24. Berlin Heidelberg New York: Springer 1979.
Cash, J. R.: On a class of implicit Runge-Kutta procedures. J. Inst. Math. Appl.19, 455–470 (1977).
Chipman, F. H.: The implementation of Runge-Kutta implicit processes. BIT13, 391–393 (1973).
Collings, A. G., Tee, G. J.: An analysis of Euler and implicit Runge-Kutta numerical integration schemes for structural dynamic problems. Proc. Sixth Australasian Conf. on the Mechanics of Structures and Materials 1977,1, 147–154.
Cooper, G. J.: On the implementation of singly implicit Runge-Kutta methods. Math. Comp.57, 663–672 (1991).
Cooper, G. J.: On the existence of solutions for algebraically stable Runge-Kutta methods. IMA J. Numer. Anal.6, 325–330 (1986).
Cooper, G. J., Butcher, J. C.: An iteration scheme for implicit Runge-Kutta methods. IMA J. Numer. Anal.3, 127–140 (1983).
Cooper, G. J., Sayfy, A.: Semi-explicitA-stable Runge-Kutta methods. Math. Comp.33, 541–556 (1979).
Cooper, G. J., Vignesvaran, R.: A scheme for the implementation of implicit Runge-Kutta methods. Computing45, 321–332 (1990).
Cooper, G. J., Vignesvaran, R.: Some schemes for the implementation of implicit Runge-Kutta methods. J. Comp. Appl. Math.45, 213–225 (1993).
Crouzeix, M., Hundsdorfer, W. H., Spijker, M. N.: On the existence of solutions to the algebraic equations in implicit Runge-Kutta methods. BIT23, 84–91 (1983).
Davis, H. T.: Introduction to non-linear differential and integral equations. New York: Dover 1962.
Dekker, K.: Error bounds for the solution to the algebraic equations in Runge-Kutta methods. BIT24, 347–356 (1984).
Dekker, K., Verwer, J. G.: Stability of Runge-Kutta methods for stiff nonlinear differential equations. Amsterdam. North-Holland, 1984.
Enright, W. H.: Improving the efficiency of matrix operations in the numerical solution of stiff ordinary differential equations. ACM Trans. Math. Software4, 127–136 (1978).
Frank, R., Ueberhuber, C. W.: Iterated defect correction for the efficient solution of stiff systems of ordinary differential equations. BIT17, 146–159 (1977).
Gear, C. W.: The potential for parallelsim in ordinary differential equations. Tech. Rep. UIUC-DCS-R-86-1246, Computer Science Dept., Univ. of Illinois 1986.
Gear, C. W.: Parallel methods in ordinary differential equations, Tech. Rep. UIUC-DCS-R-87-1396. Computer Science Dept., Univ. of Illinois, 1987.
Gear, C. W.: The automatic integration of stiff ordinary differential equations. Proc. IFIP Congress 1968, 187–193.
Hairer, E., Norsett, S. P. Wanner, G.: Solving ordinary differential equations I. Nonstiff problems. Berlin Heidelberg New York Tokyo: Springer 1987.
Hairer, E., Wanner, G.: Algebraically stable and implementalbe Runge-Kutta methods of high order. SIAM J. Numer. Anal.18, 1098–1108 (1981).
van der Houwen, P. J., Sommeijer, B. P.: Iterated Runge-Kutta methods on parallel computers. SIAM J. Sci. Comput.12, 1000–1028 (1991).
Iserles, A., Norsett, S. P.: On the theory of parallel Runge-Kutta methods. IMA J. Numer. Anal.10, 463–488 (1990).
Jackson, K. R., Norsett, S. P.: The potential for parallelism in Runge-Kutta methods. Part 1: RK formulas in standard form, Tech. Rep. 239/90, Computer Science Dept., Univ. of Toronto 1990.
Karakashian, O. A., Rust, W.: On the parallel implementation of implicit Runge-Kutta methods. SIAM J. Sci. Stat. Comput.9, 1085–1090 (1988).
Keeling, S. L.: On implicit Runge-Kutta methods with a stability function having real poles. BIT29, 91–109 (1989).
Kraaijevanger, J. F. B. M., Schneid, J.: On the unique solvability of the Runge-Kutta equations. Numer. Anal. Rep. 90-6, Mathematics Dept., Massachusetts Institute of Technology 1990.
Kreiss, H. O.: Difference methods for stiff ordinary differential equations. SIAM J. Numer. Anal.15, 21–58 (1978).
Norsett, S. P.: Runge-Kutta methods with a multiple real eigenvalue only. BIT16, 388–393 (1976).
Norsett, S. P., Simonsen, H. H.: Aspects of parallel Runge-Kutta methods. In: Numerical methods for ordinary differential equations. Proc. L'Aquila 1987, pp. 103–107. Berlin Heidelberg New York Tokyo: Springer 1989.
Wanner, G., Hairer, E., Norsett, S. P.: Order stars and stability theorems. BIT18, 475–489 (1978).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Cooper, G.J., Vignesvaran, R. On the use of parallel processors for implicit Runge-Kutta methods. Computing 51, 135–150 (1993). https://doi.org/10.1007/BF02243848
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02243848