staple_SF.h

Go to the documentation of this file.

 
#ifndef STAPLE_SF_INCLUDED
#define STAPLE_SF_INCLUDED
 
#include "lib/Field/field_SF.h"
#include "lib/Field/shiftField_lex.h"
#include "lib/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;
 
  Field_G *m_Umu;
  Field_G *m_Unu;
  Field_G *m_vt;
  Field_G *m_wt;
 
  Field_G *m_stpl_upper;
  Field_G *m_stpl_lower;
 
 public:
 
  Staple_SF()
    : m_wk(CommonParameters::Nc()),
      m_wkpr(CommonParameters::Nc()),
      m_i_omega0(CommonParameters::Nc())
  { init(); }
 
  // optional
  Staple_SF(const Parameters& params)
    : m_vl(CommonParameters::Vlevel()),
    m_wk(CommonParameters::Nc()),
    m_wkpr(CommonParameters::Nc()),
    m_i_omega0(CommonParameters::Nc())
  { init(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&);
 
 private:
 
  void init(const Parameters& params);
 
  void init();
 
  void tidyup();
 
 
};
#endif