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

Fixed bug in LapEvec if machine running spread-out in time

Browse Source
This commit is contained in:
Michael Marshall 2019-03-26 13:24:39 +00:00
parent ae565b006a
commit 9fce1263be
2 changed files with 9 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Hadrons/Modules/MDistil/Distil.hpp
View File

@ -250,6 +250,12 @@ inline GridCartesian * MakeLowerDimGrid( GridCartesian * gridHD )
(Disable big-endian by setting Endian_Scalar_Size=1)
IndexNames contains one name for each index, and IndexNames are validated on load.
(NB: Indices of dimension 1 are not saved, and not validated on load)
WHAT TO SAVE / VALIDATE ON LOAD (Override to warn instead of assert on load)
Ensemble string
Configuration number
Noise unique string
Distillation parameters
******************************************************************************/
template<typename Scalar_, int NumIndices_, uint16_t Endian_Scalar_Size = sizeof(Scalar_)>

4
Hadrons/Modules/MDistil/LapEvec.hpp
View File

@ -328,7 +328,9 @@ void TLapEvec<GImpl>::execute(void)
for (int i=0;i<LPar.Nvec;i++){
std::cout << "Inserting Timeslice " << t << " into vector " << i << std::endl;
InsertSliceLocal(eig[t].evec[i],eig4d.evec[i],0,t,3);
eig4d.eval[i] = eig[0].eval[i];
// TODO: Discuss: is this needed? Is there a better way?
if(t==0)
eig4d.eval[i] = eig[t].eval[i];
}
}