wilsonLoop.h

Go to the documentation of this file.

#ifndef WILSONLOOP_INCLUDED
#define WILSONLOOP_INCLUDED

#include "defs.h"
#include "parameters.h"
#include "field_G.h"
#include "shiftField_eo.h"

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


//- parameters class
class Parameters_WilsonLoop : virtual public Parameters
{
 public:
  Parameters_WilsonLoop();
};
//- end

class WilsonLoop
{
 protected:
  Bridge::VerboseLevel m_vl;

 private:
  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::valarray<int>   unit_vec;
  std::valarray<unit_vec> m_Nunit;
  std::valarray<int>      m_Nmax;

 public:

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

  virtual ~WilsonLoop() {}

  virtual void set_parameters(const Parameters& params);

  void set_parameters(int Nspc_size, int Ntmp_size, int Ntype);

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

  double measure(Field_G& U);

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

 private:

  void init();

  double calc_wloop(Field_G& Uspc, int t_ext);

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

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

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