portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Hadrons: Z2 source code fix

Browse Source
This commit is contained in:
Antonin Portelli 2016-05-01 14:49:45 -07:00
parent 5e4b58ac40
commit 7905c5b8e5
1 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
programs/Hadrons/MSource.cc
View File

@ -119,17 +119,15 @@ void MSource::execute(Environment &env)
int tb = std::stoi(par_.arguments[1]); int tb = std::stoi(par_.arguments[1]);
Lattice<iScalar<vInteger>> t(env.getGrid()); Lattice<iScalar<vInteger>> t(env.getGrid());
LatticeComplex eta(env.getGrid()); LatticeComplex eta(env.getGrid());
LatticeComplex z(env.getGrid()), shift(env.getGrid());
LatticeFermion phi(env.getGrid()); LatticeFermion phi(env.getGrid());
Complex shift(1., 1.);
LatticeCoordinate(t, Tp); LatticeCoordinate(t, Tp);
bernoulli(*env.get4dRng(), eta); bernoulli(*env.get4dRng(), eta);
shift = ComplexD(1., 1.); eta = (2.*eta - shift)*(1./::sqrt(2.));
z = zero; eta = where((t >= ta) and (t <= tb), eta, 0.*eta);
eta = (2.*eta - shift)*(1./::sqrt(2.)); *src_ = 1.;
eta = where((t >= ta) and (t <= tb), eta, z); *src_ = (*src_)*eta;
*src_ = 1.;
*src_ *= eta;
break; break;
} }