portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

bugfix in prambLight

Browse Source
This commit is contained in:
ferben 2019-01-25 13:08:56 +00:00
parent 2568f5b925
commit 7496da0987
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Hadrons/Modules/MDistil/PerambLight.hpp
View File

@ -102,10 +102,11 @@ void TPerambLight<FImpl>::setup(void)
int TI=par().TI; int TI=par().TI;
int nnoise=par().nnoise; int nnoise=par().nnoise;
int Nt=par().Nt; int Nt=par().Nt;
int Nt_inv=par().Nt_inv;
int nvec=par().nvec; int nvec=par().nvec;
envCreate(Perambulator<SpinVector>, getName() + "_perambulator_light", 1, envCreate(Perambulator<SpinVector>, getName() + "_perambulator_light", 1,
LI*SI*TI*nnoise*nvec*Nt); Nt,nvec,LI,nnoise,Nt_inv,SI);
envCreate(std::vector<Complex>, getName() + "_noise", 1, envCreate(std::vector<Complex>, getName() + "_noise", 1,
nvec*Ns*Nt*nnoise); nvec*Ns*Nt*nnoise);
@ -278,8 +279,10 @@ void TPerambLight<FImpl>::execute(void)
// current_sink[inoise+nnoise*(dk+LI*(dt+Nt_inv*ds))] = result; // current_sink[inoise+nnoise*(dk+LI*(dt+Nt_inv*ds))] = result;
std::cout << "Contraction of perambulator from noise " << inoise << " and dilution component (d_k,d_t,d_alpha) : (" << dk << ","<< dt << "," << ds << ")" << std::endl; std::cout << "Contraction of perambulator from noise " << inoise << " and dilution component (d_k,d_t,d_alpha) : (" << dk << ","<< dt << "," << ds << ")" << std::endl;
for (int is = 0; is < Ns; is++) { for (int is = 0; is < Ns; is++) {
std::cout << "is" << is << std::endl;
result_nospin = peekSpin(result,is); result_nospin = peekSpin(result,is);
for (int t = Ntfirst; t < Ntfirst + Ntlocal; t++) { for (int t = Ntfirst; t < Ntfirst + Ntlocal; t++) {
std::cout << "t" << t << std::endl;
ExtractSliceLocal(result_3d,result_nospin,0,t-Ntfirst,Grid::QCD::Tdir); ExtractSliceLocal(result_3d,result_nospin,0,t-Ntfirst,Grid::QCD::Tdir);
for (int ivec = 0; ivec < nvec; ivec++) { for (int ivec = 0; ivec < nvec; ivec++) {
std::cout << "1" << std::endl; std::cout << "1" << std::endl;