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

changed some naming conditions to resemble rare-kaons

Browse Source
This commit is contained in:
ferben 2019-12-05 16:56:54 +00:00
parent 1a5e562bde
commit 58b6a0d8d1
2 changed files with 47 additions and 47 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp
View File

@ -43,29 +43,29 @@ BEGIN_HADRONS_NAMESPACE
/* /*
* Sigma-to-nucleon 3-pt diagrams, eye topologies. * Sigma-to-nucleon 3-pt diagrams, eye topologies.
* *
* Schematics: qq_loop | * Schematics: qqLoop |
* /->-¬ | * /->-¬ |
* / \ | qs_xi G qd_xf * / \ | qsTi G qdTf
* \ / | /---->------*------>----¬ * \ / | /---->------*------>----¬
* qs_xi \ / qd_xf | / /-*-¬ \ * qsTi \ / qdTf | / /-*-¬ \
* /----->-----* *----->----¬ | / / G \ \ * /----->-----* *----->----¬ | / / G \ \
* * G G * | * \ / qq_loop * * * G G * | * \ / qqLoop *
* |\ /| | |\ \-<-/ /| * |\ /| | |\ \-<-/ /|
* | \ / | | | \ / | * | \ / | | | \ / |
* | \---------->---------/ | | | \----------->----------/ | * | \---------->---------/ | | | \----------->----------/ |
* \ qu_spec / | \ qu_spec / * \ quSpec / | \ quSpec /
* \ / | \ / * \ / | \ /
* \---------->---------/ | \----------->----------/ * \---------->---------/ | \----------->----------/
* qu_spec | qu_spec * quSpec | quSpec
* *
* analogously to the rare-kaon naming, the left diagram is named 'one-trace' and * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and
* the diagram on the right 'two-trace' * the diagram on the right 'two-trace'
* *
* Propagators: * Propagators:
* * qq_loop * * qqLoop
* * qu_spec, source at xi * * quSpec, source at ti
* * qd_xf, source at xf * * qdTf, source at tf
* * qs_xi, source at xi * * qsTi, source at ti
*/ */
BEGIN_MODULE_NAMESPACE(MContraction) BEGIN_MODULE_NAMESPACE(MContraction)
@ -73,11 +73,11 @@ class SigmaToNucleonEyePar: Serializable
{ {
public: public:
GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonEyePar, GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonEyePar,
std::string, qq_loop, std::string, qqLoop,
std::string, qu_spec, std::string, quSpec,
std::string, qd_xf, std::string, qdTf,
std::string, qs_xi, std::string, qsTi,
unsigned int, xf, unsigned int, tf,
std::string, sink, std::string, sink,
std::string, output); std::string, output);
}; };
@ -134,7 +134,7 @@ TSigmaToNucleonEye<FImpl>::TSigmaToNucleonEye(const std::string name)
template <typename FImpl> template <typename FImpl>
std::vector<std::string> TSigmaToNucleonEye<FImpl>::getInput(void) std::vector<std::string> TSigmaToNucleonEye<FImpl>::getInput(void)
{ {
std::vector<std::string> input = {par().qq_loop, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; std::vector<std::string> input = {par().qqLoop, par().quSpec, par().qdTf, par().qsTi, par().sink};
return input; return input;
} }
@ -175,17 +175,17 @@ void TSigmaToNucleonEye<FImpl>::execute(void)
r.info.gammaA_nucl = Id.g; r.info.gammaA_nucl = Id.g;
r.info.gammaB_nucl = GammaB.g; r.info.gammaB_nucl = GammaB.g;
auto &qq_loop = envGet(PropagatorField, par().qq_loop); auto &qqLoop = envGet(PropagatorField, par().qqLoop);
auto &qu_spec = envGet(SlicedPropagator, par().qu_spec); auto &quSpec = envGet(SlicedPropagator, par().quSpec);
auto &qd_xf = envGet(PropagatorField, par().qd_xf); auto &qdTf = envGet(PropagatorField, par().qdTf);
auto &qs_xi = envGet(PropagatorField, par().qs_xi); auto &qsTi = envGet(PropagatorField, par().qsTi);
auto qut = qu_spec[par().xf]; auto qut = quSpec[par().tf];
for (auto &G: Gamma::gall) for (auto &G: Gamma::gall)
{ {
r.info.gamma_H = G.g; r.info.gamma_H = G.g;
//Operator Q1, equivalent to the two-trace case in the rare-kaons module //Operator Q1, equivalent to the two-trace case in the rare-kaons module
c=Zero(); c=Zero();
BaryonUtils<FIMPL>::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q1",c); BaryonUtils<FIMPL>::Sigma_to_Nucleon_Eye(qqLoop,qut,qdTf,qsTi,G,GammaB,GammaB,"Q1",c);
sliceSum(c,buf,Tp); sliceSum(c,buf,Tp);
r.corr.clear(); r.corr.clear();
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)
@ -196,7 +196,7 @@ void TSigmaToNucleonEye<FImpl>::execute(void)
result.push_back(r); result.push_back(r);
//Operator Q2, equivalent to the one-trace case in the rare-kaons module //Operator Q2, equivalent to the one-trace case in the rare-kaons module
c=Zero(); c=Zero();
BaryonUtils<FIMPL>::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q2",c); BaryonUtils<FIMPL>::Sigma_to_Nucleon_Eye(qqLoop,qut,qdTf,qsTi,G,GammaB,GammaB,"Q2",c);
sliceSum(c,buf,Tp); sliceSum(c,buf,Tp);
r.corr.clear(); r.corr.clear();
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)

48
Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp
View File

@ -44,7 +44,7 @@ BEGIN_HADRONS_NAMESPACE
* Sigma-to-Nucleon 3-pt diagrams, non-eye topologies. * Sigma-to-Nucleon 3-pt diagrams, non-eye topologies.
* *
* Schematic: * Schematic:
* qs_xi qu_xf | qs_xi qd_xf * qsTi quTf | qsTi qdTf
* /-->--¬ /-->--¬ | /-->--¬ /-->--¬ * /-->--¬ /-->--¬ | /-->--¬ /-->--¬
* / \ / \ | / \ / \ * / \ / \ | / \ / \
* / \ / \ | / \ / \ * / \ / \ | / \ / \
@ -54,22 +54,22 @@ BEGIN_HADRONS_NAMESPACE
* | \ / \ /| | | \ / \ / | * | \ / \ /| | | \ / \ / |
* | \ / \ / | | | \ / \ / | * | \ / \ / | | | \ / \ / |
* | \ / \ / | | | \-->--/ \-->--/ | * | \ / \ / | | | \-->--/ \-->--/ |
* \ \-->--/ \-->--/ / | \ qu_xi qu_xf / * \ \-->--/ \-->--/ / | \ quTi quTf /
* \ qu_xi qd_xf / | \ / * \ quTi qdTf / | \ /
* \ / | \ / * \ / | \ /
* \--------->----------/ | \--------->-----------/ * \--------->----------/ | \--------->-----------/
* qu_spec | qu_spec * quSpec | quSpec
* *
* *
* analogously to the rare-kaon naming, the left diagram is named 'one-trace' and * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and
* the diagram on the right 'two-trace' * the diagram on the right 'two-trace'
* *
* Propagators: * Propagators:
* * qu_xi, source at xi * * quTi, source at ti
* * qu_xf, source at xf * * quTf, source at tf
* * qu_spec, source at xi * * quSpec, source at ti
* * qd_xf, source at xf * * qdTf, source at tf
* * qs_xi, source at xi * * qsTi, source at ti
*/ */
BEGIN_MODULE_NAMESPACE(MContraction) BEGIN_MODULE_NAMESPACE(MContraction)
@ -77,12 +77,12 @@ class SigmaToNucleonNonEyePar: Serializable
{ {
public: public:
GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonNonEyePar, GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonNonEyePar,
std::string, qu_xi, std::string, quTi,
std::string, qu_xf, std::string, quTf,
std::string, qu_spec, std::string, quSpec,
std::string, qd_xf, std::string, qdTf,
std::string, qs_xi, std::string, qsTi,
unsigned int, xf, unsigned int, tf,
std::string, sink, std::string, sink,
std::string, output); std::string, output);
}; };
@ -139,7 +139,7 @@ TSigmaToNucleonNonEye<FImpl>::TSigmaToNucleonNonEye(const std::string name)
template <typename FImpl> template <typename FImpl>
std::vector<std::string> TSigmaToNucleonNonEye<FImpl>::getInput(void) std::vector<std::string> TSigmaToNucleonNonEye<FImpl>::getInput(void)
{ {
std::vector<std::string> input = {par().qu_xi, par().qu_xf, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; std::vector<std::string> input = {par().quTi, par().quTf, par().quSpec, par().qdTf, par().qsTi, par().sink};
return input; return input;
} }
@ -180,18 +180,18 @@ void TSigmaToNucleonNonEye<FImpl>::execute(void)
r.info.gammaA_nucl = Id.g; r.info.gammaA_nucl = Id.g;
r.info.gammaB_nucl = GammaB.g; r.info.gammaB_nucl = GammaB.g;
auto &qu_xi = envGet(PropagatorField, par().qu_xi); auto &quTi = envGet(PropagatorField, par().quTi);
auto &qu_xf = envGet(PropagatorField, par().qu_xf); auto &quTf = envGet(PropagatorField, par().quTf);
auto &qu_spec = envGet(SlicedPropagator, par().qu_spec); auto &quSpec = envGet(SlicedPropagator, par().quSpec);
auto &qd_xf = envGet(PropagatorField, par().qd_xf); auto &qdTf = envGet(PropagatorField, par().qdTf);
auto &qs_xi = envGet(PropagatorField, par().qs_xi); auto &qsTi = envGet(PropagatorField, par().qsTi);
auto qut = qu_spec[par().xf]; auto qut = quSpec[par().tf];
for (auto &G: Gamma::gall) for (auto &G: Gamma::gall)
{ {
r.info.gamma_H = G.g; r.info.gamma_H = G.g;
//Operator Q1, equivalent to the two-trace case in the rare-kaons module //Operator Q1, equivalent to the two-trace case in the rare-kaons module
c=Zero(); c=Zero();
BaryonUtils<FIMPL>::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q1",c); BaryonUtils<FIMPL>::Sigma_to_Nucleon_NonEye(quTi,quTf,qut,qdTf,qsTi,G,GammaB,GammaB,"Q1",c);
sliceSum(c,buf,Tp); sliceSum(c,buf,Tp);
r.corr.clear(); r.corr.clear();
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)
@ -202,7 +202,7 @@ void TSigmaToNucleonNonEye<FImpl>::execute(void)
result.push_back(r); result.push_back(r);
//Operator Q2, equivalent to the one-trace case in the rare-kaons module //Operator Q2, equivalent to the one-trace case in the rare-kaons module
c=Zero(); c=Zero();
BaryonUtils<FIMPL>::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q2",c); BaryonUtils<FIMPL>::Sigma_to_Nucleon_NonEye(quTi,quTf,qut,qdTf,qsTi,G,GammaB,GammaB,"Q2",c);
sliceSum(c,buf,Tp); sliceSum(c,buf,Tp);
r.corr.clear(); r.corr.clear();
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)