I need a generalized
eigenvector solver (A*x = lambda*B*x) for the case where the matrix B is not
necessarily positive definate. The routine used by
vnl_generalized_eigensystem requires the matrix B be positive
definate.

Lapack version 3 has a routine
(dggev) for solving this case. This is also the Lapack routine used by
MATLAB. It looks like vnl uses Lapack version 2. Are there plans to
upgrade vnl to use version 3 of Lapack? Version 3 came out in 1999.

I tried grabbing just this file
from Clapack but as you might expect it does not compile with the rest of
version 2.

There is also an Eispack
routine for this problem rgg but vnl does not seem to have that routine
either.

**Jim Miller**
*_____________________________________**Visualization &
Computer Vision**GE
Research**Bldg. KW, Room C218B**P.O. Box 8, Schenectady NY
12301*__millerjv@research.ge.com__

__james.miller@research.ge.com__*(518) 387-4005, Dial
Comm: 8*833-4005, **Cell: (518) 505-7065,
Fax: (518) 387-6981*