portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-25 13:15:55 +01:00
Code Releases Activity

added some TODO comments; needs discussion

Browse Source
This commit is contained in:
Felix Erben 2019-01-31 16:54:11 +00:00
parent b6b267fd4b
commit 7cc13f48d5
2 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Hadrons/Modules/MDistil/DistilVectors.hpp
View File

@ -156,7 +156,7 @@ void TDistilVectors<FImpl>::execute(void)
int nnoise=par().nnoise; int nnoise=par().nnoise;
int LI=par().LI; int LI=par().LI;
int Ns=par().Ns; int Ns=par().Ns;
int Nt_inv=par().Nt_inv; int Nt_inv=par().Nt_inv; // TODO: No input, but define through Nt, TI
int Nt=par().Nt; int Nt=par().Nt;
int TI=par().TI; int TI=par().TI;
int nvec=par().nvec; int nvec=par().nvec;

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

@ -105,7 +105,6 @@ std::vector<std::string> TPerambLight<GImpl>::getInput(void)
{ {
std::vector<std::string> in; std::vector<std::string> in;
//in.push_back(par().noise);
in.push_back(par().eigenPack); in.push_back(par().eigenPack);
return in; return in;
@ -175,7 +174,7 @@ void TPerambLight<GImpl>::execute(void)
const int TI{Distil.TI}; const int TI{Distil.TI};
const int nnoise{Distil.nnoise}; const int nnoise{Distil.nnoise};
const int Nt{Distil.Nt}; const int Nt{Distil.Nt};
const int Nt_inv{Distil.Nt_inv}; const int Nt_inv{Distil.Nt_inv}; // TODO: PROBABLY BETTER: if (full_tdil) Nt_inv=1; else Nt_inv = TI;
const int tsrc{Distil.tsrc}; const int tsrc{Distil.tsrc};
const int Ns{Distil.Ns}; const int Ns{Distil.Ns};
@ -242,7 +241,7 @@ void TPerambLight<GImpl>::execute(void)
//Create Noises //Create Noises
//std::cout << pszGaugeConfigFile << std::endl; //std::cout << pszGaugeConfigFile << std::endl;
//GridSerialRNG sRNG; sRNG.SeedUniqueString(std::string(pszGaugeConfigFile)); //GridSerialRNG sRNG; sRNG.SeedUniqueString(std::string(pszGaugeConfigFile));
GridSerialRNG sRNG; sRNG.SeedUniqueString("unique_string"); GridSerialRNG sRNG; sRNG.SeedUniqueString("unique_string"); // TODO: Proper unique string. Include quark mass, gauge field? Maybe also nvec, but in a way that more nvec would only add noises, not change all of them???
Real rn; Real rn;
for (int inoise=0;inoise<nnoise;inoise++) { for (int inoise=0;inoise<nnoise;inoise++) {
@ -296,7 +295,7 @@ void TPerambLight<GImpl>::execute(void)
if (full_tdil) t_inv = tsrc; else t_inv = it; if (full_tdil) t_inv = tsrc; else t_inv = it;
if( t_inv >= Ntfirst && t_inv < Ntfirst + Ntlocal ) { if( t_inv >= Ntfirst && t_inv < Ntfirst + Ntlocal ) {
for (int ik = dk; ik < nvec; ik += LI){ for (int ik = dk; ik < nvec; ik += LI){
for (int is = ds; is < Ns; is += Ns){ //at the moment, full spin dilution is enforced for (int is = ds; is < Ns; is += Ns){ // TODO: Also allow non-full spin dilution (re-define exact_distillation?)
ExtractSliceLocal(evec3d,epack.evec[ik],0,t_inv,3); ExtractSliceLocal(evec3d,epack.evec[ik],0,t_inv,3);
tmp3d_nospin = evec3d * noise[inoise + nnoise*(t_inv + Nt*(ik+nvec*is))]; tmp3d_nospin = evec3d * noise[inoise + nnoise*(t_inv + Nt*(ik+nvec*is))];
tmp3d=zero; tmp3d=zero;
@ -315,6 +314,7 @@ void TPerambLight<GImpl>::execute(void)
Dop.ImportPhysicalFermionSource(dist_source,src5); Dop.ImportPhysicalFermionSource(dist_source,src5);
SchurSolver(Dop,src5,sol5); SchurSolver(Dop,src5,sol5);
Dop.ExportPhysicalFermionSolution(sol5,result); //These are the meson sinks Dop.ExportPhysicalFermionSolution(sol5,result); //These are the meson sinks
// TODO: Can we inherit something from MContraction to compute the fourier-transformed sinks???
//if (compute_current_sink) //if (compute_current_sink)
// 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;