Hadrons: 4-quark non-eye 3-pt contractions

Hadrons/Modules.hpp

@@ -4,6 +4,7 @@
 #include <Hadrons/Modules/MContraction/Meson.hpp>
 #include <Hadrons/Modules/MContraction/DiscLoop.hpp>
 #include <Hadrons/Modules/MContraction/Gamma3pt.hpp>
+#include <Hadrons/Modules/MContraction/WeakNonEye3pt.hpp>
 #include <Hadrons/Modules/MFermion/FreeProp.hpp>
 #include <Hadrons/Modules/MFermion/GaugeProp.hpp>
 #include <Hadrons/Modules/MSource/SeqGamma.hpp>

Hadrons/Modules/MContraction/WeakNonEye3pt.cc

@@ -0,0 +1,7 @@
+#include <Hadrons/Modules/MContraction/WeakNonEye3pt.hpp>
+
+using namespace Grid;
+using namespace Hadrons;
+using namespace MContraction;
+
+template class Grid::Hadrons::MContraction::TWeakNonEye3pt<FIMPL>;

Hadrons/Modules/MContraction/WeakNonEye3pt.hpp

@@ -0,0 +1,169 @@
+#ifndef Hadrons_MContraction_WeakNonEye3pt_hpp_
+#define Hadrons_MContraction_WeakNonEye3pt_hpp_
+
+#include <Hadrons/Global.hpp>
+#include <Hadrons/Module.hpp>
+#include <Hadrons/ModuleFactory.hpp>
+
+BEGIN_HADRONS_NAMESPACE
+
+/*
+ * Weak Hamiltonian meson 3-pt diagrams, non-eye-topologies.
+ * 
+ * 
+ * Schematic:     
+ *            qbl           qbr            |            qbl             qbr
+ *          /--<--¬       /--<--¬          |          /--<--¬         /--<--¬       
+ *         /       \     /       \         |         /       \       /       \      
+ *        /         \   /         \        |        /         \     /         \     
+ *       /           \ /           \       |       /           \   /           \    
+ *  gIn *             * G           * gOut |  gIn *           G * * G           * gOut
+ *      \             * G           |      |       \           /   \           /
+ *       \           / \           /       |        \         /     \         /    
+ *        \         /   \         /        |         \       /       \       /  
+ *         \       /     \       /         |          \-->--/         \-->--/      
+ *          \-->--/       \-->--/          |            ql              qr 
+ *            ql            qr             |
+ *               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)
+ * 
+ */
+
+BEGIN_MODULE_NAMESPACE(MContraction)
+
+class WeakNonEye3ptPar: Serializable
+{
+public:
+    GRID_SERIALIZABLE_CLASS_MEMBERS(WeakNonEye3ptPar,
+                                    std::string,    qLeft,
+                                    std::string,    qBarLeft,
+                                    std::string,    qRight,
+                                    std::string,    qBarRight,
+                                    Gamma::Algebra, gammaIn,
+                                    Gamma::Algebra, gammaOut,
+                                    std::string,    output);
+};
+
+template <typename FImpl>
+class TWeakNonEye3pt: public Module<WeakNonEye3ptPar>
+{
+public:
+    FERM_TYPE_ALIASES(FImpl,);
+    class Metadata: Serializable
+    {
+    public:
+        GRID_SERIALIZABLE_CLASS_MEMBERS(Metadata,
+                                        Gamma::Algebra, in,
+                                        Gamma::Algebra, out,
+                                        Gamma::Algebra, op,
+                                        unsigned int,   trace);
+    };
+    typedef Correlator<Metadata> Result;
+public:
+    // constructor
+    TWeakNonEye3pt(const std::string name);
+    // destructor
+    virtual ~TWeakNonEye3pt(void) {};
+    // dependency relation
+    virtual std::vector<std::string> getInput(void);
+    virtual std::vector<std::string> getOutput(void);
+    // setup
+    virtual void setup(void);
+    // execution
+    virtual void execute(void);
+};
+
+MODULE_REGISTER_TMP(WeakNonEye3pt, TWeakNonEye3pt<FIMPL>, MContraction);
+
+/******************************************************************************
+ *                 TWeakNonEye3pt implementation                             *
+ ******************************************************************************/
+// constructor /////////////////////////////////////////////////////////////////
+template <typename FImpl>
+TWeakNonEye3pt<FImpl>::TWeakNonEye3pt(const std::string name)
+: Module<WeakNonEye3ptPar>(name)
+{}
+
+// dependencies/products ///////////////////////////////////////////////////////
+template <typename FImpl>
+std::vector<std::string> TWeakNonEye3pt<FImpl>::getInput(void)
+{
+    std::vector<std::string> in = {par().qLeft, par().qBarLeft, 
+                                   par().qRight, par().qBarRight};
+    
+    return in;
+}
+
+template <typename FImpl>
+std::vector<std::string> TWeakNonEye3pt<FImpl>::getOutput(void)
+{
+    std::vector<std::string> out = {};
+    
+    return out;
+}
+
+// setup ///////////////////////////////////////////////////////////////////////
+template <typename FImpl>
+void TWeakNonEye3pt<FImpl>::setup(void)
+{
+    envTmpLat(ComplexField, "corr");
+}
+
+// execution ///////////////////////////////////////////////////////////////////
+template <typename FImpl>
+void TWeakNonEye3pt<FImpl>::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) << "ql  : " << par().qLeft << std::endl;
+    LOG(Message) << "qbl : " << par().qBarLeft << std::endl;
+    LOG(Message) << "qr  : " << par().qRight << std::endl;
+    LOG(Message) << "qbr : " << par().qBarRight << std::endl;
+
+    std::vector<Result> result;
+    Result              r;
+    auto                &ql  = envGet(PropagatorField, par().qLeft);
+    auto                &qbl = envGet(PropagatorField, par().qBarLeft);
+    auto                &qr  = envGet(PropagatorField, par().qRight);
+    auto                &qbr = envGet(PropagatorField, par().qBarRight);
+    Gamma               gIn(par().gammaIn), gOut(par().gammaOut);
+    Gamma               g5(Gamma::Algebra::Gamma5);
+
+    envGetTmp(ComplexField, corr);
+    r.info.in  = par().gammaIn;
+    r.info.out = par().gammaOut;
+    for (auto &G: Gamma::gall)
+    {
+        SlicedComplex buf;
+
+        r.info.op = G.g;
+        // one trace
+        corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G*qbr*gOut*g5*adj(qr)*g5*G);
+        sliceSum(corr, buf, Tp);
+        for (unsigned int t = 0; t < buf.size(); ++t)
+        {
+            r.corr.push_back(TensorRemove(buf[t]));
+        }
+        r.info.trace = 1;
+        result.push_back(r);
+        // two traces
+        corr = trace(ql*adj(gIn)*g5*adj(qbl)*g5*G)*trace(qbr*gOut*g5*adj(qr)*g5*G);
+        sliceSum(corr, buf, Tp);
+        for (unsigned int t = 0; t < buf.size(); ++t)
+        {
+            r.corr.push_back(TensorRemove(buf[t]));
+        }
+        r.info.trace = 2;
+        result.push_back(r);
+    }
+    saveResult(par().output, "weakNonEye3pt", result);
+}
+
+END_MODULE_NAMESPACE
+
+END_HADRONS_NAMESPACE
+
+#endif // Hadrons_MContraction_WeakNonEye3pt_hpp_

Hadrons/modules.inc

@@ -1,6 +1,7 @@
 modules_cc =\
   Modules/MContraction/Baryon.cc \
   Modules/MContraction/Meson.cc \
+  Modules/MContraction/WeakNonEye3pt.cc \
   Modules/MContraction/A2AAslashField.cc \
   Modules/MContraction/A2AMesonField.cc \
   Modules/MContraction/DiscLoop.cc \