From 78bae9417c7da1df6677582227894208064edafb Mon Sep 17 00:00:00 2001
From: Felix Erben <felix.erben@ed.ac.uk>
Date: Tue, 27 Jun 2023 14:38:19 +0100
Subject: [PATCH] returning Nstop vectors even if not all meet true convergence
 criterion

---
 Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h | 9 +++------
 1 file changed, 3 insertions(+), 6 deletions(-)

diff --git a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
index d2f95734..d39701f4 100644
--- a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
+++ b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
@@ -284,7 +284,6 @@ until convergence
     int iter;
     for(iter = 0; iter<MaxIter; ++iter){
       
-  notYetConverged:
       OrthoTime=0.;
 
       std::cout<< GridLogMessage <<" **********************"<< std::endl;
@@ -422,15 +421,13 @@ until convergence
 
       if ( Nconv < Nstop ) {
 	std::cout << GridLogIRL << "Nconv ("<<Nconv<<") < Nstop ("<<Nstop<<")"<<std::endl;
-	std::cout << GridLogIRL << "starting one more iteration"<<std::endl;
-	iter++;
- 	goto notYetConverged;
+	std::cout << GridLogIRL << "returning Nstop vectors, the last "<< Nstop-Nconv << "of which might meet convergence criterion only approximately" <<std::endl;
       }
       eval=eval2;
       
       //Keep only converged
-      eval.resize(Nconv);// Nstop?
-      evec.resize(Nconv,grid);// Nstop?
+      eval.resize(Nstop);// was Nconv
+      evec.resize(Nstop,grid);// was Nconv
       basisSortInPlace(evec,eval,reverse);
       
     }