107 void mesons2(
const LatticePropagator& quark_prop_1,
108 const LatticePropagator& quark_prop_2,
121 LatticePropagator anti_quark_prop = Gamma(
G5) * quark_prop_2 * Gamma(
G5);
131 XMLArrayWriter xml_gamma(xml,Ns*Ns);
132 push(xml_gamma, xml_group);
134 for (
int gamma_value=0; gamma_value < (Ns*Ns); ++gamma_value)
137 write(xml_gamma,
"gamma_value", gamma_value);
140 LatticeComplex corr_fn;
141 corr_fn = trace(adj(anti_quark_prop) * (Gamma(gamma_value) *
142 quark_prop_1 * Gamma(gamma_value)));
144 multi2d<DComplex> hsum;
145 hsum = phases.
sft(corr_fn);
148 XMLArrayWriter xml_sink_mom(xml_gamma,phases.
numMom());
149 push(xml_sink_mom,
"momenta");
151 for (
int sink_mom_num=0; sink_mom_num < phases.
numMom(); ++sink_mom_num)
154 write(xml_sink_mom,
"sink_mom_num", sink_mom_num);
155 write(xml_sink_mom,
"sink_mom", phases.
numToMom(sink_mom_num));
157 multi1d<DComplex> mesprop(length);
158 for (
int t=0;
t < length; ++
t)
160 int t_eff = (
t -
t0 + length) % length;
161 mesprop[t_eff] = hsum[sink_mom_num][
t];
164 write(xml_sink_mom,
"mesprop", mesprop);
Primary include file for CHROMA library code.
Fourier transform phase factor support.
int numSubsets() const
Number of subsets - length in decay direction.
multi1d< int > numToMom(int mom_num) const
Convert momenta id to actual array of momenta.
multi2d< DComplex > sft(const LatticeComplex &cf) const
Do a sumMulti(cf*phases,getSet())
int numMom() const
Number of momenta.
void write(XMLWriter &xml, const std::string &path, const AsqtadFermActParams ¶m)
Writer parameters.
Asqtad Staggered-Dirac operator.
push(xml_out,"Condensates")
void mesons2(const LatticePropagator &quark_prop_1, const LatticePropagator &quark_prop_2, const SftMom &phases, int t0, XMLWriter &xml, const std::string &xml_group)
Meson 2-pt functions.
Fourier transform phase factor support.