F08SEF (DSYGST) reduces a real symmetric-definite generalized eigenproblem
,
or
to the standard form
, where
is a real symmetric matrix and
has been factorized by
F07FDF (DPOTRF).
To reduce the real symmetric-definite generalized eigenproblem
,
or
to the standard form
, F08SEF (DSYGST) must be preceded by a call to
F07FDF (DPOTRF) which computes the Cholesky factorization of
;
must be positive definite.
The different problem types are specified by the parameter
ITYPE, as indicated in the table below. The table shows how
is computed by the routine, and also how the eigenvectors
of the original problem can be recovered from the eigenvectors of the standard form.
ITYPE |
Problem |
UPLO |
|
|
|
|
|
'U'
'L' |
|
|
|
|
|
'U'
'L' |
|
|
|
|
|
'U'
'L' |
|
|
|
Forming the reduced matrix
is a stable procedure. However it involves implicit multiplication by
(if
) or
(if
or
). When F08SEF (DSYGST) is used as a step in the computation of eigenvalues and eigenvectors of the original problem, there may be a significant loss of accuracy if
is ill-conditioned with respect to inversion.
See the document for
F08SAF (DSYGV) for further details.
The complex analogue of this routine is
F08SSF (ZHEGST).
This example computes all the eigenvalues of
, where
Here
is symmetric positive definite and must first be factorized by
F07FDF (DPOTRF). The program calls F08SEF (DSYGST) to reduce the problem to the standard form
; then
F08FEF (DSYTRD) to reduce
to tridiagonal form, and
F08JFF (DSTERF) to compute the eigenvalues.