index_eo.h

Go to the documentation of this file.

 
#ifndef INDEX_EO_INCLUDED
#define INDEX_EO_INCLUDED
 
#include <vector>
 
#include "Field/index_lex.h"
#include "Field/field.h"
 
#include "Communicator/communicator.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class Index_eo {
 public:
  static const std::string class_name;
 
 private:
  Bridge::VerboseLevel m_vl;
 
  int m_Nx, m_Ny, m_Nz, m_Nt, m_Nvol;
  int m_Nx2, m_Nvol2;
  int m_node_eo;   
 
  std::vector<int> m_yzt_eo;
  std::vector<int> m_site_up;
  std::vector<int> m_site_dn;
 
  Index_lex m_index_lex;
 
 public:
 
  Index_eo() { init(); }
 
  int leo(const int y, const int z, const int t) const
  { return m_yzt_eo[y + m_Ny * (z + m_Nz * t)]; }
 
  int site(const int x2, const int y, const int z, const int t,
           const int ieo) const
  { return x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo; }
 
  int site(const int is, const int ieo) const
  { return is + m_Nvol2 * ieo; }
 
  int site_up(const int x2, const int y, const int z, const int t,
              const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_up[s] + m_Nvol2 * (1 - ieo);
  }
 
  int site_xup(const int x2, const int y, const int z, const int t,
               const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_up[s] + m_Nvol2 * (1 - ieo);
  }
 
  int site_dn(const int x2, const int y, const int z, const int t,
              const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_dn[s] + m_Nvol2 * (1 - ieo);
  }
 
  int site_xdn(const int x2, const int y, const int z, const int t,
               const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_dn[s] + m_Nvol2 * (1 - ieo);
  }
 
  int siteh(const int x2, const int y, const int z, const int t)
  const
  { return x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)); }
 
  int siteh_up(const int x2, const int y, const int z, const int t,
               const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_up[s];
  }
 
  int siteh_xup(const int x2, const int y, const int z, const int t,
                const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_up[s];
  }
 
  int siteh_dn(const int x2, const int y, const int z, const int t,
               const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_dn[s];
  }
 
  int siteh_xdn(const int x2, const int y, const int z, const int t,
                const int ieo) const
  {
    int s = x2 + m_Nx2 * (y + m_Ny * (z + m_Nz * t)) + m_Nvol2 * ieo;
    return m_site_dn[s];
  }
 
  void init();
 
  void convertField(Field& eo, const Field& lex);
  void convertField(Field& eo, const Field& lex, const int ieo);
 
  void reverseField(Field& lex, const Field& eo);
  void reverseField(Field& lex, const Field& eo, const int ieo);
 
  void splitField(Field& e, Field& o, const Field& eo);
 
  void mergeField(Field& eo, const Field& e, const Field& o);
};
#endif