Abstract
In certain circumstances, it is advantageous to use an optimization approach in order to solve the generalized eigenproblem, Ax = λBx, where A and B are real symmetric matrices and B is positive definite. In particular, this is the case when the matrices A and B are very large the computational cost, prohibitive, of solving, with high accuracy, systems of equations involving these matrices. Usually, the optimization approach involves optimizing the Rayleigh quotient.
We first propose alternative objective functions to solve the (generalized) eigenproblem via (unconstrained) optimization, and we describe the variational properties of these functions.
We then introduce some optimization algorithms (based on one of these formulations) designed to compute the largest eigenpair. According to preliminary numerical experiments, this work could lead the way to practical methods for computing the largest eigenpair of a (very) large symmetric matrix (pair).
Similar content being viewed by others
References
G. Auchmuty, “Unconstrained variational principles for eigenvalues of real symmetric matrices,” SIAM Journal on Mathematical Analysis, vol. 20, no. 5, pp. 1186–1207, 1989.
G. Auchmuty, “Globally and rapidly convergent algorithms for symmetric eigenproblems,” SIAM Journal on Matrix Analysis and Applications, vol. 12, pp. 690–706, 1991.
Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide. SIAM, 2000.
J.F. Bonnans, J.-C. Gilbert, C. Lemaréchal, and C. Sagastizábal, Optimisation Numérique, Springer, 1997.
R.H. Byrd, J. Nocedal, and R.B. Schnabel, “Representations of quasi-Newton matrices and their use in limited-memory,” Mathematical Programming, vol. 63, pp. 129–156, 1994.
F.H. Clarke, Optimization and Nonsmooth Analysis, John Wiley, 1983.
R.S. Dembo and T. Steihaug, “Truncated-Newton algorithms for large-scale unconstrained optimization,” Mathematical Programming, vol. 26, pp. 190–212, 1983.
J.W. Demmel, Applied Numerical Linear Algebra, SIAM, 1997.
J.E. Dennis and R.B. Schnabel, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, Prentice-Hall, 1983.
A. Edelman and S.T. Smith, “On conjugate gradient-like methods for eigen-like problems,” BIT, vol. 36, no. 3, pp. 494–508, 1996.
D.C. Liu and J. Nocedal, “On the limited memory method for large scale optimization,” Mathematical Programming, vol. 45, no. 3, pp. 503–528, 1989.
D.E. Longsine and S.F. McCormick, “Simultaneous Rayleigh-quotient minimization methods for Ax = λBx” Linear Algebra and its Applications, vol. 34, pp. 195–234, 1980.
Matrix Market, A visual repository of test data for use in comparative studies of algorithms for numerical linear algebra. See web site http://math.nist.gov/MatrixMarket/
M. Mongeau and M. Torki, “Computing eigenelements of real symmetric matrices via optimization,” Technical Report MIP 99–54, Université Paul Sabatier, Toulouse, France, 1999. Available at web site http://mip.ups-tlse.fr.
J.L. Nazareth, H. Wolkowicz, and M. Zhu, “The quasi-Cauchy relation and diagonal updating,” SIAM Journal on Optimization, vol. 9, no. 4, pp. 1192–1204, 1999.
A.H. Sameh and J.A. Wisniewski, “A trace minimization algorithm for the generalized eigenvalue problem,” SIAM Journal on Numerical Analysis, vol. 19, no. 6, pp. 1243–1259, 1982.
D.C. Sorensen, “Numerical methods for large eigenvalue problems,” Acta Numerica, vol. 11, pp. 519–584, 2002.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Mongeau, M., Torki, M. Computing Eigenelements of Real Symmetric Matrices via Optimization. Computational Optimization and Applications 29, 263–287 (2004). https://doi.org/10.1023/B:COAP.0000044182.33308.82
Issue Date:
DOI: https://doi.org/10.1023/B:COAP.0000044182.33308.82