polyakovLoop.h

Go to the documentation of this file.

 
#ifndef POLYAKOVLOOP_INCLUDED
#define POLYAKOVLOOP_INCLUDED
 
#include <cassert>
#include "Parameters/parameters.h"
#include "Field/field_G.h"
#include "Field/shiftField_lex.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class PolyakovLoop
{
 public:
  static const std::string class_name;
 
 protected:
  Bridge::VerboseLevel m_vl;
 
 private:
  std::string m_filename_output;
 
  int m_Nspc_size;  
  int m_Ntype;      
 
  Field_G m_P, m_Pcp1, m_Pcp2, m_Ut;
 
  // int m_Ntype_max;  //!< maximum size of loop-type
  // int m_Nx_ext;     //!< size of extended gauge config.
  // int m_Ny_ext;     //!< size of extended gauge config.
  // int m_Nz_ext;     //!< size of extended gauge config.
  // int m_Nt_ext;     //!< size of extended gauge config.
  // int m_Nvol_ext;   //!< volume of extended gauge config.
  //
  // typedef std::vector<int> unitvec;
  // std::vector<unitvec> m_Nunit;
  // std::vector<int> m_Nmax;
 
 public:
  PolyakovLoop()
    : m_vl(CommonParameters::Vlevel()), m_Nspc_size(0), m_Ntype(0)
  {
    init();
  }
 
  PolyakovLoop(const Parameters& params)
    : m_vl(CommonParameters::Vlevel()), m_Nspc_size(0), m_Ntype(0)
  {
    init();
    set_parameters(params);
  }
 
  virtual ~PolyakovLoop() {}
 
 private:
  // non-copyable
  PolyakovLoop(const PolyakovLoop&);
  PolyakovLoop& operator=(const PolyakovLoop&);
 
 public:
  virtual void set_parameters(const Parameters& params);
  void set_parameters(const int Nspc_size, const int Ntype);
 
  void set_parameter_verboselevel(const Bridge::VerboseLevel vl) { m_vl = vl; }
 
  virtual void get_parameters(Parameters& params) const;
 
  dcomplex measure_ploop(const Field_G& U);
 
  void calc_ploop(dcomplex& ploop, const Field_G& U);
 
  // double measure_ploop_corr(const Field_G& U);
 
 private:
 
  void init();
};
#endif