portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Merge branch 'develop' into hisq_fat_links

Browse Source
This commit is contained in:
david clarke 2023-08-11 23:08:04 -06:00
commit affff3865f
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -419,14 +419,15 @@ 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 << "returning Nstop vectors, the last "<< Nstop-Nconv << "of which might meet convergence criterion only approximately" <<std::endl;
}
eval=eval2; eval=eval2;
//Keep only converged //Keep only converged
eval.resize(Nconv);// Nstop? eval.resize(Nstop);// was Nconv
evec.resize(Nconv,grid);// Nstop? evec.resize(Nstop,grid);// was Nconv
basisSortInPlace(evec,eval,reverse); basisSortInPlace(evec,eval,reverse);
} }