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();
  }

  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; }

  dcomplex measure_ploop(const Field_G& U);

  // double measure_ploop_corr(const Field_G& U);

 private:
  void calc_ploop(dcomplex& ploop, const Field_G& U);

  void init();
};
#endif