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

Updated masked / fthmc smeared config container

Browse Source
This commit is contained in:
Peter Boyle 2023-06-01 06:23:02 -04:00
parent d2a8494044
commit ffd7301649
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -701,13 +701,16 @@ public:
{ {
if(domask) assert(Nsmear%(2*Nd)==0); // Or multiply by 8?? if(domask) assert(Nsmear%(2*Nd)==0); // Or multiply by 8??
// was resized in base class
assert(this->SmearedSet.size()==Nsmear);
GridRedBlackCartesian * UrbGrid; GridRedBlackCartesian * UrbGrid;
UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(_UGrid); UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(_UGrid);
LatticeComplex one(_UGrid); one = ComplexD(1.0,0.0); LatticeComplex one(_UGrid); one = ComplexD(1.0,0.0);
LatticeComplex tmp(_UGrid); LatticeComplex tmp(_UGrid);
for (unsigned int i = 0; i < this->smearingLevels; ++i) { for (unsigned int i = 0; i < this->smearingLevels; ++i) {
this->SmearedSet.push_back(*(new GaugeField(_UGrid)));
masks.push_back(*(new LatticeLorentzComplex(_UGrid))); masks.push_back(*(new LatticeLorentzComplex(_UGrid)));
if (domask) { if (domask) {