sparse-linear-algebra
sparse-linear-algebra copied to clipboard
Eigensolvers should work over Complex
Lift real values onto the complex plane toC r = r :+ 0
NB: Givens' rotation in the complex plane: Saad sec. 6.5.9
Current implementation of (real-valued) Givens' rotation introduces a Ord constraint due to the comparison between the abs of two values. (Whereas magnitude a < magnitude b does not).
Depends on QR factorization working over C (#27) , since both QR method (duh) and Rayleigh iteration (which uses <\> which in turn depends on GMRES) do.