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

beauty

Browse Source
This commit is contained in:
ferben 2019-11-13 15:14:51 +00:00
parent ee9dd22643
commit 500ef17143
1 changed files with 12 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
Hadrons/Modules/MDistil/Noises.hpp
View File

@ -108,13 +108,20 @@ void TNoises<FImpl>::execute(void)
// We use our own seeds so we can specify different noises per quark // We use our own seeds so we can specify different noises per quark
Real rn; Real rn;
auto &noise = envGet(NoiseTensor, getName()); auto &noise = envGet(NoiseTensor, getName());
for (int inoise = 0; inoise < dp.nnoise; inoise++) { for (int inoise = 0; inoise < dp.nnoise; inoise++)
for (int t = 0; t < Nt; t++) { {
for (int ivec = 0; ivec < dp.nvec; ivec++) { for (int t = 0; t < Nt; t++)
for (int is = 0; is < Ns; is++) { {
for (int ivec = 0; ivec < dp.nvec; ivec++)
{
for (int is = 0; is < Ns; is++)
{
if (exact_distillation) if (exact_distillation)
{
noise.tensor(inoise, t, ivec, is) = 1.; noise.tensor(inoise, t, ivec, is) = 1.;
else{ }
else
{
random(rngSerial(),rn); random(rngSerial(),rn);
// We could use a greater number of complex roots of unity // We could use a greater number of complex roots of unity
// ... but this seems to work well // ... but this seems to work well