portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Make use of GammaL class in Weak Hamiltonian contractions

Browse Source
This commit is contained in:
Lanny91 2017-02-08 12:45:39 +00:00
parent 6ebf8b12b6
commit beba824136
4 changed files with 6 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.cc
View File

@ -114,10 +114,8 @@ void TWeakHamiltonianEye::execute(void)
// Setup for S-type contractions.
for (int mu = 0; mu < ndim; ++mu)
{
S_body[mu] = MAKE_SE_BODY(q1, q2, q3, Gamma::gmu[mu]) -
MAKE_SE_BODY(q1, q2, q3, Gamma::gmu[mu]*g5);
S_loop[mu] = MAKE_SE_LOOP(q4, Gamma::gmu[mu]) -
MAKE_SE_LOOP(q4, Gamma::gmu[mu]*g5);
S_body[mu] = MAKE_SE_BODY(q1, q2, q3, GammaL(Gamma::gmu[mu]));
S_loop[mu] = MAKE_SE_LOOP(q4, GammaL(Gamma::gmu[mu]));
}
// Perform S-type contractions.

2
extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.hpp
View File

@ -46,7 +46,7 @@ enum
};
// Saucer and Eye subdiagram contractions.
#define MAKE_SE_BODY(Q_1, Q_2, Q_3, gamma) (Q_3*g5*Q_1*adj(Q_2)*(g5*gamma))
#define MAKE_SE_BODY(Q_1, Q_2, Q_3, gamma) (Q_3*g5*Q_1*adj(Q_2)*g5*gamma)
#define MAKE_SE_LOOP(Q_loop, gamma) (Q_loop*gamma)
class TWeakHamiltonianEye: public Module<WeakHamiltonianPar>

6
extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.cc
View File

@ -116,10 +116,8 @@ void TWeakHamiltonianNonEye::execute(void)
// Setup for C-type contractions.
for (int mu = 0; mu < ndim; ++mu)
{
C_i_side_loop[mu] = MAKE_CW_SUBDIAG(q1, q2, Gamma::gmu[mu]) -
MAKE_CW_SUBDIAG(q1, q2, Gamma::gmu[mu]*g5);
C_f_side_loop[mu] = MAKE_CW_SUBDIAG(q3, q4, Gamma::gmu[mu]) -
MAKE_CW_SUBDIAG(q3, q4, Gamma::gmu[mu]*g5);
C_i_side_loop[mu] = MAKE_CW_SUBDIAG(q1, q2, GammaL(Gamma::gmu[mu]));
C_f_side_loop[mu] = MAKE_CW_SUBDIAG(q3, q4, GammaL(Gamma::gmu[mu]));
}
// Perform C-type contractions.

2
extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp
View File

@ -46,7 +46,7 @@ enum
};
// Wing and Connected subdiagram contractions
#define MAKE_CW_SUBDIAG(Q_1, Q_2, gamma) (Q_1*adj(Q_2)*(g5*gamma))
#define MAKE_CW_SUBDIAG(Q_1, Q_2, gamma) (Q_1*adj(Q_2)*g5*gamma)
class TWeakHamiltonianNonEye: public Module<WeakHamiltonianPar>
{