wilsonLoop.h

Go to the documentation of this file.

#ifndef WILSONLOOP_INCLUDED
#define WILSONLOOP_INCLUDED

#include "Field/field_G.h"
#include "Field/shiftField_eo.h"

#include "Parameters/parameters.h"

#include "IO/bridgeIO.h"
using Bridge::vout;


class WilsonLoop
{
 public:
  static const std::string class_name;

 protected:
  Bridge::VerboseLevel m_vl;

 private:
  std::string m_filename_output;

  int m_Nspc_size;  
  int m_Ntmp_size;  
  int m_Ntype;      

  int m_Ntype_max;  
  int m_Nx_ext;     
  int m_Ny_ext;     
  int m_Nz_ext;     
  int m_Nt_ext;     
  int m_Nvol_ext;   

  typedef std::vector<int> unit_vec;
  std::vector<unit_vec> m_Nunit;
  std::vector<int> m_Nmax;

 public:
  WilsonLoop()
    : m_vl(CommonParameters::Vlevel())
  {
    init();
  }

  virtual ~WilsonLoop() {}

 private:
  // non-copyable
  WilsonLoop(const WilsonLoop&);
  WilsonLoop& operator=(const WilsonLoop&);

 public:
  virtual void set_parameters(const Parameters& params);
  void set_parameters(const int Nspc_size, const int Ntmp_size, const int Ntype);

  void set_parameter_verboselevel(const Bridge::VerboseLevel vl) { m_vl = vl; }

  double measure(const Field_G& U);

  int index_wloop(const int i_spc, const int i_tmp, const int i_type)
  {
    return i_spc + m_Nspc_size * (i_tmp + m_Ntmp_size * i_type);
  }

 private:
  void init();

  double calc_wloop(const Field_G& Uspc, const int t_ext);

  void redef_Uspc(Field_G& Uspc, const Field_G& Uext,
                  const int j, const int nu, const std::vector<int>& unit_v);

  void set_extfield(Field_G& Uext, const Field_G& Uorg);

  void gfix_temporal(Field_G& Uext);
};
#endif