Re: [eigen] generalized eigenvector solver bug

[Gael Guennebaud <gael.guennebaud@xxxxxxxxx>]

For the self-adjoint case, if one of the matrices is positive-definite, then the class GeneralizedSelfAdjointEigenSolver is already pretty good. It is based on a LL^T factorizarion of the positive-definite one, and on the SelfAdjointEigenSolver class to carry out the actual eigen-problem.

For the non-symmetric but real case, the class GeneralizedEigenSolver is based on the quite standard QZ algorithm (mostly from Golub's Matrix Computations book), but only the eigenvalues are returned and we're only missing the bits to extract the eigenvectors from QZ and eigenvalues. Improving the QZ implementation as suggested by Pavel won't change the fact that the aforementioned part of the puzzle (aka bug 645) will be missing and should probably be worked out first.

gael

On Sun, Mar 8, 2015 at 4:09 PM, Gareth Owen <gareth.w.owen@xxxxxxxxx> wrote:

Hi,

I've been looking over solvers for the generalized eigenvalue problem a bit this week, and I was wondering which implementation would be meant for the generalized eigensolver within Eigen - generalized symmetric definite or generalized nonsymmetric? Which might be best to try and implement first?

thanks

Gareth

On Fri, Mar 6, 2015 at 3:18 PM, Gael Guennebaud <gael.guennebaud@xxxxxxxxx> wrote:

As far as I known, nothing has been started yet.

gael

On Fri, Mar 6, 2015 at 12:48 PM, Gareth Owen <gareth.w.owen@xxxxxxxxx> wrote:

Hi,

I've noticed that the generalized eigensolver bug (645) is still open - is there any work being done on this?

thanks

Gareth