From 421a4395af5669b4fabada96437c481b14b1d8c9 Mon Sep 17 00:00:00 2001 From: ferben Date: Thu, 21 Nov 2019 17:25:37 +0000 Subject: [PATCH 01/14] Sigma to Nucleon contractions --- Grid/qcd/utils/BaryonUtils.h | 321 +++++++++++++++++- Hadrons/Global.hpp | 2 +- Hadrons/Modules.hpp | 129 +++---- .../Modules/MContraction/SigmaToNucleonEye.cc | 7 + .../MContraction/SigmaToNucleonEye.hpp | 219 ++++++++++++ .../MContraction/SigmaToNucleonNonEye.cc | 7 + .../MContraction/SigmaToNucleonNonEye.hpp | 225 ++++++++++++ Hadrons/Modules/MContraction/WeakEye3pt.hpp | 2 +- .../Modules/MContraction/WeakNonEye3pt.hpp | 2 +- Hadrons/Modules/MSource/MomentumPhase.hpp | 2 +- Hadrons/modules.inc | 256 +++++++------- tests/hadrons/Test_sigma_to_nucleon.cc | 154 +++++++++ 12 files changed, 1125 insertions(+), 201 deletions(-) create mode 100644 Hadrons/Modules/MContraction/SigmaToNucleonEye.cc create mode 100644 Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp create mode 100644 Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc create mode 100644 Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp create mode 100644 tests/hadrons/Test_sigma_to_nucleon.cc diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 73d41422..7fcc4502 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -41,9 +41,12 @@ public: typedef typename FImpl::SitePropagator pobj; typedef typename ComplexField::vector_object vobj; - + + typedef Lattice> SpinMatrixField; + typedef typename SpinMatrixField::vector_object sobj; + static constexpr int epsilon[6][3] = {{0,1,2},{1,2,0},{2,0,1},{0,2,1},{2,1,0},{1,0,2}}; - static constexpr Complex epsilon_sgn[6]= {1,1,1,-1,-1,-1}; + static constexpr int epsilon_sgn[6]= {1,1,1,-1,-1,-1}; private: template @@ -81,12 +84,70 @@ public: const char * quarks_right, const int parity, robj &result); + private: + template + static void Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result); + template + static void Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, + const mobj &Du_xf, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result); + + + template + static void Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result); + template + static void Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, + const mobj &Du_xf, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result); + public: + template + static void Sigma_to_Nucleon(const PropagatorField &qq_xi, + const PropagatorField &qq_xf, + const mobj &Du_spec, + const PropagatorField &qd_xf, + const PropagatorField &qs_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + const std::string op, + SpinMatrixField &stn_corr); }; template constexpr int BaryonUtils::epsilon[6][3]; template -constexpr Complex BaryonUtils::epsilon_sgn[6]; +constexpr int BaryonUtils::epsilon_sgn[6]; template template @@ -123,7 +184,7 @@ void BaryonUtils::baryon_site(const mobj &D1, for (int alpha_right=0; alpha_right(epsilon_sgn[ie_left] * epsilon_sgn[ie_right]) * pD1()(gamma_left,gamma_left)(c_right,c_left)*D2g()(alpha_right,beta_left)(a_right,a_left)*gD3()(alpha_right,beta_left)(b_right,b_left); }}} } //This is the \delta_{456}^{231} part @@ -132,7 +193,7 @@ void BaryonUtils::baryon_site(const mobj &D1, for (int alpha_right=0; alpha_right(epsilon_sgn[ie_left] * epsilon_sgn[ie_right]) * pD1g()(gamma_left,beta_left)(c_right,a_left)*D2()(alpha_right,beta_left)(a_right,b_left)*gD3()(alpha_right,gamma_left)(b_right,c_left); }}} } //This is the \delta_{456}^{312} part @@ -141,7 +202,7 @@ void BaryonUtils::baryon_site(const mobj &D1, for (int alpha_right=0; alpha_right(epsilon_sgn[ie_left] * epsilon_sgn[ie_right]) * pD1()(gamma_left,beta_left)(c_right,b_left)*D2()(alpha_right,gamma_left)(a_right,c_left)*gD3g()(alpha_right,beta_left)(b_right,a_left); }}} } //This is the \delta_{456}^{132} part @@ -150,7 +211,7 @@ void BaryonUtils::baryon_site(const mobj &D1, for (int alpha_right=0; alpha_right(epsilon_sgn[ie_left] * epsilon_sgn[ie_right]) * pD1()(gamma_left,gamma_left)(c_right,c_left)*D2()(alpha_right,beta_left)(a_right,b_left)*gD3g()(alpha_right,beta_left)(b_right,a_left); }}} } //This is the \delta_{456}^{321} part @@ -159,7 +220,7 @@ void BaryonUtils::baryon_site(const mobj &D1, for (int alpha_right=0; alpha_right(epsilon_sgn[ie_left] * epsilon_sgn[ie_right]) * pD1()(gamma_left,beta_left)(c_right,b_left)*D2g()(alpha_right,beta_left)(a_right,a_left)*gD3()(alpha_right,gamma_left)(b_right,c_left); }}} } //This is the \delta_{456}^{213} part @@ -168,7 +229,7 @@ void BaryonUtils::baryon_site(const mobj &D1, for (int alpha_right=0; alpha_right(epsilon_sgn[ie_left] * epsilon_sgn[ie_right]) * pD1g()(gamma_left,beta_left)(c_right,a_left)*D2()(alpha_right,gamma_left)(a_right,c_left)*gD3()(alpha_right,beta_left)(b_right,b_left); }}} } } @@ -249,4 +310,246 @@ void BaryonUtils::ContractBaryons_Sliced(const mobj &D1, result=Zero(); baryon_site(D1,D2,D3,GammaA_left,GammaB_left,GammaA_right,GammaB_right,parity,wick_contraction,result); } + +/*********************************************************************** + * End of Baryon 2pt-function code. * + * * + * The following code is for Sigma -> N rare hypeon decays * + **********************************************************************/ + +/* Dq_loop is a quark line from x_H to x_H + * Du_spec is a quark line from x_i to x_f + * Dd_xf is a quark line from x_H to x_f + * Ds_xi is a quark line from x_H to x_i */ +template +template +void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result) +{ + + Gamma g5(Gamma::Algebra::Gamma5); + + auto DuG = Du_spec * GammaB_nucl; + // Gamma^B * Ds * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) + auto GDsGDd = GammaB_sigma * Ds_xi * Gamma_H * g5 * adj(Dd_xf) * g5; + // Dq_loop * \gamma_\mu^L + auto DqG = Dq_loop * Gamma_H; + + for (int ie_n=0; ie_n < 6 ; ie_n++){ + int a_n = epsilon[ie_n][0]; //a + int b_n = epsilon[ie_n][1]; //b + int c_n = epsilon[ie_n][2]; //c + for (int ie_s=0; ie_s < 6 ; ie_s++){ + int a_s = epsilon[ie_s][0]; //a' + int b_s = epsilon[ie_s][1]; //b' + int c_s = epsilon[ie_s][2]; //c' + for (int j=0; j(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DqG()(tau2,tau2)(j,j) * DuG()(alpha_s, beta_n)(a_s,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); + result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DqG()(tau2,tau2)(j,j) * DuG()(gamma_s, beta_n)(c_s,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); + }}}}}} + } + } +} + +/* Du_xi is a quark line from x_H to x_i + * Du_xf is a quark line from x_H to x_f + * Du_spec is a quark line from x_i to x_f + * Dd_xf is a quark line from x_H to x_f + * Ds_xi is a quark line from x_H to x_i */ +template +template +void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, + const mobj &Du_xf, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result) +{ + + Gamma g5(Gamma::Algebra::Gamma5); + + auto DuG = Du_spec * GammaB_nucl; + auto adjDu = g5 * adj(Du_xf) * g5; + auto adjDuG = adjDu * GammaB_nucl; + // Gamma^B * Ds * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) + auto GDsGDd = GammaB_sigma * Ds_xi * Gamma_H * g5 * adj(Dd_xf) * g5; + // Dq_loop * \gamma_\mu^L + auto DuGH = Du_xi * Gamma_H; + + for (int ie_n=0; ie_n < 6 ; ie_n++){ + int a_n = epsilon[ie_n][0]; //a + int b_n = epsilon[ie_n][1]; //b + int c_n = epsilon[ie_n][2]; //c + for (int ie_s=0; ie_s < 6 ; ie_s++){ + int a_s = epsilon[ie_s][0]; //a' + int b_s = epsilon[ie_s][1]; //b' + int c_s = epsilon[ie_s][2]; //c' + for (int j=0; j(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(alpha_s,tau2)(a_s,j) * DuG()(gamma_s, beta_n)(c_s,a_n) * adjDu()(tau2,gamma_n)(j,c_n); + result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(gamma_s,tau2)(c_s,j) * adjDuG()(tau2, beta_n)(j,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); + result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(alpha_s,tau2)(a_s,j) * adjDuG()(tau2, beta_n)(j,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); + result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(gamma_s,tau2)(c_s,j) * DuG()(alpha_s, beta_n)(a_s,a_n) * adjDu()(tau2,gamma_n)(j,c_n); + }}}}}} + } + } +} + +//Equivalent to "One-trace" +template +template +void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result) +{ + + Gamma g5(Gamma::Algebra::Gamma5); + + auto DuG = Du_spec * GammaB_nucl; + // Gamma^B * Ds * \gamma_\mu^L + auto GDsG = GammaB_sigma * Ds_xi * Gamma_H; + // Dq_loop * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) + auto DqGDd = Dq_loop * Gamma_H * g5 * adj(Dd_xf) * g5; + + for (int ie_n=0; ie_n < 6 ; ie_n++){ + int a_n = epsilon[ie_n][0]; //a + int b_n = epsilon[ie_n][1]; //b + int c_n = epsilon[ie_n][2]; //c + for (int ie_s=0; ie_s < 6 ; ie_s++){ + int a_s = epsilon[ie_s][0]; //a' + int b_s = epsilon[ie_s][1]; //b' + int c_s = epsilon[ie_s][2]; //c' + for (int i=0; i(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DqGDd()(tau,beta_n)(i,b_n) * DuG()(alpha_s, beta_n)(a_s,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); + result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DqGDd()(tau,beta_n)(i,b_n) * DuG()(gamma_s, beta_n)(c_s,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); + }}}}}} + } + } +} + +template +template +void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, + const mobj &Du_xf, + const mobj2 &Du_spec, + const mobj &Dd_xf, + const mobj &Ds_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + robj &result) +{ + + Gamma g5(Gamma::Algebra::Gamma5); + + auto DuG = Du_spec * GammaB_nucl; + auto adjDu = g5 * adj(Du_xf) * g5; + auto adjDuG = adjDu * GammaB_nucl; + // Gamma^B * Ds * \gamma_\mu^L + auto GDsG = GammaB_sigma * Ds_xi * Gamma_H; + // Du * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) + auto DuGDd = Du_xi * Gamma_H * g5 * adj(Dd_xf) * g5; + + for (int ie_n=0; ie_n < 6 ; ie_n++){ + int a_n = epsilon[ie_n][0]; //a + int b_n = epsilon[ie_n][1]; //b + int c_n = epsilon[ie_n][2]; //c + for (int ie_s=0; ie_s < 6 ; ie_s++){ + int a_s = epsilon[ie_s][0]; //a' + int b_s = epsilon[ie_s][1]; //b' + int c_s = epsilon[ie_s][2]; //c' + for (int i=0; i(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(alpha_s,beta_n)(a_s,b_n) * DuG()(gamma_s, beta_n)(c_s,a_n) * adjDu()(tau,gamma_n)(i,c_n); + result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(gamma_s,beta_n)(c_s,b_n) * adjDuG()(tau, beta_n)(i,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); + result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(alpha_s,beta_n)(a_s,b_n) * adjDuG()(tau, beta_n)(i,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); + result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(gamma_s,beta_n)(c_s,b_n) * DuG()(alpha_s, beta_n)(a_s,a_n) * adjDu()(tau,gamma_n)(i,c_n); + }}}}}} + } + } +} + + +template +template +void BaryonUtils::Sigma_to_Nucleon(const PropagatorField &qq_xi, + const PropagatorField &qq_xf, + const mobj &Du_spec, + const PropagatorField &qd_xf, + const PropagatorField &qs_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + const std::string op, + SpinMatrixField &stn_corr) +{ + assert(parity==1 || parity == -1 && "Parity must be +1 or -1"); + + GridBase *grid = qs_xi.Grid(); + + auto vcorr= stn_corr.View(); + auto vq_xi = qq_xi.View(); + auto vq_xf = qq_xf.View(); + auto vd_xf = qd_xf.View(); + auto vs_xi = qs_xi.View(); + + // accelerator_for(ss, grid->oSites(), grid->Nsimd(), { + thread_for(ss,grid->oSites(),{ + auto Dq_xi = vq_xi[ss]; + auto Dq_xf = vq_xf[ss]; + auto Dd_xf = vd_xf[ss]; + auto Ds_xi = vs_xi[ss]; + sobj result=Zero(); + if(op == "Q1_NonEye"){ + Sigma_to_Nucleon_Q1_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + } else if(op == "Q1_Eye"){ + Sigma_to_Nucleon_Q1_Eye_site(Dq_xi,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + } else if(op == "Q2_NonEye"){ + Sigma_to_Nucleon_Q2_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + } else if(op == "Q2_Eye"){ + Sigma_to_Nucleon_Q2_Eye_site(Dq_xi,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + } else { + assert(0 && "Weak Operator not correctly specified"); + } + vcorr[ss] = result; + } );//end loop over lattice sites +} + NAMESPACE_END(Grid); diff --git a/Hadrons/Global.hpp b/Hadrons/Global.hpp index 9b592721..b82c6288 100644 --- a/Hadrons/Global.hpp +++ b/Hadrons/Global.hpp @@ -272,7 +272,7 @@ struct Correlator: Serializable { GRID_SERIALIZABLE_CLASS_MEMBERS(ARG(Correlator), Metadata, info, - std::vector, corr); + std::vector, corr); }; END_HADRONS_NAMESPACE diff --git a/Hadrons/Modules.hpp b/Hadrons/Modules.hpp index c041afa0..ebe82910 100644 --- a/Hadrons/Modules.hpp +++ b/Hadrons/Modules.hpp @@ -1,73 +1,76 @@ -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include #include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include +#include #include -#include -#include +#include #include #include #include -#include -#include -#include -#include +#include +#include +#include #include #include +#include +#include #include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include #include #include -#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc b/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc new file mode 100644 index 00000000..4e938309 --- /dev/null +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc @@ -0,0 +1,7 @@ +#include + +using namespace Grid; +using namespace Hadrons; +using namespace MContraction; + +template class Grid::Hadrons::MContraction::TSigmaToNucleonEye; diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp new file mode 100644 index 00000000..e4189f7c --- /dev/null +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -0,0 +1,219 @@ +/************************************************************************************* + +Grid physics library, www.github.com/paboyle/Grid + +Source file: Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp + +Copyright (C) 2015-2019 + +Author: Antonin Portelli +Author: Felix Erben + +This program is free software; you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation; either version 2 of the License, or +(at your option) any later version. + +This program is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + +You should have received a copy of the GNU General Public License along +with this program; if not, write to the Free Software Foundation, Inc., +51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. + +See the full license in the file "LICENSE" in the top level distribution directory +*************************************************************************************/ +/* END LEGAL */ + +#ifndef Hadrons_MContraction_SigmaToNucleonEye_hpp_ +#define Hadrons_MContraction_SigmaToNucleonEye_hpp_ + +#include +#include +#include +#include + +BEGIN_HADRONS_NAMESPACE + +/****************************************************************************** + * SigmaToNucleonEye * + ******************************************************************************/ +/* + * Sigma-to-nucleon 3-pt diagrams, eye topologies. + * + * Schematics: qq_loop | + * /-<-¬ | + * / \ | qs_xi G qd_xf + * \ / | /----<------*------<----¬ + * qs_xi \ / qd_xf | / /-*-¬ \ + * /-----<-----* *-----<----¬ | / / G \ \ + * * G G * | * \ / qq_loop * + * |\ /| | |\ \->-/ /| + * | \ / | | | \ / | + * | \---------->---------/ | | | \----------->----------/ | + * \ qu_spec / | \ qu_spec / + * \ / | \ / + * \---------->---------/ | \----------->----------/ + * qu_spec | qu_spec + * + * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and + * the diagram on the right 'two-trace' + */ +BEGIN_MODULE_NAMESPACE(MContraction) + +class SigmaToNucleonEyePar: Serializable +{ +public: + GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonEyePar, + std::string, qq_loop, + std::string, qu_spec, + std::string, qd_xf, + std::string, qs_xi, + unsigned int, tOut, + std::string, parity, + std::string, sink, + std::string, output); +}; + +template +class TSigmaToNucleonEye: public Module +{ +public: + FERM_TYPE_ALIASES(FImpl1, 1); + FERM_TYPE_ALIASES(FImpl2, 2); + FERM_TYPE_ALIASES(FImpl3, 3); + FERM_TYPE_ALIASES(FImpl3, 4); + BASIC_TYPE_ALIASES(ScalarImplCR, Scalar); + SINK_TYPE_ALIASES(Scalar); + typedef typename SpinMatrixField1::vector_object::scalar_object SpinMatrix; + class Metadata: Serializable + { + public: + GRID_SERIALIZABLE_CLASS_MEMBERS(Metadata, + Gamma::Algebra, gamma_H, + Gamma::Algebra, gammaA_sigma, + Gamma::Algebra, gammaB_sigma, + Gamma::Algebra, gammaA_nucl, + Gamma::Algebra, gammaB_nucl, + int, trace, + int, parity); + }; + typedef Correlator Result; +public: + // constructor + TSigmaToNucleonEye(const std::string name); + // destructor + virtual ~TSigmaToNucleonEye(void) {}; + // dependency relation + virtual std::vector getInput(void); + virtual std::vector getOutput(void); +protected: + // setup + virtual void setup(void); + // execution + virtual void execute(void); + // Which gamma algebra was specified + Gamma::Algebra al; +}; + +MODULE_REGISTER_TMP(SigmaToNucleonEye, ARG(TSigmaToNucleonEye), MContraction); + +/****************************************************************************** + * TSigmaToNucleonEye implementation * + ******************************************************************************/ +// constructor ///////////////////////////////////////////////////////////////// +template +TSigmaToNucleonEye::TSigmaToNucleonEye(const std::string name) +: Module(name) +{} + +// dependencies/products /////////////////////////////////////////////////////// +template +std::vector TSigmaToNucleonEye::getInput(void) +{ + std::vector input = {par().qq_loop, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; + + return input; +} + +template +std::vector TSigmaToNucleonEye::getOutput(void) +{ + std::vector out = {}; + + return out; +} + +// setup /////////////////////////////////////////////////////////////////////// +template +void TSigmaToNucleonEye::setup(void) +{ + envTmpLat(SpinMatrixField1, "c"); +} + +// execution /////////////////////////////////////////////////////////////////// +template +void TSigmaToNucleonEye::execute(void) +{ + const int parity {par().parity.size()>0 ? std::stoi(par().parity) : 1}; + const Gamma GammaB(Gamma::Algebra::SigmaXZ); // C*gamma_5 + const Gamma Id(Gamma::Algebra::Identity); // C*gamma_5 + + LOG(Message) << "Computing sigma-to-nucleon contractions '" << getName() << "'" << std::endl; + LOG(Message) << "' with (Gamma^A,Gamma^B)_sigma = ( Identity, C*gamma_5 ) and (Gamma^A,Gamma^B)_nucl = ( Identity, C*gamma_5 )" << std::endl; + LOG(Message) << "and parity " << parity << " using sink " << par().sink << "." << std::endl; + + envGetTmp(SpinMatrixField1, c); + std::vector buf; + + std::vector result; + Result r; + r.info.parity = parity; + r.info.gammaA_sigma = Id.g; + r.info.gammaB_sigma = GammaB.g; + r.info.gammaA_nucl = Id.g; + r.info.gammaB_nucl = GammaB.g; + + auto &qq_loop = envGet(PropagatorField1, par().qq_loop); + auto &qu_spec = envGet(SlicedPropagator2, par().qu_spec); + auto &qd_xf = envGet(PropagatorField3, par().qd_xf); + auto &qs_xi = envGet(PropagatorField4, par().qs_xi); + auto qut = qu_spec[par().tOut]; + for (auto &G: Gamma::gall) + { + r.info.gamma_H = G.g; + //Operator Q1, equivalent to the two-trace case in the rare-kaons module + c=Zero(); + BaryonUtils::Sigma_to_Nucleon(qq_loop,qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1_Eye",c); + sliceSum(c,buf,Tp); + r.corr.clear(); + for (unsigned int t = 0; t < buf.size(); ++t) + { + r.corr.push_back(buf[t]); + } + r.info.trace = 2; + result.push_back(r); + //Operator Q2, equivalent to the one-trace case in the rare-kaons module + c=Zero(); + BaryonUtils::Sigma_to_Nucleon(qq_loop,qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2_Eye",c); + sliceSum(c,buf,Tp); + r.corr.clear(); + for (unsigned int t = 0; t < buf.size(); ++t) + { + r.corr.push_back(buf[t]); + } + r.info.trace = 1; + result.push_back(r); + } + + saveResult(par().output, "StN_Eye", result); + +} + +END_MODULE_NAMESPACE + +END_HADRONS_NAMESPACE + +#endif // Hadrons_MContraction_SigmaToNucleonEye_hpp_ diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc new file mode 100644 index 00000000..7d5ee247 --- /dev/null +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc @@ -0,0 +1,7 @@ +#include + +using namespace Grid; +using namespace Hadrons; +using namespace MContraction; + +template class Grid::Hadrons::MContraction::TSigmaToNucleonNonEye; diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp new file mode 100644 index 00000000..15dc5eaa --- /dev/null +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -0,0 +1,225 @@ +/************************************************************************************* + +Grid physics library, www.github.com/paboyle/Grid + +Source file: Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp + +Copyright (C) 2015-2019 + +Author: Antonin Portelli +Author: Felix Erben + +This program is free software; you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation; either version 2 of the License, or +(at your option) any later version. + +This program is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + +You should have received a copy of the GNU General Public License along +with this program; if not, write to the Free Software Foundation, Inc., +51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. + +See the full license in the file "LICENSE" in the top level distribution directory +*************************************************************************************/ +/* END LEGAL */ + +#ifndef Hadrons_MContraction_SigmaToNucleonNonEye_hpp_ +#define Hadrons_MContraction_SigmaToNucleonNonEye_hpp_ + +#include +#include +#include +#include + +BEGIN_HADRONS_NAMESPACE + +/****************************************************************************** + * SigmaToNucleonNonEye * + ******************************************************************************/ +/* + * Sigma-to-Nucleon 3-pt diagrams, non-eye topologies. + * + * Schematic: + * qs_xi qd_xf | qs_xi qd_xf + * /--<--¬ /--<--¬ | /--<--¬ /--<--¬ + * / \ / \ | / \ / \ + * / \ / \ | / \ / \ + * / \ / \ | / \ / \ + * * * G * | * G * * G * + * |\ * G | | |\ / \ /| + * | \ / \ /| | | \ / \ / | + * | \ / \ / | | | \ / \ / | + * | \ / \ / | | | \-->--/ \-->--/ | + * \ \-->--/ \-->--/ / | \ qu_xi qu_xf / + * \ qu_xf qu_xf / | \ / + * \ / | \ / + * \--------->----------/ | \--------->-----------/ + * qu_spec | qu_spec + * + * + * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and + * the diagram on the right 'two-trace' + * + */ +BEGIN_MODULE_NAMESPACE(MContraction) + +class SigmaToNucleonNonEyePar: Serializable +{ +public: + GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonNonEyePar, + std::string, qu_xi, + std::string, qu_xf, + std::string, qu_spec, + std::string, qd_xf, + std::string, qs_xi, + unsigned int, tOut, + std::string, parity, + std::string, sink, + std::string, output); +}; + +template +class TSigmaToNucleonNonEye: public Module +{ +public: + FERM_TYPE_ALIASES(FImpl1, 1); + FERM_TYPE_ALIASES(FImpl2, 2); + FERM_TYPE_ALIASES(FImpl3, 3); + FERM_TYPE_ALIASES(FImpl3, 4); + BASIC_TYPE_ALIASES(ScalarImplCR, Scalar); + SINK_TYPE_ALIASES(Scalar); + typedef typename SpinMatrixField1::vector_object::scalar_object SpinMatrix; + class Metadata: Serializable + { + public: + GRID_SERIALIZABLE_CLASS_MEMBERS(Metadata, + Gamma::Algebra, gamma_H, + Gamma::Algebra, gammaA_sigma, + Gamma::Algebra, gammaB_sigma, + Gamma::Algebra, gammaA_nucl, + Gamma::Algebra, gammaB_nucl, + int, trace, + int, parity); + }; + typedef Correlator Result; +public: + // constructor + TSigmaToNucleonNonEye(const std::string name); + // destructor + virtual ~TSigmaToNucleonNonEye(void) {}; + // dependency relation + virtual std::vector getInput(void); + virtual std::vector getOutput(void); +protected: + // setup + virtual void setup(void); + // execution + virtual void execute(void); + // Which gamma algebra was specified + Gamma::Algebra al; +}; + +MODULE_REGISTER_TMP(SigmaToNucleonNonEye, ARG(TSigmaToNucleonNonEye), MContraction); + +/****************************************************************************** + * TSigmaToNucleonNonEye implementation * + ******************************************************************************/ +// constructor ///////////////////////////////////////////////////////////////// +template +TSigmaToNucleonNonEye::TSigmaToNucleonNonEye(const std::string name) +: Module(name) +{} + +// dependencies/products /////////////////////////////////////////////////////// +template +std::vector TSigmaToNucleonNonEye::getInput(void) +{ + std::vector input = {par().qu_xi, par().qu_xf, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; + + return input; +} + +template +std::vector TSigmaToNucleonNonEye::getOutput(void) +{ + std::vector out = {}; + + return out; +} + +// setup /////////////////////////////////////////////////////////////////////// +template +void TSigmaToNucleonNonEye::setup(void) +{ + envTmpLat(SpinMatrixField1, "c"); +} + +// execution /////////////////////////////////////////////////////////////////// +template +void TSigmaToNucleonNonEye::execute(void) +{ + const int parity {par().parity.size()>0 ? std::stoi(par().parity) : 1}; + const Gamma GammaB(Gamma::Algebra::SigmaXZ); // C*gamma_5 + const Gamma Id(Gamma::Algebra::Identity); // C*gamma_5 + + LOG(Message) << "Computing sigma-to-nucleon contractions '" << getName() << "'" << std::endl; + LOG(Message) << "' with (Gamma^A,Gamma^B)_sigma = ( Identity, C*gamma_5 ) and (Gamma^A,Gamma^B)_nucl = ( Identity, C*gamma_5 )" << std::endl; + LOG(Message) << "and parity " << parity << " using sink " << par().sink << "." << std::endl; + + envGetTmp(SpinMatrixField1, c); + std::vector buf; + + std::vector result; + Result r; + r.info.parity = parity; + r.info.gammaA_sigma = Id.g; + r.info.gammaB_sigma = GammaB.g; + r.info.gammaA_nucl = Id.g; + r.info.gammaB_nucl = GammaB.g; + + auto &qu_xi = envGet(PropagatorField1, par().qu_xi); + auto &qu_xf = envGet(PropagatorField1, par().qu_xf); + auto &qu_spec = envGet(SlicedPropagator2, par().qu_spec); + auto &qd_xf = envGet(PropagatorField3, par().qd_xf); + auto &qs_xi = envGet(PropagatorField4, par().qs_xi); + auto qut = qu_spec[par().tOut]; + for (auto &G: Gamma::gall) + { + r.info.gamma_H = G.g; + //Operator Q1, equivalent to the two-trace case in the rare-kaons module + c=Zero(); + BaryonUtils::Sigma_to_Nucleon(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1_NonEye",c); + sliceSum(c,buf,Tp); + r.corr.clear(); + for (unsigned int t = 0; t < buf.size(); ++t) + { + r.corr.push_back(buf[t]); + } + r.info.trace = 2; + result.push_back(r); + //Operator Q2, equivalent to the one-trace case in the rare-kaons module + c=Zero(); + BaryonUtils::Sigma_to_Nucleon(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2_NonEye",c); + sliceSum(c,buf,Tp); + r.corr.clear(); + for (unsigned int t = 0; t < buf.size(); ++t) + { + r.corr.push_back(buf[t]); + } + r.info.trace = 1; + result.push_back(r); + } + + saveResult(par().output, "StN_NonEye", result); + +} + +END_MODULE_NAMESPACE + +END_HADRONS_NAMESPACE + +#endif // Hadrons_MContraction_SigmaToNucleonNonEye_hpp_ diff --git a/Hadrons/Modules/MContraction/WeakEye3pt.hpp b/Hadrons/Modules/MContraction/WeakEye3pt.hpp index ea7ff529..b4165098 100644 --- a/Hadrons/Modules/MContraction/WeakEye3pt.hpp +++ b/Hadrons/Modules/MContraction/WeakEye3pt.hpp @@ -144,7 +144,7 @@ void TWeakEye3pt::execute(void) { LOG(Message) << "Computing mesonic weak 3pt contractions, eye topologies" << std::endl; LOG(Message) << "gIn : " << par().gammaIn << std::endl; - LOG(Message) << "gOut: " << par().gammaIn << std::endl; + LOG(Message) << "gOut: " << par().gammaOut << std::endl; LOG(Message) << "tOut: " << par().tOut << std::endl; LOG(Message) << "qbl : " << par().qBarLeft << std::endl; LOG(Message) << "qbr : " << par().qBarRight << std::endl; diff --git a/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp b/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp index 08a1188b..4068cd8c 100644 --- a/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp +++ b/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp @@ -144,7 +144,7 @@ void TWeakNonEye3pt::execute(void) { LOG(Message) << "Computing mesonic weak 3pt contractions, non-eye topologies" << std::endl; LOG(Message) << "gIn : " << par().gammaIn << std::endl; - LOG(Message) << "gOut: " << par().gammaIn << std::endl; + LOG(Message) << "gOut: " << par().gammaOut << std::endl; LOG(Message) << "ql : " << par().qLeft << std::endl; LOG(Message) << "qbl : " << par().qBarLeft << std::endl; LOG(Message) << "qr : " << par().qRight << std::endl; diff --git a/Hadrons/Modules/MSource/MomentumPhase.hpp b/Hadrons/Modules/MSource/MomentumPhase.hpp index 89bd0c8e..f2f94af1 100644 --- a/Hadrons/Modules/MSource/MomentumPhase.hpp +++ b/Hadrons/Modules/MSource/MomentumPhase.hpp @@ -140,7 +140,7 @@ void TMomentumPhase::execute(void) envGetTmp(LatticeComplex, coor); p = strToVec(par().mom); - ph = zero; + ph = Zero(); for(unsigned int mu = 0; mu < env().getNd(); mu++) { LatticeCoordinate(coor, mu); diff --git a/Hadrons/modules.inc b/Hadrons/modules.inc index f4ede1bb..bcf58d00 100644 --- a/Hadrons/modules.inc +++ b/Hadrons/modules.inc @@ -1,147 +1,153 @@ modules_cc =\ - Modules/MFermion/FreeProp.cc \ - Modules/MFermion/GaugeProp.cc \ - Modules/MFermion/EMLepton.cc \ - Modules/MIO/LoadA2AVectors.cc \ - Modules/MIO/LoadEigenPack.cc \ - Modules/MIO/LoadCosmHol.cc \ - Modules/MIO/LoadBinary.cc \ - Modules/MIO/LoadA2AMatrixDiskVector.cc \ - Modules/MIO/LoadCoarseEigenPack.cc \ - Modules/MIO/LoadNersc.cc \ - Modules/MAction/ZMobiusDWF.cc \ - Modules/MAction/ScaledDWF.cc \ - Modules/MAction/Wilson.cc \ - Modules/MAction/DWF.cc \ - Modules/MAction/WilsonClover.cc \ - Modules/MAction/MobiusDWF.cc \ Modules/MUtilities/RandomVectors.cc \ Modules/MUtilities/PrecisionCast.cc \ - Modules/MNoise/FullVolumeSpinColorDiagonal.cc \ - Modules/MNoise/TimeDilutedSpinColorDiagonal.cc \ - Modules/MContraction/Gamma3pt.cc \ - Modules/MContraction/A2AFourQuarkContraction.cc \ - Modules/MContraction/A2AAslashField.cc \ - Modules/MContraction/A2ALoop.cc \ - Modules/MContraction/WeakEye3pt.cc \ - Modules/MContraction/WeakNonEye3pt.cc \ - Modules/MContraction/A2AMesonField.cc \ - Modules/MContraction/DiscLoop.cc \ - Modules/MContraction/Baryon.cc \ - Modules/MContraction/WeakMesonDecayKl2.cc \ - Modules/MContraction/Meson.cc \ - Modules/MSink/Point.cc \ - Modules/MSink/Smear.cc \ - Modules/MScalarSUN/TrPhi.cc \ - Modules/MScalarSUN/TrMag.cc \ - Modules/MScalarSUN/TrKinetic.cc \ - Modules/MScalarSUN/TwoPoint.cc \ - Modules/MScalarSUN/Grad.cc \ - Modules/MScalarSUN/TwoPointNPR.cc \ - Modules/MScalarSUN/StochFreeField.cc \ - Modules/MScalarSUN/TransProj.cc \ - Modules/MScalarSUN/EMT.cc \ - Modules/MScalarSUN/Div.cc \ - Modules/MNPR/Amputate.cc \ - Modules/MNPR/Bilinear.cc \ - Modules/MNPR/FourQuark.cc \ - Modules/MSource/SeqGamma.cc \ + Modules/MGauge/FundtoHirep.cc \ + Modules/MGauge/UnitEm.cc \ + Modules/MGauge/StoutSmearing.cc \ + Modules/MGauge/Electrify.cc \ + Modules/MGauge/Unit.cc \ + Modules/MGauge/StochEm.cc \ + Modules/MGauge/GaugeFix.cc \ + Modules/MGauge/Random.cc \ + Modules/MSolver/MixedPrecisionRBPrecCG.cc \ + Modules/MSolver/RBPrecCG.cc \ + Modules/MSolver/LocalCoherenceLanczos.cc \ + Modules/MSolver/A2AAslashVectors.cc \ + Modules/MSolver/A2AVectors.cc \ Modules/MSource/Z2.cc \ - Modules/MSource/Convolution.cc \ - Modules/MSource/Momentum.cc \ Modules/MSource/Wall.cc \ Modules/MSource/Point.cc \ - Modules/MSource/SeqAslash.cc \ + Modules/MSource/MomentumPhase.cc \ Modules/MSource/Gauss.cc \ + Modules/MSource/Convolution.cc \ + Modules/MSource/SeqGamma.cc \ + Modules/MSource/Momentum.cc \ Modules/MSource/SeqConserved.cc \ - Modules/MGauge/StoutSmearing.cc \ - Modules/MGauge/GaugeFix.cc \ - Modules/MGauge/Electrify.cc \ - Modules/MGauge/Random.cc \ - Modules/MGauge/Unit.cc \ - Modules/MGauge/UnitEm.cc \ - Modules/MGauge/FundtoHirep.cc \ - Modules/MGauge/StochEm.cc \ - Modules/MSolver/RBPrecCG.cc \ - Modules/MSolver/A2AAslashVectors.cc \ - Modules/MSolver/MixedPrecisionRBPrecCG.cc \ - Modules/MSolver/A2AVectors.cc \ - Modules/MSolver/LocalCoherenceLanczos.cc \ + Modules/MSource/SeqAslash.cc \ + Modules/MNPR/Bilinear.cc \ + Modules/MNPR/FourQuark.cc \ + Modules/MNPR/Amputate.cc \ + Modules/MContraction/Baryon.cc \ + Modules/MContraction/A2ALoop.cc \ + Modules/MContraction/A2AFourQuarkContraction.cc \ + Modules/MContraction/DiscLoop.cc \ + Modules/MContraction/SigmaToNucleonEye.cc \ + Modules/MContraction/WeakEye3pt.cc \ + Modules/MContraction/A2AAslashField.cc \ + Modules/MContraction/Meson.cc \ + Modules/MContraction/WeakMesonDecayKl2.cc \ + Modules/MContraction/WeakNonEye3pt.cc \ + Modules/MContraction/SigmaToNucleonNonEye.cc \ + Modules/MContraction/Gamma3pt.cc \ + Modules/MContraction/A2AMesonField.cc \ + Modules/MScalar/FreeProp.cc \ Modules/MScalar/ChargedProp.cc \ - Modules/MScalar/FreeProp.cc + Modules/MSink/Point.cc \ + Modules/MSink/Smear.cc \ + Modules/MNoise/TimeDilutedSpinColorDiagonal.cc \ + Modules/MNoise/FullVolumeSpinColorDiagonal.cc \ + Modules/MIO/LoadA2AVectors.cc \ + Modules/MIO/LoadBinary.cc \ + Modules/MIO/LoadA2AMatrixDiskVector.cc \ + Modules/MIO/LoadNersc.cc \ + Modules/MIO/LoadCoarseEigenPack.cc \ + Modules/MIO/LoadCosmHol.cc \ + Modules/MIO/LoadEigenPack.cc \ + Modules/MScalarSUN/Div.cc \ + Modules/MScalarSUN/TwoPointNPR.cc \ + Modules/MScalarSUN/TrPhi.cc \ + Modules/MScalarSUN/TransProj.cc \ + Modules/MScalarSUN/TwoPoint.cc \ + Modules/MScalarSUN/Grad.cc \ + Modules/MScalarSUN/EMT.cc \ + Modules/MScalarSUN/TrKinetic.cc \ + Modules/MScalarSUN/StochFreeField.cc \ + Modules/MScalarSUN/TrMag.cc \ + Modules/MAction/MobiusDWF.cc \ + Modules/MAction/DWF.cc \ + Modules/MAction/ScaledDWF.cc \ + Modules/MAction/Wilson.cc \ + Modules/MAction/ZMobiusDWF.cc \ + Modules/MAction/WilsonClover.cc \ + Modules/MFermion/EMLepton.cc \ + Modules/MFermion/FreeProp.cc \ + Modules/MFermion/GaugeProp.cc modules_hpp =\ - Modules/MFermion/GaugeProp.hpp \ - Modules/MFermion/EMLepton.hpp \ - Modules/MFermion/FreeProp.hpp \ - Modules/MIO/LoadCosmHol.hpp \ - Modules/MIO/LoadEigenPack.hpp \ - Modules/MIO/LoadA2AVectors.hpp \ - Modules/MIO/LoadA2AMatrixDiskVector.hpp \ - Modules/MIO/LoadCoarseEigenPack.hpp \ - Modules/MIO/LoadNersc.hpp \ - Modules/MIO/LoadBinary.hpp \ - Modules/MAction/ZMobiusDWF.hpp \ - Modules/MAction/MobiusDWF.hpp \ - Modules/MAction/Wilson.hpp \ - Modules/MAction/DWF.hpp \ - Modules/MAction/WilsonClover.hpp \ - Modules/MAction/ScaledDWF.hpp \ - Modules/MUtilities/RandomVectors.hpp \ Modules/MUtilities/PrecisionCast.hpp \ - Modules/MNoise/TimeDilutedSpinColorDiagonal.hpp \ - Modules/MNoise/FullVolumeSpinColorDiagonal.hpp \ - Modules/MContraction/DiscLoop.hpp \ - Modules/MContraction/Meson.hpp \ - Modules/MContraction/WeakMesonDecayKl2.hpp \ - Modules/MContraction/A2AMesonField.hpp \ - Modules/MContraction/WeakNonEye3pt.hpp \ - Modules/MContraction/Gamma3pt.hpp \ - Modules/MContraction/A2AAslashField.hpp \ - Modules/MContraction/A2AFourQuarkContraction.hpp \ - Modules/MContraction/Baryon.hpp \ - Modules/MContraction/WeakEye3pt.hpp \ - Modules/MContraction/A2ALoop.hpp \ - Modules/MSink/Point.hpp \ - Modules/MSink/Smear.hpp \ - Modules/MScalarSUN/StochFreeField.hpp \ - Modules/MScalarSUN/EMT.hpp \ - Modules/MScalarSUN/Utils.hpp \ - Modules/MScalarSUN/TwoPoint.hpp \ - Modules/MScalarSUN/TransProj.hpp \ - Modules/MScalarSUN/TwoPointNPR.hpp \ - Modules/MScalarSUN/TrPhi.hpp \ - Modules/MScalarSUN/Grad.hpp \ - Modules/MScalarSUN/TrMag.hpp \ - Modules/MScalarSUN/Div.hpp \ - Modules/MScalarSUN/TrKinetic.hpp \ - Modules/MNPR/Amputate.hpp \ - Modules/MNPR/FourQuark.hpp \ - Modules/MNPR/Bilinear.hpp \ - Modules/MSource/SeqAslash.hpp \ - Modules/MSource/Momentum.hpp \ - Modules/MSource/Z2.hpp \ - Modules/MSource/Point.hpp \ - Modules/MSource/Gauss.hpp \ - Modules/MSource/SeqConserved.hpp \ - Modules/MSource/Wall.hpp \ - Modules/MSource/SeqGamma.hpp \ - Modules/MSource/Convolution.hpp \ - Modules/MGauge/Random.hpp \ + Modules/MUtilities/RandomVectors.hpp \ Modules/MGauge/FundtoHirep.hpp \ - Modules/MGauge/StochEm.hpp \ - Modules/MGauge/UnitEm.hpp \ + Modules/MGauge/Electrify.hpp \ Modules/MGauge/GaugeFix.hpp \ Modules/MGauge/StoutSmearing.hpp \ Modules/MGauge/Unit.hpp \ - Modules/MGauge/Electrify.hpp \ - Modules/MSolver/A2AVectors.hpp \ - Modules/MSolver/RBPrecCG.hpp \ - Modules/MSolver/LocalCoherenceLanczos.hpp \ + Modules/MGauge/UnitEm.hpp \ + Modules/MGauge/Random.hpp \ + Modules/MGauge/StochEm.hpp \ Modules/MSolver/Guesser.hpp \ Modules/MSolver/MixedPrecisionRBPrecCG.hpp \ + Modules/MSolver/A2AVectors.hpp \ + Modules/MSolver/LocalCoherenceLanczos.hpp \ Modules/MSolver/A2AAslashVectors.hpp \ + Modules/MSolver/RBPrecCG.hpp \ + Modules/MSource/Convolution.hpp \ + Modules/MSource/Point.hpp \ + Modules/MSource/Gauss.hpp \ + Modules/MSource/MomentumPhase.hpp \ + Modules/MSource/Z2.hpp \ + Modules/MSource/SeqAslash.hpp \ + Modules/MSource/SeqConserved.hpp \ + Modules/MSource/SeqGamma.hpp \ + Modules/MSource/Momentum.hpp \ + Modules/MSource/Wall.hpp \ + Modules/MNPR/Amputate.hpp \ + Modules/MNPR/Bilinear.hpp \ + Modules/MNPR/FourQuark.hpp \ + Modules/MContraction/Meson.hpp \ + Modules/MContraction/WeakMesonDecayKl2.hpp \ + Modules/MContraction/WeakNonEye3pt.hpp \ + Modules/MContraction/A2AFourQuarkContraction.hpp \ + Modules/MContraction/DiscLoop.hpp \ + Modules/MContraction/WeakEye3pt.hpp \ + Modules/MContraction/A2AAslashField.hpp \ + Modules/MContraction/Gamma3pt.hpp \ + Modules/MContraction/A2AMesonField.hpp \ + Modules/MContraction/SigmaToNucleonEye.hpp \ + Modules/MContraction/A2ALoop.hpp \ + Modules/MContraction/Baryon.hpp \ + Modules/MContraction/SigmaToNucleonNonEye.hpp \ + Modules/MScalar/FreeProp.hpp \ Modules/MScalar/ChargedProp.hpp \ Modules/MScalar/Scalar.hpp \ - Modules/MScalar/FreeProp.hpp + Modules/MSink/Point.hpp \ + Modules/MSink/Smear.hpp \ + Modules/MNoise/FullVolumeSpinColorDiagonal.hpp \ + Modules/MNoise/TimeDilutedSpinColorDiagonal.hpp \ + Modules/MIO/LoadBinary.hpp \ + Modules/MIO/LoadNersc.hpp \ + Modules/MIO/LoadEigenPack.hpp \ + Modules/MIO/LoadA2AMatrixDiskVector.hpp \ + Modules/MIO/LoadCoarseEigenPack.hpp \ + Modules/MIO/LoadA2AVectors.hpp \ + Modules/MIO/LoadCosmHol.hpp \ + Modules/MScalarSUN/TwoPointNPR.hpp \ + Modules/MScalarSUN/TransProj.hpp \ + Modules/MScalarSUN/TwoPoint.hpp \ + Modules/MScalarSUN/TrPhi.hpp \ + Modules/MScalarSUN/Utils.hpp \ + Modules/MScalarSUN/TrMag.hpp \ + Modules/MScalarSUN/EMT.hpp \ + Modules/MScalarSUN/Grad.hpp \ + Modules/MScalarSUN/StochFreeField.hpp \ + Modules/MScalarSUN/Div.hpp \ + Modules/MScalarSUN/TrKinetic.hpp \ + Modules/MAction/ZMobiusDWF.hpp \ + Modules/MAction/ScaledDWF.hpp \ + Modules/MAction/WilsonClover.hpp \ + Modules/MAction/Wilson.hpp \ + Modules/MAction/MobiusDWF.hpp \ + Modules/MAction/DWF.hpp \ + Modules/MFermion/FreeProp.hpp \ + Modules/MFermion/EMLepton.hpp \ + Modules/MFermion/GaugeProp.hpp diff --git a/tests/hadrons/Test_sigma_to_nucleon.cc b/tests/hadrons/Test_sigma_to_nucleon.cc new file mode 100644 index 00000000..621a8f7c --- /dev/null +++ b/tests/hadrons/Test_sigma_to_nucleon.cc @@ -0,0 +1,154 @@ +/************************************************************************************* + +Grid physics library, www.github.com/paboyle/Grid + +Source file: Tests/Hadrons/Test_hadrons_spectrum.cc + +Copyright (C) 2015-2018 + + Author: Antonin Portelli + +This program is free software; you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation; either version 2 of the License, or +(at your option) any later version. + +This program is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + +You should have received a copy of the GNU General Public License along +with this program; if not, write to the Free Software Foundation, Inc., +51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. + +See the full license in the file "LICENSE" in the top level distribution directory +*************************************************************************************/ +/* END LEGAL */ + +#include +#include + +using namespace Grid; +using namespace Hadrons; + +int main(int argc, char *argv[]) +{ + // initialization ////////////////////////////////////////////////////////// + Grid_init(&argc, &argv); + HadronsLogError.Active(GridLogError.isActive()); + HadronsLogWarning.Active(GridLogWarning.isActive()); + HadronsLogMessage.Active(GridLogMessage.isActive()); + HadronsLogIterative.Active(GridLogIterative.isActive()); + HadronsLogDebug.Active(GridLogDebug.isActive()); + LOG(Message) << "Grid initialized" << std::endl; + + // run setup /////////////////////////////////////////////////////////////// + Application application; + std::vector flavour = {"l", "s", "c"}; + std::vector mass = {.01, .04, .2 }; + + // global parameters + Application::GlobalPar globalPar; + globalPar.trajCounter.start = 1500; + globalPar.trajCounter.end = 1520; + globalPar.trajCounter.step = 20; + globalPar.runId = "test"; + application.setPar(globalPar); + // gauge field + application.createModule("gauge"); + // sources + MSource::Point::Par ptPar; + ptPar.position = "0 0 0 0"; + application.createModule("pt_0", ptPar); + ptPar.position = "0 0 0 4"; + application.createModule("pt_4", ptPar); + // sink + MSink::Point::Par sinkPar; + sinkPar.mom = "0 0 0"; + application.createModule("sink", sinkPar); + application.createModule("sink_spec", sinkPar); + + // set fermion boundary conditions to be periodic space, antiperiodic time. + std::string boundary = "1 1 1 -1"; + std::string twist = "0. 0. 0. 0."; + + for (unsigned int i = 0; i < flavour.size(); ++i) + { + // actions + MAction::DWF::Par actionPar; + actionPar.gauge = "gauge"; + actionPar.Ls = 12; + actionPar.M5 = 1.8; + actionPar.mass = mass[i]; + actionPar.boundary = boundary; + actionPar.twist = twist; + application.createModule("DWF_" + flavour[i], actionPar); + + // solvers + MSolver::RBPrecCG::Par solverPar; + solverPar.action = "DWF_" + flavour[i]; + solverPar.residual = 1.0e-8; + solverPar.maxIteration = 10000; + application.createModule("CG_" + flavour[i], + solverPar); + + } + + // propagators + MFermion::GaugeProp::Par quarkPar; + quarkPar.solver = "CG_l"; + quarkPar.source = "pt_0"; + application.createModule("Qpt_l_0", quarkPar); + quarkPar.source = "pt_4"; + application.createModule("Qpt_l_4", quarkPar); + quarkPar.solver = "CG_s"; + quarkPar.source = "pt_0"; + application.createModule("Qpt_s_0", quarkPar); + //This should be a loop - how do I make this? + quarkPar.solver = "CG_c"; + quarkPar.source = "pt_0"; + application.createModule("Qpt_c_loop", quarkPar); + quarkPar.solver = "CG_l"; + quarkPar.source = "pt_0"; + application.createModule("Qpt_l_loop", quarkPar); + + MSink::Smear::Par smearPar; + smearPar.q="Qpt_l_0"; + smearPar.sink = "sink_spec"; + application.createModule("Qpt_u_spec",smearPar); + + + MContraction::SigmaToNucleonEye::Par EyePar; + EyePar.output = "SigmaToNucleon/Eye_u"; + EyePar.qq_loop = "Qpt_l_loop"; + EyePar.qu_spec = "Qpt_u_spec"; + EyePar.qd_xf = "Qpt_l_4"; + EyePar.qs_xi = "Qpt_s_0"; + EyePar.tOut = 0; + EyePar.sink = "sink"; + application.createModule("SigmaToNucleonEye_u", EyePar); + EyePar.output = "SigmaToNucleon/Eye_c"; + EyePar.qq_loop = "Qpt_c_loop"; + application.createModule("SigmaToNucleonEye_c", EyePar); + MContraction::SigmaToNucleonNonEye::Par NonEyePar; + NonEyePar.output = "SigmaToNucleon/NonEye"; + NonEyePar.qu_xi = "Qpt_l_0"; + NonEyePar.qu_xf = "Qpt_l_4"; + NonEyePar.qu_spec = "Qpt_u_spec"; + NonEyePar.qd_xf = "Qpt_l_4"; + NonEyePar.qs_xi = "Qpt_s_0"; + NonEyePar.tOut = 0; + NonEyePar.sink = "sink"; + application.createModule("SigmaToNucleonNonEye", NonEyePar); + + // execution + application.saveParameterFile("stn.xml"); + application.run(); + + // epilogue + LOG(Message) << "Grid is finalizing now" << std::endl; + Grid_finalize(); + + return EXIT_SUCCESS; +} From 5fd5c25114b7fc6e66a2df8c52b1e2daa80ba18e Mon Sep 17 00:00:00 2001 From: Felix Erben Date: Tue, 26 Nov 2019 13:44:55 +0000 Subject: [PATCH 02/14] now two seperate functions for Eye and NonEye --- Grid/qcd/utils/BaryonUtils.h | 61 ++++++++++++++++--- .../MContraction/SigmaToNucleonEye.hpp | 4 +- .../MContraction/SigmaToNucleonNonEye.hpp | 4 +- 3 files changed, 57 insertions(+), 12 deletions(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 7fcc4502..580f91ef 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -131,7 +131,18 @@ public: robj &result); public: template - static void Sigma_to_Nucleon(const PropagatorField &qq_xi, + static void Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, + const mobj &Du_spec, + const PropagatorField &qd_xf, + const PropagatorField &qs_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + const std::string op, + SpinMatrixField &stn_corr); + template + static void Sigma_to_Nucleon_NonEye(const PropagatorField &qq_xi, const PropagatorField &qq_xf, const mobj &Du_spec, const PropagatorField &qd_xf, @@ -505,10 +516,48 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, } } +template +template +void BaryonUtils::Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, + const mobj &Du_spec, + const PropagatorField &qd_xf, + const PropagatorField &qs_xi, + const Gamma Gamma_H, + const Gamma GammaB_sigma, + const Gamma GammaB_nucl, + const int parity, + const std::string op, + SpinMatrixField &stn_corr) +{ + assert(parity==1 || parity == -1 && "Parity must be +1 or -1"); + + GridBase *grid = qs_xi.Grid(); + + auto vcorr= stn_corr.View(); + auto vq_loop = qq_loop.View(); + auto vd_xf = qd_xf.View(); + auto vs_xi = qs_xi.View(); + + // accelerator_for(ss, grid->oSites(), grid->Nsimd(), { + thread_for(ss,grid->oSites(),{ + auto Dq_loop = vq_loop[ss]; + auto Dd_xf = vd_xf[ss]; + auto Ds_xi = vs_xi[ss]; + sobj result=Zero(); + if(op == "Q1"){ + Sigma_to_Nucleon_Q1_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + } else if(op == "Q2"){ + Sigma_to_Nucleon_Q2_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + } else { + assert(0 && "Weak Operator not correctly specified"); + } + vcorr[ss] = result; + } );//end loop over lattice sites +} template template -void BaryonUtils::Sigma_to_Nucleon(const PropagatorField &qq_xi, +void BaryonUtils::Sigma_to_Nucleon_NonEye(const PropagatorField &qq_xi, const PropagatorField &qq_xf, const mobj &Du_spec, const PropagatorField &qd_xf, @@ -537,14 +586,10 @@ void BaryonUtils::Sigma_to_Nucleon(const PropagatorField &qq_xi, auto Dd_xf = vd_xf[ss]; auto Ds_xi = vs_xi[ss]; sobj result=Zero(); - if(op == "Q1_NonEye"){ + if(op == "Q1"){ Sigma_to_Nucleon_Q1_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); - } else if(op == "Q1_Eye"){ - Sigma_to_Nucleon_Q1_Eye_site(Dq_xi,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); - } else if(op == "Q2_NonEye"){ + } else if(op == "Q2"){ Sigma_to_Nucleon_Q2_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); - } else if(op == "Q2_Eye"){ - Sigma_to_Nucleon_Q2_Eye_site(Dq_xi,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); } else { assert(0 && "Weak Operator not correctly specified"); } diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index e4189f7c..99f74fcd 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -186,7 +186,7 @@ void TSigmaToNucleonEye::execute(void) r.info.gamma_H = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon(qq_loop,qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1_Eye",c); + BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) @@ -197,7 +197,7 @@ void TSigmaToNucleonEye::execute(void) result.push_back(r); //Operator Q2, equivalent to the one-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon(qq_loop,qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2_Eye",c); + BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index 15dc5eaa..c74a5f58 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -192,7 +192,7 @@ void TSigmaToNucleonNonEye::execute(void) r.info.gamma_H = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1_NonEye",c); + BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) @@ -203,7 +203,7 @@ void TSigmaToNucleonNonEye::execute(void) result.push_back(r); //Operator Q2, equivalent to the one-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2_NonEye",c); + BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) From 799ff0c96e8f21733dcdb679606a7527cd9b318d Mon Sep 17 00:00:00 2001 From: Felix Erben Date: Tue, 26 Nov 2019 15:28:47 +0000 Subject: [PATCH 03/14] speed-up --- Grid/qcd/utils/BaryonUtils.h | 90 +++++++++++++++++++++++------------- 1 file changed, 58 insertions(+), 32 deletions(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 580f91ef..5440e0e5 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -361,15 +361,20 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, int a_s = epsilon[ie_s][0]; //a' int b_s = epsilon[ie_s][1]; //b' int c_s = epsilon[ie_s][2]; //c' - for (int j=0; j(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DqG()(tau2,tau2)(j,j) * DuG()(alpha_s, beta_n)(a_s,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); - result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DqG()(tau2,tau2)(j,j) * DuG()(gamma_s, beta_n)(c_s,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); - }}}}}} + auto GDsGDd_ab_bb = GDsGDd()(alpha_s,beta_n)(b_s,b_n); + for (int tau2=0; tau2(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd_ab_bb * DqG_tt_jj; + for (int gamma_s=0; gamma_s::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, int a_s = epsilon[ie_s][0]; //a' int b_s = epsilon[ie_s][1]; //b' int c_s = epsilon[ie_s][2]; //c' - for (int j=0; j(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(alpha_s,tau2)(a_s,j) * DuG()(gamma_s, beta_n)(c_s,a_n) * adjDu()(tau2,gamma_n)(j,c_n); - result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(gamma_s,tau2)(c_s,j) * adjDuG()(tau2, beta_n)(j,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); - result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(alpha_s,tau2)(a_s,j) * adjDuG()(tau2, beta_n)(j,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); - result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd()(alpha_s,beta_n)(b_s,b_n) * DuGH()(gamma_s,tau2)(c_s,j) * DuG()(alpha_s, beta_n)(a_s,a_n) * adjDu()(tau2,gamma_n)(j,c_n); - }}}}}} + auto GDsGDd_ab_bb = GDsGDd()(alpha_s,beta_n)(b_s,b_n); + for (int tau2=0; tau2(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd_ab_bb * DuGH_at_aj; + for (int gamma_s=0; gamma_s(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd_ab_bb * DuGH_gt_cj; + for (int gamma_n=0; gamma_n::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, int a_s = epsilon[ie_s][0]; //a' int b_s = epsilon[ie_s][1]; //b' int c_s = epsilon[ie_s][2]; //c' - for (int i=0; i(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DqGDd()(tau,beta_n)(i,b_n) * DuG()(alpha_s, beta_n)(a_s,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); - result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DqGDd()(tau,beta_n)(i,b_n) * DuG()(gamma_s, beta_n)(c_s,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); - }}}}}} + for (int i=0; i(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG_at_bi * DqGDd_tb_ib; + for (int gamma_s=0; gamma_s::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, int a_s = epsilon[ie_s][0]; //a' int b_s = epsilon[ie_s][1]; //b' int c_s = epsilon[ie_s][2]; //c' - for (int i=0; i(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(alpha_s,beta_n)(a_s,b_n) * DuG()(gamma_s, beta_n)(c_s,a_n) * adjDu()(tau,gamma_n)(i,c_n); - result()(gamma_s,gamma_n)() -= static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(gamma_s,beta_n)(c_s,b_n) * adjDuG()(tau, beta_n)(i,a_n) * Du_spec()(alpha_s,gamma_n)(a_s,c_n); - result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(alpha_s,beta_n)(a_s,b_n) * adjDuG()(tau, beta_n)(i,a_n) * Du_spec()(gamma_s,gamma_n)(c_s,c_n); - result()(gamma_s,gamma_n)() += static_cast(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG()(alpha_s,tau)(b_s,i) * DuGDd()(gamma_s,beta_n)(c_s,b_n) * DuG()(alpha_s, beta_n)(a_s,a_n) * adjDu()(tau,gamma_n)(i,c_n); - }}}}}} + for (int i=0; i(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG_at_bi * DuGDd_ab_ab; + for (int gamma_s=0; gamma_s(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsG_at_bi * DuGDd_gb_cb; + for (int gamma_n=0; gamma_n Date: Thu, 28 Nov 2019 12:35:18 +0000 Subject: [PATCH 04/14] cleared up how exactly q_spec has to be defined --- Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp | 10 ++++++++-- Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp | 10 ++++++++-- 2 files changed, 16 insertions(+), 4 deletions(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index 99f74fcd..f9de7dae 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -60,6 +60,12 @@ BEGIN_HADRONS_NAMESPACE * * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and * the diagram on the right 'two-trace' + * + * Propagators: + * * qq_loop + * * qu_spec, source at xi + * * qd_xf, source at xf + * * qs_xi, source at xi */ BEGIN_MODULE_NAMESPACE(MContraction) @@ -71,7 +77,7 @@ public: std::string, qu_spec, std::string, qd_xf, std::string, qs_xi, - unsigned int, tOut, + unsigned int, xf, std::string, parity, std::string, sink, std::string, output); @@ -180,7 +186,7 @@ void TSigmaToNucleonEye::execute(void) auto &qu_spec = envGet(SlicedPropagator2, par().qu_spec); auto &qd_xf = envGet(PropagatorField3, par().qd_xf); auto &qs_xi = envGet(PropagatorField4, par().qs_xi); - auto qut = qu_spec[par().tOut]; + auto qut = qu_spec[par().xf]; for (auto &G: Gamma::gall) { r.info.gamma_H = G.g; diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index c74a5f58..43891c68 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -64,6 +64,12 @@ BEGIN_HADRONS_NAMESPACE * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and * the diagram on the right 'two-trace' * + * Propagators: + * * qu_xi, source at xi + * * qu_xf, source at xf + * * qu_spec, source at xi + * * qd_xf, source at xf + * * qs_xi, source at xi */ BEGIN_MODULE_NAMESPACE(MContraction) @@ -76,7 +82,7 @@ public: std::string, qu_spec, std::string, qd_xf, std::string, qs_xi, - unsigned int, tOut, + unsigned int, xf, std::string, parity, std::string, sink, std::string, output); @@ -186,7 +192,7 @@ void TSigmaToNucleonNonEye::execute(void) auto &qu_spec = envGet(SlicedPropagator2, par().qu_spec); auto &qd_xf = envGet(PropagatorField3, par().qd_xf); auto &qs_xi = envGet(PropagatorField4, par().qs_xi); - auto qut = qu_spec[par().tOut]; + auto qut = qu_spec[par().xf]; for (auto &G: Gamma::gall) { r.info.gamma_H = G.g; From 07feaf953141e5625218ee61446e717f240e8db6 Mon Sep 17 00:00:00 2001 From: ferben Date: Tue, 3 Dec 2019 11:17:35 +0000 Subject: [PATCH 05/14] updated ascii-doc preamble --- Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp | 8 ++++---- Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp | 6 +++--- 2 files changed, 7 insertions(+), 7 deletions(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index f9de7dae..3146c99c 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -44,13 +44,13 @@ BEGIN_HADRONS_NAMESPACE * Sigma-to-nucleon 3-pt diagrams, eye topologies. * * Schematics: qq_loop | - * /-<-¬ | + * /->-¬ | * / \ | qs_xi G qd_xf - * \ / | /----<------*------<----¬ + * \ / | /---->------*------>----¬ * qs_xi \ / qd_xf | / /-*-¬ \ - * /-----<-----* *-----<----¬ | / / G \ \ + * /----->-----* *----->----¬ | / / G \ \ * * G G * | * \ / qq_loop * - * |\ /| | |\ \->-/ /| + * |\ /| | |\ \-<-/ /| * | \ / | | | \ / | * | \---------->---------/ | | | \----------->----------/ | * \ qu_spec / | \ qu_spec / diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index 43891c68..370ee603 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -44,8 +44,8 @@ BEGIN_HADRONS_NAMESPACE * Sigma-to-Nucleon 3-pt diagrams, non-eye topologies. * * Schematic: - * qs_xi qd_xf | qs_xi qd_xf - * /--<--¬ /--<--¬ | /--<--¬ /--<--¬ + * qs_xi qu_xf | qs_xi qd_xf + * /-->--¬ /-->--¬ | /-->--¬ /-->--¬ * / \ / \ | / \ / \ * / \ / \ | / \ / \ * / \ / \ | / \ / \ @@ -55,7 +55,7 @@ BEGIN_HADRONS_NAMESPACE * | \ / \ / | | | \ / \ / | * | \ / \ / | | | \-->--/ \-->--/ | * \ \-->--/ \-->--/ / | \ qu_xi qu_xf / - * \ qu_xf qu_xf / | \ / + * \ qu_xi qd_xf / | \ / * \ / | \ / * \--------->----------/ | \--------->-----------/ * qu_spec | qu_spec From 9c7f269489f21d4503dab6a262510e82a8873476 Mon Sep 17 00:00:00 2001 From: ferben Date: Tue, 3 Dec 2019 11:19:54 +0000 Subject: [PATCH 06/14] typo in fimpl4 --- Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp | 2 +- Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index 3146c99c..75393b3a 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -90,7 +90,7 @@ public: FERM_TYPE_ALIASES(FImpl1, 1); FERM_TYPE_ALIASES(FImpl2, 2); FERM_TYPE_ALIASES(FImpl3, 3); - FERM_TYPE_ALIASES(FImpl3, 4); + FERM_TYPE_ALIASES(FImpl4, 4); BASIC_TYPE_ALIASES(ScalarImplCR, Scalar); SINK_TYPE_ALIASES(Scalar); typedef typename SpinMatrixField1::vector_object::scalar_object SpinMatrix; diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index 370ee603..5fb62259 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -95,7 +95,7 @@ public: FERM_TYPE_ALIASES(FImpl1, 1); FERM_TYPE_ALIASES(FImpl2, 2); FERM_TYPE_ALIASES(FImpl3, 3); - FERM_TYPE_ALIASES(FImpl3, 4); + FERM_TYPE_ALIASES(FImpl4, 4); BASIC_TYPE_ALIASES(ScalarImplCR, Scalar); SINK_TYPE_ALIASES(Scalar); typedef typename SpinMatrixField1::vector_object::scalar_object SpinMatrix; From e940f4db7ea2f3adca2f4c5500bb1dfbf0135b80 Mon Sep 17 00:00:00 2001 From: ferben Date: Tue, 3 Dec 2019 12:01:31 +0000 Subject: [PATCH 07/14] removed unused parameter parity --- Grid/qcd/utils/BaryonUtils.h | 24 ++++--------------- .../MContraction/SigmaToNucleonEye.hpp | 12 ++++------ .../MContraction/SigmaToNucleonNonEye.hpp | 12 ++++------ tests/hadrons/Test_sigma_to_nucleon.cc | 4 ++-- 4 files changed, 14 insertions(+), 38 deletions(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 5440e0e5..7c0e91a4 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -93,7 +93,6 @@ public: const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result); template static void Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, @@ -104,7 +103,6 @@ public: const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result); @@ -116,7 +114,6 @@ public: const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result); template static void Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, @@ -127,7 +124,6 @@ public: const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result); public: template @@ -138,7 +134,6 @@ public: const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, const std::string op, SpinMatrixField &stn_corr); template @@ -150,7 +145,6 @@ public: const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, const std::string op, SpinMatrixField &stn_corr); }; @@ -341,7 +335,6 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result) { @@ -394,7 +387,6 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result) { @@ -449,7 +441,6 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result) { @@ -497,7 +488,6 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, robj &result) { @@ -551,12 +541,9 @@ void BaryonUtils::Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, const std::string op, SpinMatrixField &stn_corr) { - assert(parity==1 || parity == -1 && "Parity must be +1 or -1"); - GridBase *grid = qs_xi.Grid(); auto vcorr= stn_corr.View(); @@ -571,9 +558,9 @@ void BaryonUtils::Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, auto Ds_xi = vs_xi[ss]; sobj result=Zero(); if(op == "Q1"){ - Sigma_to_Nucleon_Q1_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + Sigma_to_Nucleon_Q1_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else if(op == "Q2"){ - Sigma_to_Nucleon_Q2_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + Sigma_to_Nucleon_Q2_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else { assert(0 && "Weak Operator not correctly specified"); } @@ -591,12 +578,9 @@ void BaryonUtils::Sigma_to_Nucleon_NonEye(const PropagatorField &qq_xi, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, - const int parity, const std::string op, SpinMatrixField &stn_corr) { - assert(parity==1 || parity == -1 && "Parity must be +1 or -1"); - GridBase *grid = qs_xi.Grid(); auto vcorr= stn_corr.View(); @@ -613,9 +597,9 @@ void BaryonUtils::Sigma_to_Nucleon_NonEye(const PropagatorField &qq_xi, auto Ds_xi = vs_xi[ss]; sobj result=Zero(); if(op == "Q1"){ - Sigma_to_Nucleon_Q1_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + Sigma_to_Nucleon_Q1_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else if(op == "Q2"){ - Sigma_to_Nucleon_Q2_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,parity,result); + Sigma_to_Nucleon_Q2_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else { assert(0 && "Weak Operator not correctly specified"); } diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index 75393b3a..781c9d0a 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -78,7 +78,6 @@ public: std::string, qd_xf, std::string, qs_xi, unsigned int, xf, - std::string, parity, std::string, sink, std::string, output); }; @@ -103,8 +102,7 @@ public: Gamma::Algebra, gammaB_sigma, Gamma::Algebra, gammaA_nucl, Gamma::Algebra, gammaB_nucl, - int, trace, - int, parity); + int, trace); }; typedef Correlator Result; public: @@ -163,20 +161,18 @@ void TSigmaToNucleonEye::setup(void) template void TSigmaToNucleonEye::execute(void) { - const int parity {par().parity.size()>0 ? std::stoi(par().parity) : 1}; const Gamma GammaB(Gamma::Algebra::SigmaXZ); // C*gamma_5 const Gamma Id(Gamma::Algebra::Identity); // C*gamma_5 LOG(Message) << "Computing sigma-to-nucleon contractions '" << getName() << "'" << std::endl; LOG(Message) << "' with (Gamma^A,Gamma^B)_sigma = ( Identity, C*gamma_5 ) and (Gamma^A,Gamma^B)_nucl = ( Identity, C*gamma_5 )" << std::endl; - LOG(Message) << "and parity " << parity << " using sink " << par().sink << "." << std::endl; + LOG(Message) << " using sink " << par().sink << "." << std::endl; envGetTmp(SpinMatrixField1, c); std::vector buf; std::vector result; Result r; - r.info.parity = parity; r.info.gammaA_sigma = Id.g; r.info.gammaB_sigma = GammaB.g; r.info.gammaA_nucl = Id.g; @@ -192,7 +188,7 @@ void TSigmaToNucleonEye::execute(void) r.info.gamma_H = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1",c); + BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q1",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) @@ -203,7 +199,7 @@ void TSigmaToNucleonEye::execute(void) result.push_back(r); //Operator Q2, equivalent to the one-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2",c); + BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q2",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index 5fb62259..58c25e90 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -83,7 +83,6 @@ public: std::string, qd_xf, std::string, qs_xi, unsigned int, xf, - std::string, parity, std::string, sink, std::string, output); }; @@ -108,8 +107,7 @@ public: Gamma::Algebra, gammaB_sigma, Gamma::Algebra, gammaA_nucl, Gamma::Algebra, gammaB_nucl, - int, trace, - int, parity); + int, trace); }; typedef Correlator Result; public: @@ -168,20 +166,18 @@ void TSigmaToNucleonNonEye::setup(void) template void TSigmaToNucleonNonEye::execute(void) { - const int parity {par().parity.size()>0 ? std::stoi(par().parity) : 1}; const Gamma GammaB(Gamma::Algebra::SigmaXZ); // C*gamma_5 const Gamma Id(Gamma::Algebra::Identity); // C*gamma_5 LOG(Message) << "Computing sigma-to-nucleon contractions '" << getName() << "'" << std::endl; LOG(Message) << "' with (Gamma^A,Gamma^B)_sigma = ( Identity, C*gamma_5 ) and (Gamma^A,Gamma^B)_nucl = ( Identity, C*gamma_5 )" << std::endl; - LOG(Message) << "and parity " << parity << " using sink " << par().sink << "." << std::endl; + LOG(Message) << " using sink " << par().sink << "." << std::endl; envGetTmp(SpinMatrixField1, c); std::vector buf; std::vector result; Result r; - r.info.parity = parity; r.info.gammaA_sigma = Id.g; r.info.gammaB_sigma = GammaB.g; r.info.gammaA_nucl = Id.g; @@ -198,7 +194,7 @@ void TSigmaToNucleonNonEye::execute(void) r.info.gamma_H = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q1",c); + BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q1",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) @@ -209,7 +205,7 @@ void TSigmaToNucleonNonEye::execute(void) result.push_back(r); //Operator Q2, equivalent to the one-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,parity,"Q2",c); + BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q2",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) diff --git a/tests/hadrons/Test_sigma_to_nucleon.cc b/tests/hadrons/Test_sigma_to_nucleon.cc index 621a8f7c..0489438c 100644 --- a/tests/hadrons/Test_sigma_to_nucleon.cc +++ b/tests/hadrons/Test_sigma_to_nucleon.cc @@ -125,7 +125,7 @@ int main(int argc, char *argv[]) EyePar.qu_spec = "Qpt_u_spec"; EyePar.qd_xf = "Qpt_l_4"; EyePar.qs_xi = "Qpt_s_0"; - EyePar.tOut = 0; + EyePar.xf = 4; EyePar.sink = "sink"; application.createModule("SigmaToNucleonEye_u", EyePar); EyePar.output = "SigmaToNucleon/Eye_c"; @@ -138,7 +138,7 @@ int main(int argc, char *argv[]) NonEyePar.qu_spec = "Qpt_u_spec"; NonEyePar.qd_xf = "Qpt_l_4"; NonEyePar.qs_xi = "Qpt_s_0"; - NonEyePar.tOut = 0; + NonEyePar.xf = 4; NonEyePar.sink = "sink"; application.createModule("SigmaToNucleonNonEye", NonEyePar); From b766038810fba9ccb1b7b6d688eb8e67cddca2b9 Mon Sep 17 00:00:00 2001 From: ferben Date: Wed, 4 Dec 2019 18:08:00 +0000 Subject: [PATCH 08/14] new syntax after merge --- Grid/qcd/utils/BaryonUtils.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index f3d5787c..8e48cf4e 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -365,7 +365,7 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, for (int tau2=0; tau2(epsilon_sgn[ie_n] * epsilon_sgn[ie_s]) * GDsGDd_ab_bb * DqG_tt_jj; + auto ee_GDGDDG = epsilon_sgn[ie_n] * epsilon_sgn[ie_s] * GDsGDd_ab_bb * DqG_tt_jj; for (int gamma_s=0; gamma_s Date: Thu, 5 Dec 2019 16:46:58 +0000 Subject: [PATCH 09/14] only one FIMPL left! --- .../Modules/MContraction/SigmaToNucleonEye.cc | 2 +- .../MContraction/SigmaToNucleonEye.hpp | 43 +++++++++--------- .../MContraction/SigmaToNucleonNonEye.cc | 2 +- .../MContraction/SigmaToNucleonNonEye.hpp | 45 +++++++++---------- 4 files changed, 43 insertions(+), 49 deletions(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc b/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc index 4e938309..ca4bb341 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.cc @@ -4,4 +4,4 @@ using namespace Grid; using namespace Hadrons; using namespace MContraction; -template class Grid::Hadrons::MContraction::TSigmaToNucleonEye; +template class Grid::Hadrons::MContraction::TSigmaToNucleonEye; diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index 781c9d0a..1aa0b9b7 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -82,17 +82,14 @@ public: std::string, output); }; -template +template class TSigmaToNucleonEye: public Module { public: - FERM_TYPE_ALIASES(FImpl1, 1); - FERM_TYPE_ALIASES(FImpl2, 2); - FERM_TYPE_ALIASES(FImpl3, 3); - FERM_TYPE_ALIASES(FImpl4, 4); + FERM_TYPE_ALIASES(FImpl,); BASIC_TYPE_ALIASES(ScalarImplCR, Scalar); SINK_TYPE_ALIASES(Scalar); - typedef typename SpinMatrixField1::vector_object::scalar_object SpinMatrix; + typedef typename SpinMatrixField::vector_object::scalar_object SpinMatrix; class Metadata: Serializable { public: @@ -122,28 +119,28 @@ protected: Gamma::Algebra al; }; -MODULE_REGISTER_TMP(SigmaToNucleonEye, ARG(TSigmaToNucleonEye), MContraction); +MODULE_REGISTER_TMP(SigmaToNucleonEye, ARG(TSigmaToNucleonEye), MContraction); /****************************************************************************** * TSigmaToNucleonEye implementation * ******************************************************************************/ // constructor ///////////////////////////////////////////////////////////////// -template -TSigmaToNucleonEye::TSigmaToNucleonEye(const std::string name) +template +TSigmaToNucleonEye::TSigmaToNucleonEye(const std::string name) : Module(name) {} // dependencies/products /////////////////////////////////////////////////////// -template -std::vector TSigmaToNucleonEye::getInput(void) +template +std::vector TSigmaToNucleonEye::getInput(void) { std::vector input = {par().qq_loop, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; return input; } -template -std::vector TSigmaToNucleonEye::getOutput(void) +template +std::vector TSigmaToNucleonEye::getOutput(void) { std::vector out = {}; @@ -151,15 +148,15 @@ std::vector TSigmaToNucleonEye::get } // setup /////////////////////////////////////////////////////////////////////// -template -void TSigmaToNucleonEye::setup(void) +template +void TSigmaToNucleonEye::setup(void) { - envTmpLat(SpinMatrixField1, "c"); + envTmpLat(SpinMatrixField, "c"); } // execution /////////////////////////////////////////////////////////////////// -template -void TSigmaToNucleonEye::execute(void) +template +void TSigmaToNucleonEye::execute(void) { const Gamma GammaB(Gamma::Algebra::SigmaXZ); // C*gamma_5 const Gamma Id(Gamma::Algebra::Identity); // C*gamma_5 @@ -168,7 +165,7 @@ void TSigmaToNucleonEye::execute(void) LOG(Message) << "' with (Gamma^A,Gamma^B)_sigma = ( Identity, C*gamma_5 ) and (Gamma^A,Gamma^B)_nucl = ( Identity, C*gamma_5 )" << std::endl; LOG(Message) << " using sink " << par().sink << "." << std::endl; - envGetTmp(SpinMatrixField1, c); + envGetTmp(SpinMatrixField, c); std::vector buf; std::vector result; @@ -178,10 +175,10 @@ void TSigmaToNucleonEye::execute(void) r.info.gammaA_nucl = Id.g; r.info.gammaB_nucl = GammaB.g; - auto &qq_loop = envGet(PropagatorField1, par().qq_loop); - auto &qu_spec = envGet(SlicedPropagator2, par().qu_spec); - auto &qd_xf = envGet(PropagatorField3, par().qd_xf); - auto &qs_xi = envGet(PropagatorField4, par().qs_xi); + auto &qq_loop = envGet(PropagatorField, par().qq_loop); + auto &qu_spec = envGet(SlicedPropagator, par().qu_spec); + auto &qd_xf = envGet(PropagatorField, par().qd_xf); + auto &qs_xi = envGet(PropagatorField, par().qs_xi); auto qut = qu_spec[par().xf]; for (auto &G: Gamma::gall) { diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc index 7d5ee247..117541a3 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.cc @@ -4,4 +4,4 @@ using namespace Grid; using namespace Hadrons; using namespace MContraction; -template class Grid::Hadrons::MContraction::TSigmaToNucleonNonEye; +template class Grid::Hadrons::MContraction::TSigmaToNucleonNonEye; diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index 58c25e90..cb73af88 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -87,17 +87,14 @@ public: std::string, output); }; -template +template class TSigmaToNucleonNonEye: public Module { public: - FERM_TYPE_ALIASES(FImpl1, 1); - FERM_TYPE_ALIASES(FImpl2, 2); - FERM_TYPE_ALIASES(FImpl3, 3); - FERM_TYPE_ALIASES(FImpl4, 4); + FERM_TYPE_ALIASES(FImpl,); BASIC_TYPE_ALIASES(ScalarImplCR, Scalar); SINK_TYPE_ALIASES(Scalar); - typedef typename SpinMatrixField1::vector_object::scalar_object SpinMatrix; + typedef typename SpinMatrixField::vector_object::scalar_object SpinMatrix; class Metadata: Serializable { public: @@ -127,28 +124,28 @@ protected: Gamma::Algebra al; }; -MODULE_REGISTER_TMP(SigmaToNucleonNonEye, ARG(TSigmaToNucleonNonEye), MContraction); +MODULE_REGISTER_TMP(SigmaToNucleonNonEye, ARG(TSigmaToNucleonNonEye), MContraction); /****************************************************************************** * TSigmaToNucleonNonEye implementation * ******************************************************************************/ // constructor ///////////////////////////////////////////////////////////////// -template -TSigmaToNucleonNonEye::TSigmaToNucleonNonEye(const std::string name) +template +TSigmaToNucleonNonEye::TSigmaToNucleonNonEye(const std::string name) : Module(name) {} // dependencies/products /////////////////////////////////////////////////////// -template -std::vector TSigmaToNucleonNonEye::getInput(void) +template +std::vector TSigmaToNucleonNonEye::getInput(void) { std::vector input = {par().qu_xi, par().qu_xf, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; return input; } -template -std::vector TSigmaToNucleonNonEye::getOutput(void) +template +std::vector TSigmaToNucleonNonEye::getOutput(void) { std::vector out = {}; @@ -156,15 +153,15 @@ std::vector TSigmaToNucleonNonEye:: } // setup /////////////////////////////////////////////////////////////////////// -template -void TSigmaToNucleonNonEye::setup(void) +template +void TSigmaToNucleonNonEye::setup(void) { - envTmpLat(SpinMatrixField1, "c"); + envTmpLat(SpinMatrixField, "c"); } // execution /////////////////////////////////////////////////////////////////// -template -void TSigmaToNucleonNonEye::execute(void) +template +void TSigmaToNucleonNonEye::execute(void) { const Gamma GammaB(Gamma::Algebra::SigmaXZ); // C*gamma_5 const Gamma Id(Gamma::Algebra::Identity); // C*gamma_5 @@ -173,7 +170,7 @@ void TSigmaToNucleonNonEye::execute(void) LOG(Message) << "' with (Gamma^A,Gamma^B)_sigma = ( Identity, C*gamma_5 ) and (Gamma^A,Gamma^B)_nucl = ( Identity, C*gamma_5 )" << std::endl; LOG(Message) << " using sink " << par().sink << "." << std::endl; - envGetTmp(SpinMatrixField1, c); + envGetTmp(SpinMatrixField, c); std::vector buf; std::vector result; @@ -183,11 +180,11 @@ void TSigmaToNucleonNonEye::execute(void) r.info.gammaA_nucl = Id.g; r.info.gammaB_nucl = GammaB.g; - auto &qu_xi = envGet(PropagatorField1, par().qu_xi); - auto &qu_xf = envGet(PropagatorField1, par().qu_xf); - auto &qu_spec = envGet(SlicedPropagator2, par().qu_spec); - auto &qd_xf = envGet(PropagatorField3, par().qd_xf); - auto &qs_xi = envGet(PropagatorField4, par().qs_xi); + auto &qu_xi = envGet(PropagatorField, par().qu_xi); + auto &qu_xf = envGet(PropagatorField, par().qu_xf); + auto &qu_spec = envGet(SlicedPropagator, par().qu_spec); + auto &qd_xf = envGet(PropagatorField, par().qd_xf); + auto &qs_xi = envGet(PropagatorField, par().qs_xi); auto qut = qu_spec[par().xf]; for (auto &G: Gamma::gall) { From 58b6a0d8d1f749d818273e6235ba5cf27b690c4d Mon Sep 17 00:00:00 2001 From: ferben Date: Thu, 5 Dec 2019 16:56:54 +0000 Subject: [PATCH 10/14] changed some naming conditions to resemble rare-kaons --- .../MContraction/SigmaToNucleonEye.hpp | 46 +++++++++--------- .../MContraction/SigmaToNucleonNonEye.hpp | 48 +++++++++---------- 2 files changed, 47 insertions(+), 47 deletions(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index 1aa0b9b7..d8d11f12 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -43,29 +43,29 @@ BEGIN_HADRONS_NAMESPACE /* * 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 * | * \ / 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 * the diagram on the right 'two-trace' * * Propagators: - * * qq_loop - * * qu_spec, source at xi - * * qd_xf, source at xf - * * qs_xi, source at xi + * * qqLoop + * * quSpec, source at ti + * * qdTf, source at tf + * * qsTi, source at ti */ BEGIN_MODULE_NAMESPACE(MContraction) @@ -73,11 +73,11 @@ class SigmaToNucleonEyePar: Serializable { public: GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonEyePar, - std::string, qq_loop, - std::string, qu_spec, - std::string, qd_xf, - std::string, qs_xi, - unsigned int, xf, + std::string, qqLoop, + std::string, quSpec, + std::string, qdTf, + std::string, qsTi, + unsigned int, tf, std::string, sink, std::string, output); }; @@ -134,7 +134,7 @@ TSigmaToNucleonEye::TSigmaToNucleonEye(const std::string name) template std::vector TSigmaToNucleonEye::getInput(void) { - std::vector input = {par().qq_loop, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; + std::vector input = {par().qqLoop, par().quSpec, par().qdTf, par().qsTi, par().sink}; return input; } @@ -175,17 +175,17 @@ void TSigmaToNucleonEye::execute(void) r.info.gammaA_nucl = Id.g; r.info.gammaB_nucl = GammaB.g; - auto &qq_loop = envGet(PropagatorField, par().qq_loop); - auto &qu_spec = envGet(SlicedPropagator, par().qu_spec); - auto &qd_xf = envGet(PropagatorField, par().qd_xf); - auto &qs_xi = envGet(PropagatorField, par().qs_xi); - auto qut = qu_spec[par().xf]; + auto &qqLoop = envGet(PropagatorField, par().qqLoop); + auto &quSpec = envGet(SlicedPropagator, par().quSpec); + auto &qdTf = envGet(PropagatorField, par().qdTf); + auto &qsTi = envGet(PropagatorField, par().qsTi); + auto qut = quSpec[par().tf]; for (auto &G: Gamma::gall) { r.info.gamma_H = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q1",c); + BaryonUtils::Sigma_to_Nucleon_Eye(qqLoop,qut,qdTf,qsTi,G,GammaB,GammaB,"Q1",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) @@ -196,7 +196,7 @@ void TSigmaToNucleonEye::execute(void) result.push_back(r); //Operator Q2, equivalent to the one-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_Eye(qq_loop,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q2",c); + BaryonUtils::Sigma_to_Nucleon_Eye(qqLoop,qut,qdTf,qsTi,G,GammaB,GammaB,"Q2",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index cb73af88..5c9b6812 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -44,7 +44,7 @@ BEGIN_HADRONS_NAMESPACE * Sigma-to-Nucleon 3-pt diagrams, non-eye topologies. * * Schematic: - * qs_xi qu_xf | qs_xi qd_xf + * qsTi quTf | qsTi qdTf * /-->--¬ /-->--¬ | /-->--¬ /-->--¬ * / \ / \ | / \ / \ * / \ / \ | / \ / \ @@ -54,22 +54,22 @@ BEGIN_HADRONS_NAMESPACE * | \ / \ /| | | \ / \ / | * | \ / \ / | | | \ / \ / | * | \ / \ / | | | \-->--/ \-->--/ | - * \ \-->--/ \-->--/ / | \ qu_xi qu_xf / - * \ qu_xi qd_xf / | \ / + * \ \-->--/ \-->--/ / | \ quTi quTf / + * \ quTi qdTf / | \ / * \ / | \ / * \--------->----------/ | \--------->-----------/ - * qu_spec | qu_spec + * quSpec | quSpec * * * analogously to the rare-kaon naming, the left diagram is named 'one-trace' and * the diagram on the right 'two-trace' * * Propagators: - * * qu_xi, source at xi - * * qu_xf, source at xf - * * qu_spec, source at xi - * * qd_xf, source at xf - * * qs_xi, source at xi + * * quTi, source at ti + * * quTf, source at tf + * * quSpec, source at ti + * * qdTf, source at tf + * * qsTi, source at ti */ BEGIN_MODULE_NAMESPACE(MContraction) @@ -77,12 +77,12 @@ class SigmaToNucleonNonEyePar: Serializable { public: GRID_SERIALIZABLE_CLASS_MEMBERS(SigmaToNucleonNonEyePar, - std::string, qu_xi, - std::string, qu_xf, - std::string, qu_spec, - std::string, qd_xf, - std::string, qs_xi, - unsigned int, xf, + std::string, quTi, + std::string, quTf, + std::string, quSpec, + std::string, qdTf, + std::string, qsTi, + unsigned int, tf, std::string, sink, std::string, output); }; @@ -139,7 +139,7 @@ TSigmaToNucleonNonEye::TSigmaToNucleonNonEye(const std::string name) template std::vector TSigmaToNucleonNonEye::getInput(void) { - std::vector input = {par().qu_xi, par().qu_xf, par().qu_spec, par().qd_xf, par().qs_xi, par().sink}; + std::vector input = {par().quTi, par().quTf, par().quSpec, par().qdTf, par().qsTi, par().sink}; return input; } @@ -180,18 +180,18 @@ void TSigmaToNucleonNonEye::execute(void) r.info.gammaA_nucl = Id.g; r.info.gammaB_nucl = GammaB.g; - auto &qu_xi = envGet(PropagatorField, par().qu_xi); - auto &qu_xf = envGet(PropagatorField, par().qu_xf); - auto &qu_spec = envGet(SlicedPropagator, par().qu_spec); - auto &qd_xf = envGet(PropagatorField, par().qd_xf); - auto &qs_xi = envGet(PropagatorField, par().qs_xi); - auto qut = qu_spec[par().xf]; + auto &quTi = envGet(PropagatorField, par().quTi); + auto &quTf = envGet(PropagatorField, par().quTf); + auto &quSpec = envGet(SlicedPropagator, par().quSpec); + auto &qdTf = envGet(PropagatorField, par().qdTf); + auto &qsTi = envGet(PropagatorField, par().qsTi); + auto qut = quSpec[par().tf]; for (auto &G: Gamma::gall) { r.info.gamma_H = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q1",c); + BaryonUtils::Sigma_to_Nucleon_NonEye(quTi,quTf,qut,qdTf,qsTi,G,GammaB,GammaB,"Q1",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) @@ -202,7 +202,7 @@ void TSigmaToNucleonNonEye::execute(void) result.push_back(r); //Operator Q2, equivalent to the one-trace case in the rare-kaons module c=Zero(); - BaryonUtils::Sigma_to_Nucleon_NonEye(qu_xi,qu_xf,qut,qd_xf,qs_xi,G,GammaB,GammaB,"Q2",c); + BaryonUtils::Sigma_to_Nucleon_NonEye(quTi,quTf,qut,qdTf,qsTi,G,GammaB,GammaB,"Q2",c); sliceSum(c,buf,Tp); r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) From a54157e6829615a301e23e5c1e460a379e6791f2 Mon Sep 17 00:00:00 2001 From: ferben Date: Thu, 5 Dec 2019 17:08:09 +0000 Subject: [PATCH 11/14] more definitions changed --- Grid/qcd/utils/BaryonUtils.h | 144 ++++++++++++------------- tests/hadrons/Test_sigma_to_nucleon.cc | 24 ++--- 2 files changed, 84 insertions(+), 84 deletions(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 8e48cf4e..5a46bce1 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -88,18 +88,18 @@ public: template static void Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, robj &result); template - static void Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, - const mobj &Du_xf, + static void Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_ti, + const mobj &Du_tf, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -109,18 +109,18 @@ public: template static void Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, robj &result); template - static void Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, - const mobj &Du_xf, + static void Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_ti, + const mobj &Du_tf, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -129,19 +129,19 @@ public: template static void Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, const mobj &Du_spec, - const PropagatorField &qd_xf, - const PropagatorField &qs_xi, + const PropagatorField &qd_tf, + const PropagatorField &qs_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, const std::string op, SpinMatrixField &stn_corr); template - static void Sigma_to_Nucleon_NonEye(const PropagatorField &qq_xi, - const PropagatorField &qq_xf, + static void Sigma_to_Nucleon_NonEye(const PropagatorField &qq_ti, + const PropagatorField &qq_tf, const mobj &Du_spec, - const PropagatorField &qd_xf, - const PropagatorField &qs_xi, + const PropagatorField &qd_tf, + const PropagatorField &qs_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -327,16 +327,16 @@ void BaryonUtils::ContractBaryons_Sliced(const mobj &D1, * The following code is for Sigma -> N rare hypeon decays * **********************************************************************/ -/* Dq_loop is a quark line from x_H to x_H - * Du_spec is a quark line from x_i to x_f - * Dd_xf is a quark line from x_H to x_f - * Ds_xi is a quark line from x_H to x_i */ +/* Dq_loop is a quark line from t_H to t_H + * Du_spec is a quark line from t_i to t_f + * Dd_tf is a quark line from t_H to t_f + * Ds_ti is a quark line from t_H to t_i */ template template void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -347,7 +347,7 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, auto DuG = Du_spec * GammaB_nucl; // Gamma^B * Ds * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) - auto GDsGDd = GammaB_sigma * Ds_xi * Gamma_H * g5 * adj(Dd_xf) * g5; + auto GDsGDd = GammaB_sigma * Ds_ti * Gamma_H * g5 * adj(Dd_tf) * g5; // Dq_loop * \gamma_\mu^L auto DqG = Dq_loop * Gamma_H; @@ -377,18 +377,18 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, } } -/* Du_xi is a quark line from x_H to x_i - * Du_xf is a quark line from x_H to x_f - * Du_spec is a quark line from x_i to x_f - * Dd_xf is a quark line from x_H to x_f - * Ds_xi is a quark line from x_H to x_i */ +/* Du_ti is a quark line from t_H to t_i + * Du_tf is a quark line from t_H to t_f + * Du_spec is a quark line from t_i to t_f + * Dd_tf is a quark line from t_H to t_f + * Ds_ti is a quark line from t_H to t_i */ template template -void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, - const mobj &Du_xf, +void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_ti, + const mobj &Du_tf, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -398,12 +398,12 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_xi, Gamma g5(Gamma::Algebra::Gamma5); auto DuG = Du_spec * GammaB_nucl; - auto adjDu = g5 * adj(Du_xf) * g5; + auto adjDu = g5 * adj(Du_tf) * g5; auto adjDuG = adjDu * GammaB_nucl; // Gamma^B * Ds * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) - auto GDsGDd = GammaB_sigma * Ds_xi * Gamma_H * g5 * adj(Dd_xf) * g5; + auto GDsGDd = GammaB_sigma * Ds_ti * Gamma_H * g5 * adj(Dd_tf) * g5; // Dq_loop * \gamma_\mu^L - auto DuGH = Du_xi * Gamma_H; + auto DuGH = Du_ti * Gamma_H; for (int ie_n=0; ie_n < 6 ; ie_n++){ int a_n = epsilon[ie_n][0]; //a @@ -441,8 +441,8 @@ template template void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -453,9 +453,9 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, auto DuG = Du_spec * GammaB_nucl; // Gamma^B * Ds * \gamma_\mu^L - auto GDsG = GammaB_sigma * Ds_xi * Gamma_H; + auto GDsG = GammaB_sigma * Ds_ti * Gamma_H; // Dq_loop * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) - auto DqGDd = Dq_loop * Gamma_H * g5 * adj(Dd_xf) * g5; + auto DqGDd = Dq_loop * Gamma_H * g5 * adj(Dd_tf) * g5; for (int ie_n=0; ie_n < 6 ; ie_n++){ int a_n = epsilon[ie_n][0]; //a @@ -485,11 +485,11 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, template template -void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, - const mobj &Du_xf, +void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_ti, + const mobj &Du_tf, const mobj2 &Du_spec, - const mobj &Dd_xf, - const mobj &Ds_xi, + const mobj &Dd_tf, + const mobj &Ds_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, @@ -499,12 +499,12 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_xi, Gamma g5(Gamma::Algebra::Gamma5); auto DuG = Du_spec * GammaB_nucl; - auto adjDu = g5 * adj(Du_xf) * g5; + auto adjDu = g5 * adj(Du_tf) * g5; auto adjDuG = adjDu * GammaB_nucl; // Gamma^B * Ds * \gamma_\mu^L - auto GDsG = GammaB_sigma * Ds_xi * Gamma_H; + auto GDsG = GammaB_sigma * Ds_ti * Gamma_H; // Du * \gamma_\mu^L * (\gamma_5 * Dd^\dagger * \gamma_5) - auto DuGDd = Du_xi * Gamma_H * g5 * adj(Dd_xf) * g5; + auto DuGDd = Du_ti * Gamma_H * g5 * adj(Dd_tf) * g5; for (int ie_n=0; ie_n < 6 ; ie_n++){ int a_n = epsilon[ie_n][0]; //a @@ -541,31 +541,31 @@ template template void BaryonUtils::Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, const mobj &Du_spec, - const PropagatorField &qd_xf, - const PropagatorField &qs_xi, + const PropagatorField &qd_tf, + const PropagatorField &qs_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, const std::string op, SpinMatrixField &stn_corr) { - GridBase *grid = qs_xi.Grid(); + GridBase *grid = qs_ti.Grid(); auto vcorr= stn_corr.View(); auto vq_loop = qq_loop.View(); - auto vd_xf = qd_xf.View(); - auto vs_xi = qs_xi.View(); + auto vd_tf = qd_tf.View(); + auto vs_ti = qs_ti.View(); // accelerator_for(ss, grid->oSites(), grid->Nsimd(), { thread_for(ss,grid->oSites(),{ auto Dq_loop = vq_loop[ss]; - auto Dd_xf = vd_xf[ss]; - auto Ds_xi = vs_xi[ss]; + auto Dd_tf = vd_tf[ss]; + auto Ds_ti = vs_ti[ss]; sobj result=Zero(); if(op == "Q1"){ - Sigma_to_Nucleon_Q1_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); + Sigma_to_Nucleon_Q1_Eye_site(Dq_loop,Du_spec,Dd_tf,Ds_ti,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else if(op == "Q2"){ - Sigma_to_Nucleon_Q2_Eye_site(Dq_loop,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); + Sigma_to_Nucleon_Q2_Eye_site(Dq_loop,Du_spec,Dd_tf,Ds_ti,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else { assert(0 && "Weak Operator not correctly specified"); } @@ -575,36 +575,36 @@ void BaryonUtils::Sigma_to_Nucleon_Eye(const PropagatorField &qq_loop, template template -void BaryonUtils::Sigma_to_Nucleon_NonEye(const PropagatorField &qq_xi, - const PropagatorField &qq_xf, +void BaryonUtils::Sigma_to_Nucleon_NonEye(const PropagatorField &qq_ti, + const PropagatorField &qq_tf, const mobj &Du_spec, - const PropagatorField &qd_xf, - const PropagatorField &qs_xi, + const PropagatorField &qd_tf, + const PropagatorField &qs_ti, const Gamma Gamma_H, const Gamma GammaB_sigma, const Gamma GammaB_nucl, const std::string op, SpinMatrixField &stn_corr) { - GridBase *grid = qs_xi.Grid(); + GridBase *grid = qs_ti.Grid(); auto vcorr= stn_corr.View(); - auto vq_xi = qq_xi.View(); - auto vq_xf = qq_xf.View(); - auto vd_xf = qd_xf.View(); - auto vs_xi = qs_xi.View(); + auto vq_ti = qq_ti.View(); + auto vq_tf = qq_tf.View(); + auto vd_tf = qd_tf.View(); + auto vs_ti = qs_ti.View(); // accelerator_for(ss, grid->oSites(), grid->Nsimd(), { thread_for(ss,grid->oSites(),{ - auto Dq_xi = vq_xi[ss]; - auto Dq_xf = vq_xf[ss]; - auto Dd_xf = vd_xf[ss]; - auto Ds_xi = vs_xi[ss]; + auto Dq_ti = vq_ti[ss]; + auto Dq_tf = vq_tf[ss]; + auto Dd_tf = vd_tf[ss]; + auto Ds_ti = vs_ti[ss]; sobj result=Zero(); if(op == "Q1"){ - Sigma_to_Nucleon_Q1_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); + Sigma_to_Nucleon_Q1_NonEye_site(Dq_ti,Dq_tf,Du_spec,Dd_tf,Ds_ti,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else if(op == "Q2"){ - Sigma_to_Nucleon_Q2_NonEye_site(Dq_xi,Dq_xf,Du_spec,Dd_xf,Ds_xi,Gamma_H,GammaB_sigma,GammaB_nucl,result); + Sigma_to_Nucleon_Q2_NonEye_site(Dq_ti,Dq_tf,Du_spec,Dd_tf,Ds_ti,Gamma_H,GammaB_sigma,GammaB_nucl,result); } else { assert(0 && "Weak Operator not correctly specified"); } diff --git a/tests/hadrons/Test_sigma_to_nucleon.cc b/tests/hadrons/Test_sigma_to_nucleon.cc index 0489438c..b4e79649 100644 --- a/tests/hadrons/Test_sigma_to_nucleon.cc +++ b/tests/hadrons/Test_sigma_to_nucleon.cc @@ -121,24 +121,24 @@ int main(int argc, char *argv[]) MContraction::SigmaToNucleonEye::Par EyePar; EyePar.output = "SigmaToNucleon/Eye_u"; - EyePar.qq_loop = "Qpt_l_loop"; - EyePar.qu_spec = "Qpt_u_spec"; - EyePar.qd_xf = "Qpt_l_4"; - EyePar.qs_xi = "Qpt_s_0"; - EyePar.xf = 4; + EyePar.qqLoop = "Qpt_l_loop"; + EyePar.quSpec = "Qpt_u_spec"; + EyePar.qdTf = "Qpt_l_4"; + EyePar.qsTi = "Qpt_s_0"; + EyePar.tf = 4; EyePar.sink = "sink"; application.createModule("SigmaToNucleonEye_u", EyePar); EyePar.output = "SigmaToNucleon/Eye_c"; - EyePar.qq_loop = "Qpt_c_loop"; + EyePar.qqLoop = "Qpt_c_loop"; application.createModule("SigmaToNucleonEye_c", EyePar); MContraction::SigmaToNucleonNonEye::Par NonEyePar; NonEyePar.output = "SigmaToNucleon/NonEye"; - NonEyePar.qu_xi = "Qpt_l_0"; - NonEyePar.qu_xf = "Qpt_l_4"; - NonEyePar.qu_spec = "Qpt_u_spec"; - NonEyePar.qd_xf = "Qpt_l_4"; - NonEyePar.qs_xi = "Qpt_s_0"; - NonEyePar.xf = 4; + NonEyePar.quTi = "Qpt_l_0"; + NonEyePar.quTf = "Qpt_l_4"; + NonEyePar.quSpec = "Qpt_u_spec"; + NonEyePar.qdTf = "Qpt_l_4"; + NonEyePar.qsTi = "Qpt_s_0"; + NonEyePar.tf = 4; NonEyePar.sink = "sink"; application.createModule("SigmaToNucleonNonEye", NonEyePar); From 164ed9c434ecb050a2bf93b03169c4e21bcd62bb Mon Sep 17 00:00:00 2001 From: ferben Date: Thu, 5 Dec 2019 17:38:00 +0000 Subject: [PATCH 12/14] Naming conventon also applied to metadata --- .../MContraction/SigmaToNucleonEye.hpp | 22 +++++++++---------- .../MContraction/SigmaToNucleonNonEye.hpp | 20 ++++++++--------- 2 files changed, 21 insertions(+), 21 deletions(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp index d8d11f12..b65d0689 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonEye.hpp @@ -94,11 +94,11 @@ public: { public: GRID_SERIALIZABLE_CLASS_MEMBERS(Metadata, - Gamma::Algebra, gamma_H, - Gamma::Algebra, gammaA_sigma, - Gamma::Algebra, gammaB_sigma, - Gamma::Algebra, gammaA_nucl, - Gamma::Algebra, gammaB_nucl, + Gamma::Algebra, gammaH, + Gamma::Algebra, gammaASigma, + Gamma::Algebra, gammaBSigma, + Gamma::Algebra, gammaANucl, + Gamma::Algebra, gammaBNucl, int, trace); }; typedef Correlator Result; @@ -170,10 +170,10 @@ void TSigmaToNucleonEye::execute(void) std::vector result; Result r; - r.info.gammaA_sigma = Id.g; - r.info.gammaB_sigma = GammaB.g; - r.info.gammaA_nucl = Id.g; - r.info.gammaB_nucl = GammaB.g; + r.info.gammaASigma = Id.g; + r.info.gammaBSigma = GammaB.g; + r.info.gammaANucl = Id.g; + r.info.gammaBNucl = GammaB.g; auto &qqLoop = envGet(PropagatorField, par().qqLoop); auto &quSpec = envGet(SlicedPropagator, par().quSpec); @@ -182,7 +182,7 @@ void TSigmaToNucleonEye::execute(void) auto qut = quSpec[par().tf]; for (auto &G: Gamma::gall) { - r.info.gamma_H = G.g; + r.info.gammaH = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); BaryonUtils::Sigma_to_Nucleon_Eye(qqLoop,qut,qdTf,qsTi,G,GammaB,GammaB,"Q1",c); @@ -207,7 +207,7 @@ void TSigmaToNucleonEye::execute(void) result.push_back(r); } - saveResult(par().output, "StN_Eye", result); + saveResult(par().output, "stnEye", result); } diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index 5c9b6812..cbf3665d 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -99,11 +99,11 @@ public: { public: GRID_SERIALIZABLE_CLASS_MEMBERS(Metadata, - Gamma::Algebra, gamma_H, - Gamma::Algebra, gammaA_sigma, - Gamma::Algebra, gammaB_sigma, - Gamma::Algebra, gammaA_nucl, - Gamma::Algebra, gammaB_nucl, + Gamma::Algebra, gammaH, + Gamma::Algebra, gammaASigma, + Gamma::Algebra, gammaBSigma, + Gamma::Algebra, gammaANucl, + Gamma::Algebra, gammaBNucl, int, trace); }; typedef Correlator Result; @@ -175,10 +175,10 @@ void TSigmaToNucleonNonEye::execute(void) std::vector result; Result r; - r.info.gammaA_sigma = Id.g; - r.info.gammaB_sigma = GammaB.g; - r.info.gammaA_nucl = Id.g; - r.info.gammaB_nucl = GammaB.g; + r.info.gammaASigma = Id.g; + r.info.gammaBSigma = GammaB.g; + r.info.gammaANucl = Id.g; + r.info.gammaBNucl = GammaB.g; auto &quTi = envGet(PropagatorField, par().quTi); auto &quTf = envGet(PropagatorField, par().quTf); @@ -213,7 +213,7 @@ void TSigmaToNucleonNonEye::execute(void) result.push_back(r); } - saveResult(par().output, "StN_NonEye", result); + saveResult(par().output, "stnNonEye", result); } From 7ce77690b875c7004e62a8aaaba7c80ac89f8fab Mon Sep 17 00:00:00 2001 From: ferben Date: Thu, 5 Dec 2019 17:38:43 +0000 Subject: [PATCH 13/14] Naming conventon also applied to metadata --- Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp index cbf3665d..b0c86e2e 100644 --- a/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp +++ b/Hadrons/Modules/MContraction/SigmaToNucleonNonEye.hpp @@ -188,7 +188,7 @@ void TSigmaToNucleonNonEye::execute(void) auto qut = quSpec[par().tf]; for (auto &G: Gamma::gall) { - r.info.gamma_H = G.g; + r.info.gammaH = G.g; //Operator Q1, equivalent to the two-trace case in the rare-kaons module c=Zero(); BaryonUtils::Sigma_to_Nucleon_NonEye(quTi,quTf,qut,qdTf,qsTi,G,GammaB,GammaB,"Q1",c); From f7698b93ca57ea3aa4d72b133ad9ca5d1e703661 Mon Sep 17 00:00:00 2001 From: ferben Date: Fri, 6 Dec 2019 09:46:52 +0000 Subject: [PATCH 14/14] corrected comments about quark line directions --- Grid/qcd/utils/BaryonUtils.h | 21 +++++++++++++++------ 1 file changed, 15 insertions(+), 6 deletions(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 5a46bce1..d65b9176 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -329,8 +329,8 @@ void BaryonUtils::ContractBaryons_Sliced(const mobj &D1, /* Dq_loop is a quark line from t_H to t_H * Du_spec is a quark line from t_i to t_f - * Dd_tf is a quark line from t_H to t_f - * Ds_ti is a quark line from t_H to t_i */ + * Dd_tf is a quark line from t_f to t_H + * Ds_ti is a quark line from t_i to t_H */ template template void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, @@ -377,11 +377,11 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_Eye_site(const mobj &Dq_loop, } } -/* Du_ti is a quark line from t_H to t_i - * Du_tf is a quark line from t_H to t_f +/* Du_ti is a quark line from t_i to t_H + * Du_tf is a quark line from t_f to t_H * Du_spec is a quark line from t_i to t_f - * Dd_tf is a quark line from t_H to t_f - * Ds_ti is a quark line from t_H to t_i */ + * Dd_tf is a quark line from t_f to t_H + * Ds_ti is a quark line from t_i to t_H */ template template void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_ti, @@ -437,6 +437,10 @@ void BaryonUtils::Sigma_to_Nucleon_Q1_NonEye_site(const mobj &Du_ti, } //Equivalent to "One-trace" +/* Dq_loop is a quark line from t_H to t_H + * Du_spec is a quark line from t_i to t_f + * Dd_tf is a quark line from t_f to t_H + * Ds_ti is a quark line from t_i to t_H */ template template void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, @@ -483,6 +487,11 @@ void BaryonUtils::Sigma_to_Nucleon_Q2_Eye_site(const mobj &Dq_loop, } } +/* Du_ti is a quark line from t_i to t_H + * Du_tf is a quark line from t_f to t_H + * Du_spec is a quark line from t_i to t_f + * Dd_tf is a quark line from t_f to t_H + * Ds_ti is a quark line from t_i to t_H */ template template void BaryonUtils::Sigma_to_Nucleon_Q2_NonEye_site(const mobj &Du_ti,