From 014704856fc0ca5025989f852d00bb434d7731b9 Mon Sep 17 00:00:00 2001
From: Felix Erben <felix.erben@ed.ac.uk>
Date: Tue, 27 Jun 2023 11:33:30 +0100
Subject: [PATCH 1/3] do one more iteration if not all vectors converged

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

diff --git a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
index 3d0a2a75..12245a2d 100644
--- a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
+++ b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
@@ -284,6 +284,7 @@ until convergence
     int iter;
     for(iter = 0; iter<MaxIter; ++iter){
       
+  notYetConverged:
       OrthoTime=0.;
 
       std::cout<< GridLogMessage <<" **********************"<< std::endl;
@@ -419,9 +420,12 @@ until convergence
 	}
       }
 
-      if ( Nconv < Nstop )
+      if ( Nconv < Nstop ) {
 	std::cout << GridLogIRL << "Nconv ("<<Nconv<<") < Nstop ("<<Nstop<<")"<<std::endl;
-
+	std::cout << GridLogIRL << "starting one more iteration"<<std::endl;
+    iter++;
+    goto notYetConverged;
+      }
       eval=eval2;
       
       //Keep only converged

From dd170ead0170550297a9cec8ce26518ded14d1b6 Mon Sep 17 00:00:00 2001
From: Felix Erben <felix.erben@ed.ac.uk>
Date: Tue, 27 Jun 2023 11:37:01 +0100
Subject: [PATCH 2/3] whitespace

---
 Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
index 12245a2d..d2f95734 100644
--- a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
+++ b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
@@ -423,8 +423,8 @@ 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;
+	iter++;
+ 	goto notYetConverged;
       }
       eval=eval2;
       

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 3/3] 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);
       
     }