Generalized symmetric-definite eigenvalue problems are as follows: find the eigenvalues

λ

and the corresponding eigenvectors

z

that satisfy one of these equations:

Az

=

λ

Bz

,

ABz

=

λ

z

, or

B

A

z

=

λ

z

,

where

A

is an

n

-by-

n

symmetric or Hermitian matrix, and

B

is an

n

-by-

n

symmetric positive-definite or Hermitian positive-definite matrix.

In these problems, there exist

n

real eigenvectors corresponding to real eigenvalues (even for complex Hermitian matrices

A

and

B

).

Routines described in this

topic

allow you to reduce the above generalized problems to standard symmetric eigenvalue problem

Cy

=

λ

y

, which you can solve by calling LAPACK routines described earlier in this chapter (see Symmetric Eigenvalue Problems).

Different routines allow the matrices to be stored either conventionally or in packed storage. Prior to reduction, the positive-definite matrix

B

must first be factorized using either potrf or pptrf.

The reduction routine for the banded matrices

A

and

B

uses a split Cholesky factorization for which a specific routine pbstf is provided. This refinement halves the amount of work required to form matrix

C

.

Table

"Computational Routines for Reducing Generalized Eigenproblems to Standard Problems"

lists LAPACK routines that can be used to solve generalized symmetric-definite eigenvalue problems.