From 90d6d2854769f1867d080d4de12c06de26f47244 Mon Sep 17 00:00:00 2001 From: Antonin Portelli Date: Tue, 5 Feb 2019 11:35:10 +0000 Subject: [PATCH] Hadrons: non-eye weak 3pt fix --- Hadrons/Modules/MContraction/WeakNonEye3pt.hpp | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp b/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp index dc74192b..b755bf88 100644 --- a/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp +++ b/Hadrons/Modules/MContraction/WeakNonEye3pt.hpp @@ -8,8 +8,7 @@ BEGIN_HADRONS_NAMESPACE /* - * Weak Hamiltonian meson 3-pt diagrams, non-eye-topologies. - * + * Weak Hamiltonian meson 3-pt diagrams, non-eye topologies. * * Schematic: * qbl qbr | qbl qbr @@ -24,10 +23,10 @@ BEGIN_HADRONS_NAMESPACE * \ / \ / | \-->--/ \-->--/ * \-->--/ \-->--/ | 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) - * Two traces: tr(ql*adj(gIn)*g5*adj(qbl)*g5*G)*tr(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) * */ @@ -143,6 +142,7 @@ void TWeakNonEye3pt::execute(void) // one trace corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G*qbr*gOut*g5*adj(qr)*g5*G); sliceSum(corr, buf, Tp); + r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) { r.corr.push_back(TensorRemove(buf[t])); @@ -152,6 +152,7 @@ void TWeakNonEye3pt::execute(void) // two traces corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G)*trace(qbr*gOut*g5*adj(qr)*g5*G); sliceSum(corr, buf, Tp); + r.corr.clear(); for (unsigned int t = 0; t < buf.size(); ++t) { r.corr.push_back(TensorRemove(buf[t]));