portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 19:25:56 +01:00
Code Releases Activity

do one more iteration if not all vectors converged

Browse Source
This commit is contained in:
Felix Erben 2023-06-27 11:33:30 +01:00
parent ee92e08edb
commit 014704856f
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h
View File

@ -284,6 +284,7 @@ until convergence
int iter; int iter;
for(iter = 0; iter<MaxIter; ++iter){ for(iter = 0; iter<MaxIter; ++iter){
notYetConverged:
OrthoTime=0.; OrthoTime=0.;
std::cout<< GridLogMessage <<" **********************"<< std::endl; 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 << "Nconv ("<<Nconv<<") < Nstop ("<<Nstop<<")"<<std::endl;
std::cout << GridLogIRL << "starting one more iteration"<<std::endl;
iter++;
goto notYetConverged;
}
eval=eval2; eval=eval2;
//Keep only converged //Keep only converged