portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-09 23:45:36 +00:00
Code Releases Activity

Ta project

Browse Source
This commit is contained in:
Peter Boyle 2023-06-13 11:56:11 -04:00
parent 723eadbb5c
commit e3e1cc1962
1 changed files with 18 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
Grid/qcd/smearing/GaugeConfigurationMasked.h
View File

@ -696,10 +696,10 @@ private:
public: public:
/* Standard constructor */ /* Standard constructor */
SmearedConfigurationMasked(GridCartesian* _UGrid, unsigned int Nsmear, Smear_Stout<Gimpl>& Stout,bool domask=false) SmearedConfigurationMasked(GridCartesian* _UGrid, unsigned int Nsmear, Smear_Stout<Gimpl>& Stout)
: SmearedConfiguration<Gimpl>(_UGrid, Nsmear,Stout) : SmearedConfiguration<Gimpl>(_UGrid, Nsmear,Stout)
{ {
if(domask) assert(Nsmear%(2*Nd)==0); // Or multiply by 8?? assert(Nsmear%(2*Nd)==0); // Or multiply by 8??
// was resized in base class // was resized in base class
assert(this->SmearedSet.size()==Nsmear); assert(this->SmearedSet.size()==Nsmear);
@ -712,26 +712,20 @@ public:
for (unsigned int i = 0; i < this->smearingLevels; ++i) { for (unsigned int i = 0; i < this->smearingLevels; ++i) {
masks.push_back(*(new LatticeLorentzComplex(_UGrid))); masks.push_back(*(new LatticeLorentzComplex(_UGrid)));
if (domask) {
int mu= (i/2) %Nd; int mu= (i/2) %Nd;
int cb= (i%2); int cb= (i%2);
LatticeComplex tmpcb(UrbGrid); LatticeComplex tmpcb(UrbGrid);
masks[i]=Zero(); masks[i]=Zero();
//////////////////// ////////////////////
// Setup the mask // Setup the mask
//////////////////// ////////////////////
tmp = Zero(); tmp = Zero();
pickCheckerboard(cb,tmpcb,one); pickCheckerboard(cb,tmpcb,one);
setCheckerboard(tmp,tmpcb); setCheckerboard(tmp,tmpcb);
PokeIndex<LorentzIndex>(masks[i],tmp, mu); PokeIndex<LorentzIndex>(masks[i],tmp, mu);
} else {
for(int mu=0;mu<Nd;mu++){
PokeIndex<LorentzIndex>(masks[i],one, mu);
}
}
} }
delete UrbGrid; delete UrbGrid;
} }
@ -764,10 +758,14 @@ public:
tmp_mu = peekLorentz(*this->ThinLinks, mu) * peekLorentz(force, mu); tmp_mu = peekLorentz(*this->ThinLinks, mu) * peekLorentz(force, mu);
pokeLorentz(SigmaTilde, tmp_mu, mu); pokeLorentz(SigmaTilde, tmp_mu, mu);
} }
double end = usecond(); double end = usecond();
double time = (end - start)/ 1e3; double time = (end - start)/ 1e3;
std::cout << GridLogMessage << " GaugeConfigurationMasked: Smeared Force chain rule took " << time << " ms" << std::endl; std::cout << GridLogMessage << " GaugeConfigurationMasked: Smeared Force chain rule took " << time << " ms" << std::endl;
} // if smearingLevels = 0 do nothing } // if smearingLevels = 0 do nothing
SigmaTilde=Gimpl::projectForce(SigmaTilde); // Ta
} }
}; };