noiseVector_Z2.cpp

Go to the documentation of this file.

#include "noiseVector_Z2.h"

//====================================================================
void NoiseVector_Z2::set(Field& v)
{
  // This implementation assumes the given field v is complex field.

  m_rand->uniform_lex_global(v);

  int    Nex  = v.nex();
  int    Nvol = v.nvol();
  int    Nin  = v.nin();
  double rn, rn2, rn3;
  double RF2 = 1.0 / sqrt(2.0);

  for (int ex = 0; ex < Nex; ++ex) {
    for (int site = 0; site < Nvol; ++site) {
      for (int in = 0; in < Nin; ++in) {
        rn  = v.cmp(in, site, ex);
        rn2 = floor(2.0 * rn);
        rn3 = (2.0 * rn2 - 1.0) * RF2;
        v.set(in, site, ex, rn3);
      }
    }
  }
}


//====================================================================
//============================================================END=====