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

Hadrons: more contractions in the spectrum test

Browse Source
This commit is contained in:
Antonin Portelli 2016-12-05 17:41:58 +09:00
parent aa016f61b9
commit e6effcfd95
1 changed files with 19 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
tests/hadrons/Test_hadrons_mesons.cc
View File

@ -84,18 +84,26 @@ int main(int argc, char *argv[])
// propagators // propagators
Quark::Par quarkPar; Quark::Par quarkPar;
quarkPar.solver = "CG_" + flavour[i]; quarkPar.solver = "CG_" + flavour[i];
quarkPar.source = (flavour[i][0] == 'c') ? "z2" : "pt"; quarkPar.source = "pt";
application.createModule<Quark>("Q_" + flavour[i], quarkPar); application.createModule<Quark>("Qpt_" + flavour[i], quarkPar);
quarkPar.source = "z2";
application.createModule<Quark>("QZ2_" + flavour[i], quarkPar);
} }
for (unsigned int i = 0; i < flavour.size(); ++i) for (unsigned int i = 0; i < flavour.size(); ++i)
for (unsigned int j = i; j < flavour.size(); ++j) for (unsigned int j = i; j < flavour.size(); ++j)
{ {
MContraction::Meson::Par mesPar; MContraction::Meson::Par mesPar;
mesPar.output = "mesons/" + flavour[i] + flavour[j]; mesPar.output = "mesons/pt_" + flavour[i] + flavour[j];
mesPar.q1 = "Q_" + flavour[i]; mesPar.q1 = "Qpt_" + flavour[i];
mesPar.q2 = "Q_" + flavour[j]; mesPar.q2 = "Qpt_" + flavour[j];
application.createModule<MContraction::Meson>("meson_" application.createModule<MContraction::Meson>("meson_pt_"
+ flavour[i] + flavour[j],
mesPar);
mesPar.output = "mesons/Z2_" + flavour[i] + flavour[j];
mesPar.q1 = "QZ2_" + flavour[i];
mesPar.q2 = "QZ2_" + flavour[j];
application.createModule<MContraction::Meson>("meson_Z2_"
+ flavour[i] + flavour[j], + flavour[i] + flavour[j],
mesPar); mesPar);
} }
@ -105,12 +113,12 @@ int main(int argc, char *argv[])
{ {
MContraction::Baryon::Par barPar; MContraction::Baryon::Par barPar;
barPar.output = "baryons/" + flavour[i] + flavour[j] + flavour[k]; barPar.output = "baryons/pt_" + flavour[i] + flavour[j] + flavour[k];
barPar.q1 = "Q_" + flavour[i]; barPar.q1 = "Qpt_" + flavour[i];
barPar.q2 = "Q_" + flavour[j]; barPar.q2 = "Qpt_" + flavour[j];
barPar.q3 = "Q_" + flavour[k]; barPar.q3 = "Qpt_" + flavour[k];
application.createModule<MContraction::Baryon>( application.createModule<MContraction::Baryon>(
"baryon_" + flavour[i] + flavour[j] + flavour[k], barPar); "baryon_pt_" + flavour[i] + flavour[j] + flavour[k], barPar);
} }
// execution // execution