portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 13:40:46 +01:00
Code Releases Activity

QedFVol: Fix bug in ScalarVP.cc due to double use of temporary object. Still getting mpi3 errors when configured with enable-comms=mpi[-auto].

Browse Source
This commit is contained in:
James Harrison 2018-01-27 15:15:25 +00:00
parent cdf550845f
commit 4e0cf0cc28
1 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
extras/Hadrons/Modules/MScalar/ScalarVP.cc
View File

@ -270,17 +270,19 @@ void TScalarVP::execute(void)
tmpProp = Cshift(prop0, nu, -1); // S_0(0|x-a\hat{\nu}) tmpProp = Cshift(prop0, nu, -1); // S_0(0|x-a\hat{\nu})
// = S_0(a\hat{\nu}|x) // = S_0(a\hat{\nu}|x)
Usrc = Complex(1.0,0.0); Usrc = Complex(1.0,0.0);
vpContraction(buf, prop0, tmpProp, Usrc, mu); vpContraction(result, prop0, tmpProp, Usrc, mu);
*vpTensor[mu][nu] = buf; *vpTensor[mu][nu] = result;
// Output if necessary // Output if necessary
if (!par().output.empty()) if (!par().output.empty())
{ {
writeVP(writer, buf, writeVP(writer, result,
"Pi_free_"+std::to_string(mu)+"_"+std::to_string(nu)); "Pi_free_"+std::to_string(mu)+"_"+std::to_string(nu));
} }
tmpProp = result; // Just using tmpProp as a temporary ScalarField
// here (buf is modified by calls to writeVP())
// srcT // srcT
result = buf * (-0.5)*q*q*Anu0*Anu0; result = tmpProp * (-0.5)*q*q*Anu0*Anu0;
*vpTensor[mu][nu] += result; *vpTensor[mu][nu] += result;
// Output if necessary // Output if necessary
if (!par().output.empty()) if (!par().output.empty())
@ -290,7 +292,7 @@ void TScalarVP::execute(void)
} }
// snkT // snkT
result = buf * (-0.5)*q*q*Amu*Amu; result = tmpProp * (-0.5)*q*q*Amu*Amu;
*vpTensor[mu][nu] += result; *vpTensor[mu][nu] += result;
// Output if necessary // Output if necessary
if (!par().output.empty()) if (!par().output.empty())