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

fix rngs issue; block Lanczos is working

Browse Source
This commit is contained in:
Yong-Chull Jang 2020-03-25 15:45:50 -04:00
parent 2db7e6f8ab
commit 9266b89ad8
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Grid/algorithms/iterative/ImplicitlyRestartedBlockLanczos.h
View File

@ -180,7 +180,7 @@ public:
Glog<<"orthogonalize after: "<<j<<" of "<<k<<" "<< ip <<std::endl; Glog<<"orthogonalize after: "<<j<<" of "<<k<<" "<< ip <<std::endl;
} }
} }
normalize(w,if_print); assert(normalize(w,if_print) != 0);
} }
void reorthogonalize(Field& w, std::vector<Field>& evec, int k) void reorthogonalize(Field& w, std::vector<Field>& evec, int k)
{ {

2
tests/lanczos/Test_dwf_block_lanczos.cc
View File

@ -363,7 +363,7 @@ int main (int argc, char ** argv)
std::vector<RealD> eval(JP.Nm); std::vector<RealD> eval(JP.Nm);
std::vector<FermionField> src(JP.Nu,FrbGrid); std::vector<FermionField> src(JP.Nu,FrbGrid);
if (1) if (0)
{ {
std::cout << GridLogMessage << "Using RNG5"<<std::endl; std::cout << GridLogMessage << "Using RNG5"<<std::endl;
FermionField src_tmp(FGrid); FermionField src_tmp(FGrid);