portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Hadrons: non-eye weak 3pt fix

Browse Source
This commit is contained in:
Antonin Portelli 2019-02-05 11:35:10 +00:00
parent 9c31305b8d
commit 90d6d28547
1 changed files with 6 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
Hadrons/Modules/MContraction/WeakNonEye3pt.hpp
View File

@ -8,8 +8,7 @@
BEGIN_HADRONS_NAMESPACE BEGIN_HADRONS_NAMESPACE
/* /*
* Weak Hamiltonian meson 3-pt diagrams, non-eye-topologies. * Weak Hamiltonian meson 3-pt diagrams, non-eye topologies.
*
* *
* Schematic: * Schematic:
* qbl qbr | qbl qbr * qbl qbr | qbl qbr
@ -24,10 +23,10 @@ BEGIN_HADRONS_NAMESPACE
* \ / \ / | \-->--/ \-->--/ * \ / \ / | \-->--/ \-->--/
* \-->--/ \-->--/ | ql qr * \-->--/ \-->--/ | ql qr
* ql qr | * ql qr |
* One trace | Two traces * one trace | two traces
* *
* One trace : tr(ql*adj(gIn)*g5*adj(qbl)*g5*G*qbr*gOut*g5*adj(qr)*g5*G) * one trace : tr(ql*adj(gIn)*g5*adj(qbl)*g5*G*qbr*gOut*g5*adj(qr)*g5*G)
* Two traces: tr(ql*adj(gIn)*g5*adj(qbl)*g5*G)*tr(qbr*gOut*g5*adj(qr)*g5*G) * two traces: tr(ql*adj(gIn)*g5*adj(qbl)*g5*G)*tr(qbr*gOut*g5*adj(qr)*g5*G)
* *
*/ */
@ -143,6 +142,7 @@ void TWeakNonEye3pt<FImpl>::execute(void)
// one trace // one trace
corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G*qbr*gOut*g5*adj(qr)*g5*G); corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G*qbr*gOut*g5*adj(qr)*g5*G);
sliceSum(corr, buf, Tp); sliceSum(corr, buf, Tp);
r.corr.clear();
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)
{ {
r.corr.push_back(TensorRemove(buf[t])); r.corr.push_back(TensorRemove(buf[t]));
@ -152,6 +152,7 @@ void TWeakNonEye3pt<FImpl>::execute(void)
// two traces // two traces
corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G)*trace(qbr*gOut*g5*adj(qr)*g5*G); corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G)*trace(qbr*gOut*g5*adj(qr)*g5*G);
sliceSum(corr, buf, Tp); sliceSum(corr, buf, Tp);
r.corr.clear();
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)
{ {
r.corr.push_back(TensorRemove(buf[t])); r.corr.push_back(TensorRemove(buf[t]));