nekStab / LightKrylov

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

Krylov-Schur restarting procedure #58

Closed loiseaujc closed 7 months ago

loiseaujc commented 7 months ago

This PR adds two main utilities to LightKrylov :

Note that, although the code compiles and runs for restarting the block Arnoldi, it did not always passed the test. It seems to be related to how many eigenvalues have been selected relative to the blocksize. For the moment, if a block factorization is used, the routine will raise an automatic error.