nekStab / LightKrylov

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

Add default `select` procedure for `eigs` #90

Closed loiseaujc closed 4 months ago

loiseaujc commented 4 months ago

At the moment, eigs returns tentative eigenvectors and eigenvalues even if not converged. This is a design problem.

By default, there are not selection procedure and thus no krylov-schur restart. The basis of eigenvectors is pre-allocated. If kdim is not passed, it performs 4*nev iterations before stopping. If kdim is passed, it performs kdim iterations.

Two options:

loiseaujc commented 4 months ago

Closed with #95