portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-26 05:35:55 +01:00
Code Releases Activity

relocate the printing block of converged eigenvalues

Browse Source
This commit is contained in:
Yong-Chull Jang 2018-01-31 12:40:33 -05:00
parent fbe1209f7e
commit bc5ba39278
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/algorithms/iterative/ImplicitlyRestartedBlockLanczos.h
View File

@ -286,10 +286,6 @@ until convergence
clog <<" #modes converged: "<<Nconv<<std::endl; clog <<" #modes converged: "<<Nconv<<std::endl;
if( Nconv>=Nstop ){
goto converged;
}
for(int i=0; i<Nconv; ++i){ for(int i=0; i<Nconv; ++i){
std::cout.precision(13); std::cout.precision(13);
clog << "[" << std::setw(3)<< std::setiosflags(std::ios_base::right) <<Iconv[i]<<"] "; clog << "[" << std::setw(3)<< std::setiosflags(std::ios_base::right) <<Iconv[i]<<"] ";
@ -297,6 +293,10 @@ until convergence
std::cout << " |H B[i] - eval_conv[i]B[i]|^2 "<< std::setw(25)<< std::setiosflags(std::ios_base::right)<< resid[Iconv[i]]<< std::endl; std::cout << " |H B[i] - eval_conv[i]B[i]|^2 "<< std::setw(25)<< std::setiosflags(std::ios_base::right)<< resid[Iconv[i]]<< std::endl;
} }
if( Nconv>=Nstop ){
goto converged;
}
if ( iter < MaxIter-1 ) { if ( iter < MaxIter-1 ) {
if ( Nu == 1 ) { if ( Nu == 1 ) {
// reconstruct initial vector for additional pole space // reconstruct initial vector for additional pole space