Abstract
We present a practical modification of the recent divide-and-conquer algorithms of [3] for approximating the eigenvalues of a real symmetric tridiagonal matrix. In this modified version, we avoid the numerical stability problems of the algorithms of [3] but preserve their insensitivity to clustering the eigenvalues and the possibility to give a priori upper bounds on their computational cost for any input matrix. We confirm the theoretical effectiveness of our algorithms by numerical experiments.
Zusammenfassung
Wir schlagen eine Modifikation der derzeitigen Divide-and-Conquer-Algorithmen aus [3] zur näherungsweisen Berechnung der Eigenwerte einer reelen symmetrischen Matrix vor. Dabei vermeiden wir die numerischen Stabilitätsprobleme dieser Algorithmen, erhalten aber ihre Unempfindlichkeit gegen das Auftreten von Eigenwerthaufen; wir können auch für unsere Modifikation a-priori Schranken für den Berechnungsaufwand für beliebige Matrizen angeben. Die theoretische Effizienz unseres Algorithmus bestätigen wir durch numerische Experimente.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Ben-Or, M., Tiwari, P.: Simple algorithm for approximating all roots of a polynomial with real roots. J. Complexity6(4), 417–442 (1990).
Bini, D., Gemignani, L.: An iteration scheme for the divide-and-conquer eigenvalue solvers. Manuscript, Dipartimento di Matematica, Università di Pisa, Pisa, 1991.
Bini, D., Pan, V.: Parallel complexity of tridiagonal symmetric eigenvalue problem. Proc. 2nd Annual ACM-SIAM Symp. on Discrete Algorithms, pp. 384–393, 1991.
Bunch, J. R., Nielsen C. P., Sorensen, D. C.: Rank-one modificationof the symmetric eigenproblem. Numer. Math.31, 31–48 (1978).
Cuppen, J. J. M.: A divide and conquer method for the symmetric tridiagonal eigenproblem. Numer. Math.36, 177–195 (1981).
Dongarra, J. J., Sorensen, D. C.: A fully parallel algorithm for the symmetric eigenvalue problem. SIAM J. on Sci. and Stat. Computing8 (2), 139–154 (1987).
Gill, D., Tadmor, E., AnO(n2) method for computing the eigensystem ofn×n symmetric tridiagonal matrices by the divide and conquer approach. SIAM J. Sci. Stat. Comput.11, 1, 161–173 (1990).
Golub, G. H., Van Loan, C. F.: Matrix computations. Baltimore, MD: Johns Hopkins Univ. Press, 1989.
Krishnakumar, A. S., Morf, M.: Eigenvalues of a symmetric tridiagonal matrix: a divide-and-conquer approach. Numer. Math.48, 349–368 (1986).
Parlett, B. N.: The symmetric eigenvalue problem. Englewood Cliffs, NJ: Prentice-Hall, 1980.
Pereyra, V., Scherer, G.: Eigenvalues of symmetric tridiagonal matrices: a fast, accurate and reliable algorithm. J. Inst. Maths Applics.12, 209–222 (1973).
Press, W. H., Flannery, B. P., Tenkolsky, S. A., Vetterling W. T.: Numerical recipes. Cambridge: Cambridge University Press, 1986.
Renegar, J.: On the worst-case arithmetic complexity of approximating zeros of polynomials. J. Complexity3 (2), 90–113 (1987).
Author information
Authors and Affiliations
Additional information
Dedicated to Professor Willard L. Miranker on the occasion of his 60th birthday
Supported by NSF Grants CCR-8805782 and CCR-9020690 and by 40% and 60% funds of MURST.
Supported by NSF Grants CCR-8805782 and CCR 9020690 and PSC-CUNY Award 661340.
Rights and permissions
About this article
Cite this article
Bini, D., Pan, V. Practical improvement of the divide-and-conquer eigenvalue algorithms. Computing 48, 109–123 (1992). https://doi.org/10.1007/BF02241709
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02241709