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

Bug fix in the coarse restore...

Browse Source 
Think this is nearly there
This commit is contained in:
paboyle 2017-10-27 10:29:34 +01:00
parent f96c800d25
commit aa66f41c69
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
tests/lanczos/Test_dwf_compressed_lanczos_reorg.cc
View File

@ -109,7 +109,7 @@ public:
void checkpointCoarseRestore(std::string evecs_file,std::string evals_file,int nvec) void checkpointCoarseRestore(std::string evecs_file,std::string evals_file,int nvec)
{ {
std::cout << " resizing to " << nvec<< std::endl; std::cout << "resizing coarse vecs to " << nvec<< std::endl;
this->evals_coarse.resize(nvec); this->evals_coarse.resize(nvec);
this->evec_coarse.resize(nvec,this->_CoarseGrid); this->evec_coarse.resize(nvec,this->_CoarseGrid);
std::cout << GridLogIRL<< "checkpointCoarseRestore: Reading evals from "<<evals_file<<std::endl; std::cout << GridLogIRL<< "checkpointCoarseRestore: Reading evals from "<<evals_file<<std::endl;
@ -121,7 +121,7 @@ public:
std::cout << GridLogIRL<< "checkpointCoarseRestore: Reading evecs from "<<evecs_file<<std::endl; std::cout << GridLogIRL<< "checkpointCoarseRestore: Reading evecs from "<<evecs_file<<std::endl;
Grid::QCD::ScidacReader RD ; Grid::QCD::ScidacReader RD ;
RD.open(evecs_file); RD.open(evecs_file);
for(int k=0;k<nbasis;k++) { for(int k=0;k<nvec;k++) {
RD.readScidacFieldRecord(this->evec_coarse[k],record); RD.readScidacFieldRecord(this->evec_coarse[k],record);
} }
RD.close(); RD.close();