portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 14:04:32 +00:00
Code Releases Activity

Hadrons: Momentum inserted as multiples of 2*pi/L

Browse Source
This commit is contained in:
Lanny91
2017-01-30 17:14:33 +00:00
parent c4d3672720
commit 651e1a7cbc
3 changed files with 9 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
extras/Hadrons/Modules/MSink/Wall.hpp
View File

@@ -4,7 +4,7 @@ Grid physics library, www.github.com/paboyle/Grid
Source file: extras/Hadrons/Modules/MSink/Wall.hpp
Copyright (C) 2016
Copyright (C) 2017
Author: Andrew Lawson <andrew.lawson1991@gmail.com>
@@ -131,9 +131,9 @@ void TWall<FImpl>::execute(void)
for(unsigned int mu = 0; mu < Nd; mu++)
{
LatticeCoordinate(coor, mu);
ph = ph + p[mu]*coor;
ph = ph + p[mu]*coor*((1./(env().getGrid()->_fdimensions[mu])));
}
ph = exp(-i*ph);
ph = exp(-2*M_PI*i*ph);
sliceSum<SitePropagator>(ph*q, prop, Tp);
}

5
extras/Hadrons/Modules/MSource/SeqGamma.hpp
View File

@@ -6,6 +6,7 @@ Source file: extras/Hadrons/Modules/MSource/SeqGamma.hpp
Copyright (C) 2015
Copyright (C) 2016
Copyright (C) 2017
Author: Antonin Portelli <antonin.portelli@me.com>
@@ -150,9 +151,9 @@ void TSeqGamma<FImpl>::execute(void)
for(unsigned int mu = 0; mu < env().getNd(); mu++)
{
LatticeCoordinate(coor, mu);
ph = ph + p[mu]*coor;
ph = ph + p[mu]*coor*((1./(env().getGrid()->_fdimensions[mu])));
}
ph = exp(i*ph);
ph = exp(2*M_PI*i*ph);
LatticeCoordinate(t, Tp);
src = where((t >= par().tA) and (t <= par().tB), g*ph*q, 0.*q);
}

6
extras/Hadrons/Modules/MSource/Wall.hpp
View File

@@ -4,7 +4,7 @@ Grid physics library, www.github.com/paboyle/Grid
Source file: extras/Hadrons/Modules/MSource/Wall.hpp
Copyright (C) 2016
Copyright (C) 2017
Author: Andrew Lawson <andrew.lawson1991@gmail.com>
@@ -132,9 +132,9 @@ void TWall<FImpl>::execute(void)
for(unsigned int mu = 0; mu < Nd; mu++)
{
LatticeCoordinate(coor, mu);
ph = ph + p[mu]*coor;
ph = ph + p[mu]*coor*((1./(env().getGrid()->_fdimensions[mu])));
}
ph = exp(i*ph);
ph = exp(2*M_PI*i*ph);
LatticeCoordinate(t, Tp);
src = 1.;
src = where((t == par().tW), src*ph, 0.*src);