portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Added momentum choice for wall source

Browse Source
This commit is contained in:
Lanny91 2017-01-13 15:54:19 +00:00
parent f22b79da8f
commit b7f90aa011
1 changed files with 32 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -7,6 +7,19 @@
BEGIN_HADRONS_NAMESPACE BEGIN_HADRONS_NAMESPACE
/*
Wall source
-----------------------------
* src_x = theta(x_3 - tA) * theta(tB - x_3) * exp(i x.mom)
* options:
- tA: begin timeslice (integer)
- tB: end timeslice (integer)
- mom: momentum insertion, space-separated float sequence (e.g ".1 .2 1. 0.")
*/
/****************************************************************************** /******************************************************************************
* Wall * * Wall *
******************************************************************************/ ******************************************************************************/
@ -16,7 +29,8 @@ class WallPar: Serializable
{ {
public: public:
GRID_SERIALIZABLE_CLASS_MEMBERS(WallPar, GRID_SERIALIZABLE_CLASS_MEMBERS(WallPar,
unsigned int, tW); unsigned int, tW,
std::string, mom);
}; };
template <typename FImpl> template <typename FImpl>
@ -76,14 +90,27 @@ void TWall<FImpl>::setup(void)
// execution /////////////////////////////////////////////////////////////////// // execution ///////////////////////////////////////////////////////////////////
template <typename FImpl> template <typename FImpl>
void TWall<FImpl>::execute(void) void TWall<FImpl>::execute(void)
{ {
Lattice<iScalar<vInteger>> t(env().getGrid()); LOG(Message) << "Generating wall source at t = " << par().tW
<< " with momentum " << par().mom << std::endl;
LOG(Message) << "Generating wall source at t = " << par().tW << std::endl;
PropagatorField &src = *env().template createLattice<PropagatorField>(getName()); PropagatorField &src = *env().template createLattice<PropagatorField>(getName());
Lattice<iScalar<vInteger>> t(env().getGrid());
LatticeComplex ph(env().getGrid()), coor(env().getGrid());
std::vector<Real> p;
Complex i(0.0,1.0);
p = strToVec<Real>(par().mom);
ph = zero;
for(unsigned int mu = 0; mu < Nd; mu++)
{
LatticeCoordinate(coor, mu);
ph = ph + p[mu]*coor;
}
ph = exp(i*ph);
LatticeCoordinate(t, Tp); LatticeCoordinate(t, Tp);
src = 1.; src = 1.;
src = where((t == par().tW), src, 0.*src); src = where((t == par().tW), src*ph, 0.*src);
} }
END_MODULE_NAMESPACE END_MODULE_NAMESPACE