From bdd2765461d839f6d167d368b1a2c4cc311bf45f Mon Sep 17 00:00:00 2001 From: Lanny91 Date: Tue, 7 Feb 2017 13:06:42 +0000 Subject: [PATCH] Added missing allocation of Weak Hamiltonian result vector --- extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.cc | 2 +- extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.hpp | 3 ++- extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.cc | 2 +- extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp | 3 ++- 4 files changed, 6 insertions(+), 4 deletions(-) diff --git a/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.cc b/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.cc index 0df6ec2f..6c13aee5 100644 --- a/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.cc +++ b/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.cc @@ -101,7 +101,7 @@ void TWeakHamiltonianEye::execute(void) Gamma g5 = Gamma(Gamma::Algebra::Gamma5); LatticeComplex expbuf(env().getGrid()); std::vector corrbuf; - std::vector result; + std::vector result(n_eye_diag); unsigned int ndim = env().getNd(); PropagatorField tmp1(env().getGrid()); diff --git a/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.hpp b/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.hpp index 6f5a95c0..a07d814f 100644 --- a/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.hpp +++ b/extras/Hadrons/Modules/MContraction/WeakHamiltonianEye.hpp @@ -41,7 +41,8 @@ BEGIN_MODULE_NAMESPACE(MContraction) enum { S_diag = 0, - E_diag = 1 + E_diag = 1, + n_eye_diag = 2 }; // Saucer and Eye subdiagram contractions. diff --git a/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.cc b/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.cc index 148b9afa..98a5e104 100644 --- a/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.cc +++ b/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.cc @@ -103,7 +103,7 @@ void TWeakHamiltonianNonEye::execute(void) Gamma g5 = Gamma(Gamma::Algebra::Gamma5); LatticeComplex expbuf(env().getGrid()); std::vector corrbuf; - std::vector result; + std::vector result(n_noneye_diag); unsigned int ndim = env().getNd(); PropagatorField tmp1(env().getGrid()); diff --git a/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp b/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp index f0d26b54..a594cfda 100644 --- a/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp +++ b/extras/Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp @@ -41,7 +41,8 @@ BEGIN_MODULE_NAMESPACE(MContraction) enum { W_diag = 0, - C_diag = 1 + C_diag = 1, + n_noneye_diag = 2 }; // Wing and Connected subdiagram contractions