From c504b4dbad611b8f36599fb5d6202a85b465134d Mon Sep 17 00:00:00 2001 From: Lanny91 Date: Mon, 5 Jun 2017 15:56:43 +0100 Subject: [PATCH] Code cleaning --- extras/Hadrons/Modules/MContraction/WardIdentity.hpp | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/extras/Hadrons/Modules/MContraction/WardIdentity.hpp b/extras/Hadrons/Modules/MContraction/WardIdentity.hpp index 7fc7d15d..fb2ea173 100644 --- a/extras/Hadrons/Modules/MContraction/WardIdentity.hpp +++ b/extras/Hadrons/Modules/MContraction/WardIdentity.hpp @@ -128,12 +128,11 @@ void TWardIdentity::execute(void) LOG(Message) << "Performing Ward Identity checks for quark '" << par().q << "'." << std::endl; - PropagatorField psi(env().getGrid()), tmp(env().getGrid()), - vector_WI(env().getGrid()); + PropagatorField tmp(env().getGrid()), vector_WI(env().getGrid()); PropagatorField &q = *env().template getObject(par().q); FMat &act = *(env().template getObject(par().action)); Gamma g5(Gamma::Algebra::Gamma5); - + // Compute D_mu V_mu, D here is backward derivative. vector_WI = zero; for (unsigned int mu = 0; mu < Nd; ++mu) @@ -149,6 +148,7 @@ void TWardIdentity::execute(void) if (par().test_axial) { + PropagatorField psi(env().getGrid()); LatticeComplex PP(env().getGrid()), axial_defect(env().getGrid()), PJ5q(env().getGrid()); std::vector axial_buf; @@ -159,7 +159,7 @@ void TWardIdentity::execute(void) { act.ContractConservedCurrent(q, q, tmp, Current::Axial, mu); tmp -= Cshift(tmp, mu, -1); - axial_defect += 2.*trace(g5*tmp); + axial_defect += trace(g5*tmp); } // Get for 5D (zero for 4D) and . @@ -191,7 +191,7 @@ void TWardIdentity::execute(void) LOG(Message) << "|PJ5q|^2 = " << norm2(PJ5q) << std::endl; LOG(Message) << "Axial Ward Identity defect Delta_mu A_mu = " << norm2(axial_defect) << std::endl; - + // Axial defect by timeslice. axial_defect -= 2.*(par().mass*PP + PJ5q); LOG(Message) << "Check Axial defect by timeslice" << std::endl;