nekStab / LightKrylov

Lightweight implementation of Krylov subspace techniques in Fortran.
BSD 3-Clause "New" or "Revised" License
18 stars 0 forks source link

Two sided arnoldi #37

Closed loiseaujc closed 8 months ago

loiseaujc commented 8 months ago

A relatively simple implementation of the two-sided Arnoldi factorization is implemented. It works out of the box although it could benefit eventually from some upgrades, most notably some restarting procedure.

loiseaujc commented 8 months ago

Note also that I moved all the lapack wrappers (i.e. evd, hevd, svd, lstsq) from IterativeSolvers.f90 to Utils.f90 for re-usability.