field_G_SF.cpp

Go to the documentation of this file.

#include "field_G_SF.h"

//====================================================================

void Field_G_SF::set_boundary_wk(const Mat_SU_N& U)
{
  const int mn = 0;

  if (Communicator::ipe(3) == 0) {
    for (int site = 0; site < m_Svol; ++site) {
      for (int cc = 0; cc < m_Nc * m_Nc; ++cc) {
        field[myindex(2 * cc, site, mn)]     = U.r(cc);
        field[myindex(2 * cc + 1, site, mn)] = U.i(cc);
      }
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_wkpr(const Mat_SU_N& U)
{
  const int mn = 0;

  if (Communicator::ipe(3) == m_NPEt - 1) {
    for (int site = m_Nvol - m_Svol; site < m_Nvol; ++site) {
      for (int cc = 0; cc < m_Nc * m_Nc; ++cc) {
        field[myindex(2 * cc, site, mn)]     = U.r(cc);
        field[myindex(2 * cc + 1, site, mn)] = U.i(cc);
      }
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_zero()
{
  const int mn = 0;

  if (Communicator::ipe(3) == 0) {
    for (int site = 0; site < m_Svol; ++site) {
      for (int cc = 0; cc < 2 * m_Nc * m_Nc; ++cc) {
        field[myindex(cc, site, mn)] = 0.0;
      }
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_spatial_link_zero()
{
  if (Communicator::ipe(3) == 0) {
    for (int mn = 0; mn < 3; ++mn) {
      for (int site = 0; site < m_Svol; ++site) {
        for (int cc = 0; cc < 2 * m_Nc * m_Nc; ++cc) {
          field[myindex(cc, site, mn)] = 0.0;
        }
      }
    }
  }
}


//====================================================================

void Field_G_SF::mult_ct_boundary(const int t, const double ct)
{
  const int mn  = 0;
  const int ini = m_Nx * m_Ny * m_Nz * t;
  const int fin = m_Nx * m_Ny * m_Nz * (t + 1);

  if (((Communicator::ipe(3) == 0) && ((t == 0) || (t == 1))) || ((Communicator::ipe(3) == m_NPEt - 1) && (t == m_Nt - 1))) {
    for (int site = ini; site < fin; ++site) {
      for (int cc = 0; cc < 2 * m_Nc * m_Nc; ++cc) {
        field[myindex(cc, site, mn)] *= ct;
      }
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_wk(Field_G& f)
{
  const int mn = 0;

  assert(f.nex() == 1);
  if (Communicator::ipe(3) == 0) {
    for (int site = 0; site < m_Svol; ++site) {
      f.set_mat(site, mn, m_wk);
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_wkpr(Field_G& f)
{
  const int mn = 0;

  assert(f.nex() == 1);
  if (Communicator::ipe(3) == m_NPEt - 1) {
    for (int site = m_Nvol - m_Svol; site < m_Nvol; ++site) {
      f.set_mat(site, mn, m_wkpr);
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_zero(Field_G& f)
{
  const int mn = 0;

  assert(f.nex() == 1);
  if (Communicator::ipe(3) == 0) {
    for (int site = 0; site < m_Svol; ++site) {
      for (int cc = 0; cc < m_Nc * m_Nc; ++cc) {
        f.set_ri(cc, site, mn, 0.0, 0.0);
      }
    }
  }
}


//====================================================================

void Field_G_SF::set_boundary_spatial_link_zero(Field_G& f)
{
  assert(f.nex() >= 3);
  if (Communicator::ipe(3) == 0) {
    for (int mn = 0; mn < 3; ++mn) {
      for (int site = 0; site < m_Svol; ++site) {
        for (int cc = 0; cc < m_Nc * m_Nc; ++cc) {
          f.set_ri(cc, site, mn, 0.0, 0.0);
        }
      }
    }
  }
}


//====================================================================
void Field_G_SF::set_parameters(const std::vector<double>& phi,
                                const std::vector<double>& phipr)
{
  double aphi[3];
  double aphipr[3];

  for (int i = 0; i < 3; ++i) {
    aphi[i]   = phi[i];
    aphipr[i] = phipr[i];
  }

  set_parameters(aphi, aphipr);
}


//====================================================================
void Field_G_SF::set_parameters(double *phi, double *phipr)
{
  const int Lx = CommonParameters::Lx();

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

  m_wk.zero();
  m_wk.set(0, 0, c0r, c0i);
  m_wk.set(1, 1, c1r, c1i);
  m_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);

  m_wkpr.zero();
  m_wkpr.set(0, 0, c0r, c0i);
  m_wkpr.set(1, 1, c1r, c1i);
  m_wkpr.set(2, 2, c2r, c2i);
}


//====================================================================