Abstract
It is explained, how Householder reflections can be replaced by fast Givens rotations in the first and second step of the QZ algorithm. A count of the required operations for the two steps shows, that the modified first step does need less operations, but that the modified second step needs only insignificantly less operations in comparison with the original QZ algorithm. On the base of this fact a modification of the second step is proposed. The resulting new algorithm is called FQZ algorithm, which has the same numerical properties as the QZ algorithm, if a correct scaling is applied.
Zusammenfassung
Es wird dargestellt, auf welche Weise im ersten und zweiten Schritt des QZ-Algorithmus Householder-Transformationen durch schnelle Givens-Rotationen ersetzt werden können. Die Bestimmung des Rechenaufwandes für diese beiden Schritte zeigt, daß sich im ersten Schritt gegenüber dem ursprünglichen Algorithmus Operationen einsparen lassen, im zweiten Schritt jedoch nur eine unwesentliche Einsparung möglich ist. Auf Grund dieser Tatsache wird der zweite Schritt modifiziert und der daraus enstehende neue Algorithmus als FQZ-Algorithmus bezeichnet. Falls korrekt skaliert wird, entsprechen die numerischen Eigenschaften des FQZ-Algorithmus denjenigen des originalen QZ-Algorithmus.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Bunse, W., Bunse-Gerstner, A.: Numerische lineare Algebra. Stuttgart: B. G. Teubner 1985.
Gentleman, W. M.: Least squares computations by Givens transformations without square roots. J. Inst. Maths. Applics.12, 329–336 (1973).
Golub, G. H., van Loan, C.: Matrix Computations. Baltimore, Maryland: Johns Hopkins University Press 1983.
Hammerling, S.: A note on modifications to the Givens plane rotation. J. Inst. Maths. Applics13, 215–218 (1974).
Kaufman, L.: Some thoughts on the QZ algorithm for solving the generalized eigenvalue problem. ACM Trans. Math. Soft., Vol. 3, No. 1, 65–75 (1977).
Märchy, H. P.: Eine Modifikation des QZ-Algorithmus mit schnellen Givens-Rotationen Zürcher Dissertation (1986).
Moler, C. B., Stewart, G. W.: An algorithm for generalized matrix eigenvalue problems. SIAM J. Num. Anal.10, 241–256 (1973).
Rath, W.: Orthogonale Ähnlichkeitstransformationen mit schnellen Givensrotationen und ihre Anwendung für das Jacobiverfahren und den QR-Algorithmus. Diplomarbeit, Bielefeld (1980).
Rath, W.: Fast Givens rotations for orthogonal similarity transformations. Numer. Math.40, 47–56 (1982).
Schwarz, H. R.: Numerische Mathematik. Stuttgart: B. G. Teubner 1986.
Ward, R. C.: An extension of the QZ algorithm for solving the generalized matrix eigenvalue problem. NASA Tech. Note, NASA TN D-7305 (1973).
Ward, R. C.: The combination shift QZ algorithm. SIAM J. Num. Anal.12, 835–853 (1975).
Wilkinson, J. H.: The Algebraic Eigenvalue Problem. Oxford: Clarendon Press 1965.
Wilkinson, J. H.: Some recent advances in numerical linear algebra in the state of the art in numerical analysis (Jacobs, D. A. H., ed.), pp. 10–18. New York: Academic Press 1977.
IMSL: International Mathematical and Statistical Libraries Inc. (Edition 9) (1982).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Märchy, H.P. On a modification of the QZ algorithm with fast givens rotations. Computing 38, 247–259 (1987). https://doi.org/10.1007/BF02240099
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02240099