Hi Jessica,
Resuming iterations is already possible now. The Lanczos algorithm allows you to continue iterating with a new iteration object. In your case I would not iterate until convergence but only for a fixed number of steps.
Indeed, all you have to store are the member variables of the Lanczos (and underlying T-matrix) object, as well as the two last vectors of the iteration (not the starting vector!). It is not hard to add that.
Matthias
On 22 Jul 2011, at 04:22, jessica.alfonsi@unipd.it wrote:
Hi all, I'd like to modify the Lanczos *.cpp examples in IETL example folder to be able to restart the iteration from any interruption point during the Lanczos iteration. I suppose first one has to dump some information, for instance the cofficients of T-matrix and the starting vector, however the more difficult part would be the resuming of the iteration. Can somebody give me some hints on how accomplishing this task? Thanks in advance...
Best regards,
Jessica Alfonsi