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

test works up to perambulators now

Browse Source
This commit is contained in:
ferben 2019-01-25 13:44:19 +00:00
parent 7496da0987
commit 9ef6f9878e
2 changed files with 3 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -279,17 +279,12 @@ 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;
ExtractSliceLocal(evec3d,epack.evec[ivec],0,t,3); ExtractSliceLocal(evec3d,epack.evec[ivec],0,t,3);
std::cout << "2" << std::endl;
pokeSpin(perambulator(t, ivec, dk, inoise,dt,ds),innerProduct(evec3d, result_3d),is); pokeSpin(perambulator(t, ivec, dk, inoise,dt,ds),innerProduct(evec3d, result_3d),is);
std::cout << "3" << std::endl;
} }
} }
} }

9
tests/hadrons/Test_hadrons_distil.cc
View File

@ -89,10 +89,10 @@ void test_Perambulators(Application &application)
PerambPar.eigenPack="LapEvec"; PerambPar.eigenPack="LapEvec";
PerambPar.tsrc = 0; PerambPar.tsrc = 0;
PerambPar.nnoise = 1; PerambPar.nnoise = 1;
PerambPar.LI=6; PerambPar.LI=5;
PerambPar.SI=4; PerambPar.SI=4;
PerambPar.TI=64; PerambPar.TI=64;
PerambPar.nvec=6; PerambPar.nvec=5;
PerambPar.Ns=4; PerambPar.Ns=4;
PerambPar.Nt=64; PerambPar.Nt=64;
PerambPar.Nt_inv=1; PerambPar.Nt_inv=1;
@ -113,7 +113,7 @@ void test_DistilVectors(Application &application)
MDistil::DistilVectors::Par DistilVecPar; MDistil::DistilVectors::Par DistilVecPar;
DistilVecPar.noise="Peramb_noise"; DistilVecPar.noise="Peramb_noise";
DistilVecPar.perambulator="Peramb_perambulator_light"; DistilVecPar.perambulator="Peramb_perambulator_light";
DistilVecPar.eigenPack="LapEvec_eigenPack"; DistilVecPar.eigenPack="LapEvec";
DistilVecPar.tsrc = 0; DistilVecPar.tsrc = 0;
DistilVecPar.nnoise = 1; DistilVecPar.nnoise = 1;
DistilVecPar.LI=6; DistilVecPar.LI=6;
@ -123,9 +123,6 @@ void test_DistilVectors(Application &application)
DistilVecPar.Ns=4; DistilVecPar.Ns=4;
DistilVecPar.Nt=64; DistilVecPar.Nt=64;
DistilVecPar.Nt_inv=1; DistilVecPar.Nt_inv=1;
// gauge field
//application.createModule<MGauge::Unit>("gauge");
// Now make an instance of the DistilVectors object
application.createModule<MDistil::DistilVectors>("DistilVectorsInstance",DistilVecPar); application.createModule<MDistil::DistilVectors>("DistilVectorsInstance",DistilVecPar);
} }