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();
  }
 
  WilsonLoop(const Parameters& params)
    : m_vl(CommonParameters::Vlevel())
  {
    init();
    set_parameters(params);
  }
 
  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; }
 
  void get_parameters(Parameters& params) const;
 
  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