nk027/sanic

Consider supporting eigenvalue problems

nk027 opened this issue · 2 comments

nk027 commented

Also iteratively, Arnoldi and Lanzcos should be implementable.

nk027 commented

Arnoldi and Lanczos (+ tridiagonalisation and Hessenberg solvers & Eigen methods for transforming) seem to be working. Atm Lanczos can be run without (which means we find ghost eigenvalues) or full orthogonalisation. It would be interesting to implement a selective procedure (see https://books.google.at/books?id=mlOa7wPX6OYC&lpg=PA470&hl=de&pg=PA483 or Parlett and Scott 1979).

nk027 commented

As of bbf4930, this is implemented (partly based on the Spectra library, which is more specialised for this task).