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

Gamma input for meson contraction with round brackets

Browse Source
This commit is contained in:
Vera Guelpers 2018-03-12 18:02:12 +00:00
parent 70ec2faa98
commit 8462bbfe63
2 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
extras/Hadrons/Modules/MContraction/Meson.hpp
View File

@ -45,8 +45,8 @@ BEGIN_HADRONS_NAMESPACE
- q1: input propagator 1 (string) - q1: input propagator 1 (string)
- q2: input propagator 2 (string) - q2: input propagator 2 (string)
- gammas: gamma products to insert at sink & source, pairs of gamma matrices - gammas: gamma products to insert at sink & source, pairs of gamma matrices
(space-separated strings) in angled brackets (i.e. <g_sink g_src>), (space-separated strings) in round brackets (i.e. (g_sink g_src)),
in a sequence (e.g. "<Gamma5 Gamma5><Gamma5 GammaT>"). in a sequence (e.g. "(Gamma5 Gamma5)(Gamma5 GammaT)").
Special values: "all" - perform all possible contractions. Special values: "all" - perform all possible contractions.
- sink: module to compute the sink to use in contraction (string). - sink: module to compute the sink to use in contraction (string).

8
tests/hadrons/Test_QED.cc
View File

@ -189,7 +189,7 @@ int main(int argc, char *argv[])
mesPar.output = "QED/pt_" + flavour[i] + flavour[j]; mesPar.output = "QED/pt_" + flavour[i] + flavour[j];
mesPar.q1 = "Qpt_" + flavour[i]; mesPar.q1 = "Qpt_" + flavour[i];
mesPar.q2 = "Qpt_" + flavour[j]; mesPar.q2 = "Qpt_" + flavour[j];
mesPar.gammas = "<Gamma5 Gamma5>"; mesPar.gammas = "(Gamma5 Gamma5)";
mesPar.sink = "sink"; mesPar.sink = "sink";
application.createModule<MContraction::Meson>("meson_pt_" application.createModule<MContraction::Meson>("meson_pt_"
+ flavour[i] + flavour[j], + flavour[i] + flavour[j],
@ -203,7 +203,7 @@ int main(int argc, char *argv[])
+ flavour[i] + "__" + flavour[j]; + flavour[i] + "__" + flavour[j];
mesPar_seq_T.q1 = "Qpt_" + flavour[i] + "_seq_T" + flavour[i]; mesPar_seq_T.q1 = "Qpt_" + flavour[i] + "_seq_T" + flavour[i];
mesPar_seq_T.q2 = "Qpt_" + flavour[j]; mesPar_seq_T.q2 = "Qpt_" + flavour[j];
mesPar_seq_T.gammas = "<Gamma5 Gamma5>"; mesPar_seq_T.gammas = "(Gamma5 Gamma5)";
mesPar_seq_T.sink = "sink"; mesPar_seq_T.sink = "sink";
application.createModule<MContraction::Meson>("meson_tadpole_pt_" + application.createModule<MContraction::Meson>("meson_tadpole_pt_" +
flavour[i] + "_seq_T" flavour[i] + "_seq_T"
@ -219,7 +219,7 @@ int main(int argc, char *argv[])
+ flavour[j]; + flavour[j];
mesPar_seq_E.q1 = "Qpt_" + flavour[i] + "_seq_V_ph_" + flavour[i]; mesPar_seq_E.q1 = "Qpt_" + flavour[i] + "_seq_V_ph_" + flavour[i];
mesPar_seq_E.q2 = "Qpt_" + flavour[j] + "_seq_V_ph_" + flavour[j]; mesPar_seq_E.q2 = "Qpt_" + flavour[j] + "_seq_V_ph_" + flavour[j];
mesPar_seq_E.gammas = "<Gamma5 Gamma5>"; mesPar_seq_E.gammas = "(Gamma5 Gamma5)";
mesPar_seq_E.sink = "sink"; mesPar_seq_E.sink = "sink";
application.createModule<MContraction::Meson>("meson_exchange_pt_" application.createModule<MContraction::Meson>("meson_exchange_pt_"
+ flavour[i] + "_seq_V_ph_" + flavour[i] + flavour[i] + "_seq_V_ph_" + flavour[i]
@ -236,7 +236,7 @@ int main(int argc, char *argv[])
mesPar_seq_S.q1 = "Qpt_" + flavour[i] + "_seq_V_ph_" + flavour[i] mesPar_seq_S.q1 = "Qpt_" + flavour[i] + "_seq_V_ph_" + flavour[i]
+ "_seq_V_ph_" + flavour[i]; + "_seq_V_ph_" + flavour[i];
mesPar_seq_S.q2 = "Qpt_" + flavour[j]; mesPar_seq_S.q2 = "Qpt_" + flavour[j];
mesPar_seq_S.gammas = "<Gamma5 Gamma5>"; mesPar_seq_S.gammas = "(Gamma5 Gamma5)";
mesPar_seq_S.sink = "sink"; mesPar_seq_S.sink = "sink";
application.createModule<MContraction::Meson>("meson_selfenergy_pt_" application.createModule<MContraction::Meson>("meson_selfenergy_pt_"
+ flavour[i] + "_seq_V_ph_" + flavour[i] + "_seq_V_ph_"