Abstract
An algorithm is presented, which produces highly accurate and automatically verified bounds for the matrix exponential function. Our computational approach involves iterative defect correction, interval analysis and advanced computer arithmetic. The algorithm presented is based on the “scaling and squaring” scheme, utilizing Padé approximations and safe error monitoring. A PASCAL-SC program is reported and numerical results are discussed.
Zusammenfassung
Es wird ein Algorithmus vorgestellt, der hochgenaue und automatisch verifizierte Grenzen fuer die Exponentialfunktion einer Matrix liefert. Unser Verfahren benuetzt iterative Defektkorrektur, Intervall-Analysis und eine erweiterte Rechnerarithmetik. Der dargestellte Algorithmus basiert auf dem “scaling and squaring” Schema und benutzt Padé-Approximationen und safe-error-monitoring. Es werden ein PASCAL-SC Programm vorgestellt und numerische Resultate diskutiert.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
ACRITH General Information Manual IBM Publication No. GC33-6163-01 (Version 1, Release 2), 1984.
Auzinger, W., Stetter, H. J.: Accurate arithmetic results for decimal data on non-decimal computers. Computing 35, 141–151, 1985.
Boehmer, K., Hemker, P., Stetter, H.: The Defect Correction Approach. Computing Suppl. 5, 1984, p. 1–32.
HIFICOMP Subroutine library for highly efficient numerical computations. Methodological Guide, ed. S. Markov. Bulgarian Academy of Sciences, Sofia 1987.
Kulisch, U., Miranker, W.: Computer Arithmetic in Theory and Practice. Academic Press, 1981.
Markov, S.: Mathematical fundamentals of numerical computation. Proc. 17th spring conference of UMB “Mathematics and education in mathematics”, Publ. house of the Bulg. Acad. of Sci., Sofia 1988.
Moler, C., van Loan, C.: Nineteen Dubious Ways to Compute The Exponential of a Matrix. SIAM Review,20, 4, 1978, p. 801–836.
Moore, R., E.: Methods and Applications of the Interval Analysis. SIAM, Phyladelfia, 1979.
PASCAL-SC Information Manual and Floppy Disks U. Kulisch—editor, John Willey & Sons, 1987.
Rump, S., Boehm, H.: Least significant Bit Evaluation of Arithmetic Expressions in Single Precision. Computing30, 1983, p. 189–199.
Stetter, H.: Sequential Defect Correction for High-accuracy Floating-point Algorithms. Lecture Notes in Mathematics1006, 1984, p. 186–202.
Ward, R.: Numerical Computation of the Matrix Exponential With Accuracy Estimate. SIAM J. Num. Anal.14, 4, 1977, p. 600–610.
Author information
Authors and Affiliations
Additional information
The present research is partially supported by the Committee of Science, Sofia, Bulgaria, according to contract No. 755/87
Rights and permissions
About this article
Cite this article
Bochev, P., Markov, S. A self-validating numerical method for the matrix exponential. Computing 43, 59–72 (1989). https://doi.org/10.1007/BF02243806
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02243806