Computing Eigenelements of Real Symmetric Matrices via Optimization

  • Authors:
  • M. Mongeau;M. Torki

  • Affiliations:
  • Labo MIP, Université Paul Sabatier, 31062 Toulouse Cedex 04, France. mongeau@mip.ups-tlse.fr;Labo MIP, Université Paul Sabatier, 31062 Toulouse Cedex 04, France. torki@iup.univ-avignon.fr

  • Venue:
  • Computational Optimization and Applications
  • Year:
  • 2004

Quantified Score

Hi-index 0.00

Visualization

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 and 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).