Computing generalized eigenvalues

[daanhb]

This is the PR mentioned in #126 and based on code provided to me by @thijssteel.

It implements eigval(A, B) for generic matrices A and B using the QZ algorithm. Copying a comment from the code:

Two orthogonal matrices Q and Z are constructed to reduce the "pencil" (A,B) to another form Q'*(A,B)*Z:

• in a first step A and B are reduced to upper Hessenberg and upper triangular form, respectively (see: hesstriangular)
• next a generalized Schur decomposition further reduces both matrices to upper triangular form
• the generalized eigenvalues are then given by the ratios of the diagonals of those two triangular matrices.

[daanhb]

I have yet to add tests. First I'm wondering whether this package is the right place to add this functionality and whether it's deemed useful. It is at least to me :-)

[codecov]

Codecov Report

Attention: Patch coverage is 0% with 525 lines in your changes missing coverage. Please review.

Project coverage is 72.74%. Comparing base (27a78e0) to head (6a5beeb).
Report is 17 commits behind head on master.

Files with missing lines	Patch %	Lines
src/qz.jl	0.00%	525 Missing ⚠️

Additional details and impacted files

@@             Coverage Diff             @@
##           master     #127       +/-   ##
===========================================
- Coverage   96.08%   72.74%   -23.35%     
===========================================
  Files          11       12        +1     
  Lines        1636     2161      +525     
===========================================
  Hits         1572     1572               
- Misses         64      589      +525     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[daanhb]

Looking at this code, at other parts of GenericLinearAlgebra and at LinearAlgebra, I see many opportunities for better integration (and less code in total). For this PR:

• reuse LinearAlgebra.floatmin2
• better share code for Givens rotations, Householder reflections and eigenvalues of 2x2 blocks
• make computation of eigenvalues more consistent. LinearAlgebra returns real eigenvalues when they are real, GenericLinearAlgebra always returns complex vectors (it seems?) and I've followed LinearAlgebra

It might make the implementation quite a bit shorter.

[andreasnoack]

I'm sorry for the silence here. This contributions is highly appreciated. I'm a bit stretched for the time being but will try to get this reviewed over the summer.

[daanhb]

Closing this PR in favour of an improved version in the future, explanation in the issue: #126