field_G_SF.h

Go to the documentation of this file.

#ifndef FIELD_G_SF_INCLUDED
#define FIELD_G_SF_INCLUDED

#include "Parameters/parameters.h"
#include "field_G.h"


class Field_G_SF : public Field_G
{
 private:
  int m_Nc;
  int m_Ndf;
  int m_Nvol;
  int m_Nex;

  int Nx, Ny, Nz, Nt, Svol;
  int NPEt;

  Mat_SU_N wk;
  Mat_SU_N wkpr;

 public:

  //  Field_G_SF(){}
  Field_G_SF(const int Nvol = CommonParameters::Nvol(), const int Nex = 1) :
    m_Nc(CommonParameters::Nc()), m_Nvol(Nvol), m_Nex(Nex),
    wk(m_Nc), wkpr(m_Nc)
  {
    Nx    = CommonParameters::Nx();
    Ny    = CommonParameters::Ny();
    Nz    = CommonParameters::Nz();
    Nt    = CommonParameters::Nt();
    NPEt  = CommonParameters::NPEt();
    Svol  = Nx * Ny * Nz;
    m_Ndf = 2 * m_Nc * m_Nc;
    Field::reset(m_Ndf, m_Nvol, m_Nex);
  }

  Field_G_SF(const Field& x) :
    Field_G(x),
    m_Nc(CommonParameters::Nc()), m_Ndf(x.nin()), m_Nvol(x.nvol()), m_Nex(x.nex()),
    wk(m_Nc), wkpr(m_Nc)
  {
    assert(m_Ndf == 2 * m_Nc * m_Nc);
    Nx   = CommonParameters::Nx();
    Ny   = CommonParameters::Ny();
    Nz   = CommonParameters::Nz();
    Nt   = CommonParameters::Nt();
    NPEt = CommonParameters::NPEt();
    Svol = Nx * Ny * Nz;
  }

  Field_G_SF(double *phi, double *phipr) :
    m_Nc(CommonParameters::Nc()), wk(m_Nc), wkpr(m_Nc)
  {
    int Lx = CommonParameters::Lx();

    m_Nex  = 1;
    m_Nvol = CommonParameters::Nvol();
    Nx     = CommonParameters::Nx();
    Ny     = CommonParameters::Ny();
    Nz     = CommonParameters::Nz();
    Nt     = CommonParameters::Nt();
    NPEt   = CommonParameters::NPEt();
    Svol   = Nx * Ny * Nz;
    m_Ndf  = 2 * m_Nc * m_Nc;

    double c0r, c0i, c1r, c1i, c2r, c2i;
    c0r = cos(phi[0] / Lx);
    c0i = sin(phi[0] / Lx);
    c1r = cos(phi[1] / Lx);
    c1i = sin(phi[1] / Lx);
    c2r = cos(phi[2] / Lx);
    c2i = sin(phi[2] / Lx);
    wk.zero();
    wk.set(0, 0, c0r, c0i);
    wk.set(1, 1, c1r, c1i);
    wk.set(2, 2, c2r, c2i);

    c0r = cos(phipr[0] / Lx);
    c0i = sin(phipr[0] / Lx);
    c1r = cos(phipr[1] / Lx);
    c1i = sin(phipr[1] / Lx);
    c2r = cos(phipr[2] / Lx);
    c2i = sin(phipr[2] / Lx);
    wkpr.zero();
    wkpr.set(0, 0, c0r, c0i);
    wkpr.set(1, 1, c1r, c1i);
    wkpr.set(2, 2, c2r, c2i);
  }

  Field_G_SF& operator=(const Field_G_SF& v) { copy(*this, v); return *this; }

  void set_parameters(const std::vector<double>& phi,
                      const std::vector<double>& phipr);

  void set_parameters(double *phi, double *phipr);

  void set_boundary_wk(const Mat_SU_N& U);

  void set_boundary_wkpr(const Mat_SU_N& U);

  void set_boundary_zero();

  void set_boundary_spatial_link_zero();

  void mult_ct_boundary(int t, double ct);


  void set_boundary_wk(Field_G& f);

  void set_boundary_wkpr(Field_G& f);

  void set_boundary_zero(Field_G& f);

  void set_boundary_spatial_link_zero(Field_G& f);
};
#endif