Merge pull request #439 from felixerben/bugfix/IRL_convergence

Bugfix/irl convergence
returning Nstop vectors even if not all meet true convergence criterion
2024-11-09 23:45:36 +00:00 · 2023-07-12 16:32:26 -04:00 · 2023-06-27 14:38:19 +01:00 · 2023-06-27 11:37:01 +01:00 · 2023-06-27 11:33:30 +01:00
1 changed files with 5 additions and 4 deletions
--- a/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
+++ b/Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
@ -419,14 +419,15 @@ until convergence
 	}
      }

-      if ( Nconv < Nstop )
+      if ( Nconv < Nstop ) {
 	std::cout << GridLogIRL << "Nconv ("<<Nconv<<") < Nstop ("<<Nstop<<")"<<std::endl;
-
+	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);
      
    }
Author	SHA1	Message	Date
Peter Boyle	994512048e	Merge pull request #439 from felixerben/bugfix/IRL_convergence Bugfix/irl convergence	2023-07-12 16:32:26 -04:00
Felix Erben	78bae9417c	returning Nstop vectors even if not all meet true convergence criterion	2023-06-27 14:38:19 +01:00
Felix Erben	dd170ead01	whitespace	2023-06-27 11:37:01 +01:00
Felix Erben	014704856f	do one more iteration if not all vectors converged	2023-06-27 11:33:30 +01:00