field_F_SF.h

Go to the documentation of this file.

#ifndef FIELD_F_SF_INCLUDED
#define FIELD_F_SF_INCLUDED

#include "Parameters/commonParameters.h"
#include "Communicator/communicator.h"
//#include "field_F.h"
//#include "field_G.h"
#include "Tools/vec_SU_N.h"


class Field_F_SF {
 private:
  int m_Nc2;  // num of the double color elements
  int m_Nvol; // lattice volume
  int m_Nd;   // num of the spinor elements
  int m_Nin;  // internal d.o.f.
  int m_Svol;

 public:
  Field_F_SF() :
    m_Nvol(CommonParameters::Nvol()),
    m_Nd(CommonParameters::Nd())
  {
    m_Nc2  = 2 * CommonParameters::Nc();
    m_Svol = m_Nvol / CommonParameters::Nt();
  }

  void set_boundary_zero(Field& f)
  {
    if (Communicator::ipe(3) == 0) {
      for (int site = 0; site < m_Svol; ++site) {
        for (int s = 0; s < m_Nd; ++s) {
          for (int cc = 0; cc < m_Nc2; ++cc) {
            f.set(cc + m_Nc2 * s, site, 0, 0.0);
          }
        }
      }
    }
  }
};
#endif