RE: [Vxl-users] generalized eigenvectors

[Miller, J. V (Research) <mil...@cr...>]

I was able to get the dggev routine from Lapack 3 to build with the rest of vnl.  I needed to grab
the blaswrap.h from Lapack 3 so dggev would compile. (Don't know whether it works yet :)
 
I hesitate to check these in since we'd be mixing Lapack 3 files with Lapack 2 files which could be
confusing to the users and maintainers.
 
Jim

-----Original Message-----
From: Miller, James V (Research) 
Sent: Thursday, June 20, 2002 2:01 PM
To: Vxl-users (E-mail)
Subject: [Vxl-users] generalized eigenvectors


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

mil...@re... <mailto:mil...@re...> 

jam...@re...
(518) 387-4005, Dial Comm: 8*833-4005, 
Cell: (518) 505-7065, Fax: (518) 387-6981