staple_SF.h

Go to the documentation of this file.

 
#ifndef STAPLE_SF_INCLUDED
#define STAPLE_SF_INCLUDED
 
#include "Field/field_SF.h"
#include "Field/shiftField_lex.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class Staple_SF
{
 public:
  static const std::string class_name;
 
 protected:
  Bridge::VerboseLevel m_vl;
 
 private:
  Index_lex m_index;
  ShiftField_lex m_shift;
 
  Mat_SU_N m_wk, m_wkpr, m_i_omega0;
  int m_initialized;
 
  std::vector<double> m_phi, m_phipr, m_p_omega;
 
 public:
  Staple_SF()
    : m_vl(CommonParameters::Vlevel()),
    m_wk(CommonParameters::Nc()),
    m_wkpr(CommonParameters::Nc()),
    m_i_omega0(CommonParameters::Nc()),
    m_initialized(0) {}
 
  // optional
  Staple_SF(const Parameters& params)
    : m_vl(CommonParameters::Vlevel()),
    m_wk(CommonParameters::Nc()),
    m_wkpr(CommonParameters::Nc()),
    m_i_omega0(CommonParameters::Nc()),
    m_initialized(0)
  {
    set_parameters(params);
  }
 
 private:
  // non-copyable
  Staple_SF(const Staple_SF&);
  Staple_SF& operator=(const Staple_SF&);
 
 public:
  void set_parameters(const Parameters& params);
 
  //  void set_parameters(const double *phi, const double *phipr);
  void set_parameters(const std::vector<double>& phi,
                      const std::vector<double>& phipr);
 
  void set_parameters(const double *phi, const double *phipr, const double *pomega);
  void set_parameters(const std::vector<double>& phi, const std::vector<double>& phipr, const std::vector<double>& p_omega);
 
  void set_parameter_verboselevel(const Bridge::VerboseLevel vl) { m_vl = vl; }
 
  void get_parameters(Parameters& params) const;
 
  void upper(Field_G&, const Field_G&, const int, const int);
  void lower(Field_G&, const Field_G&, const int, const int);
 
  double plaq_s(const Field_G&);
  double plaq_t(const Field_G&);
  double plaq_t_ct(const Field_G&, const double ct);
  double plaquette(const Field_G&);
  double plaquette_ct(const Field_G&, const double ct);
 
  double sf_coupling_plaq(const Field_G&, const double ct);
  double sf_coupling_rect(const Field_G&, const double ctr);
 
  void staple(Field_G&, const Field_G&, const int);
  void staple_ct(Field_G&, const Field_G&, const int, const double ct);
 
  void print_plaquette(const Field_G&);
};
#endif