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 PMTo:Vxl-users (E-mail)Subject:[Vxl-users] generalized eigenvectorsI 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 VisionGE ResearchBldg. KW, Room C218BP.O. Box 8, Schenectady NY 12301millerjv@research.ge.com

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