portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Merge remote-tracking branch 'origin/develop' into develop

Browse Source
This commit is contained in:
david clarke 2023-08-11 23:06:42 -06:00
commit 9c22655b5a
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);
} }