polyakovLoop.h

Go to the documentation of this file.

#ifndef POLYAKOVLOOP_INCLUDED
#define POLYAKOVLOOP_INCLUDED

#include <cassert>
#include "parameters.h"
#include "field_G.h"
#include "shiftField_eo.h"

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


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

class PolyakovLoop
{
 protected:
  Bridge::VerboseLevel m_vl;

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

 public:

  PolyakovLoop()
    : m_vl(CommonParameters::Vlevel()), m_Nspc_size(0), m_Ntype(0)
  {
    init();
  }

  virtual ~PolyakovLoop() {}

  virtual void set_parameters(const Parameters& params);

  void set_parameters(int Nspc_size, int Ntype);

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

  double measure_ploop(Field_G& U);

  double measure_ploop_corr(Field_G& U);

 private:

  void calc_ploop(Field_G& P, Field_G& U);

  void init();
};
#endif