topologicalCharge.h

Go to the documentation of this file.

 
#ifndef TOPOLOGICALCHARGE_INCLUDED
#define TOPOLOGICALCHARGE_INCLUDED
 
#include "fieldStrength.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class TopologicalCharge
{
 public:
  static const std::string class_name;
 
 protected:
  Bridge::VerboseLevel m_vl;
 
 private:
  std::string m_filename_output;
 
  double m_c_plaq;
  double m_c_rect;
 
  int m_max_mom;
 
  FieldStrength m_field_strength;
 
  std::vector<Field_G> m_Fmunu_plaq;
  std::vector<Field_G> m_Fmunu_1x1;
  std::vector<Field_G> m_Fmunu_1x2;
  bool m_flag_field_set;
 
 public:
  TopologicalCharge()
    : m_vl(CommonParameters::Vlevel())
  {
    m_filename_output = "stdout";
    m_flag_field_set  = false;
  }
 
  TopologicalCharge(const Parameters& params)
    : m_vl(CommonParameters::Vlevel())
  {
    m_filename_output = "stdout";
    m_flag_field_set  = false;
    set_parameters(params);
  }
 
  virtual ~TopologicalCharge() {}
 
 private:
  // non-copyable
  TopologicalCharge(const TopologicalCharge&);
  TopologicalCharge& operator=(const TopologicalCharge&);
 
 public:
  virtual void set_parameters(const Parameters& params);
  void set_parameters(const double c_plaq, const double c_rect, const int max_mom);
 
  void set_parameter_verboselevel(const Bridge::VerboseLevel vl) { m_vl = vl; }
 
  virtual void get_parameters(Parameters& params) const;
 
  double measure(const Field_G& U);
 
  void measure_topological_charge(const double tt);
 
  void measure_topological_density_t(const double tt);
 
  void measure_topological_density_t_FT(const double tt);
 
  void measure_topological_density_x(const double tt);
 
  void measure_topological_density_x_FT(const double tt);
 
  void measure_topological_density_y(const double tt);
 
  void measure_topological_density_y_FT(const double tt);
 
  void measure_topological_density_z(const double tt);
 
  void measure_topological_density_z_FT(const double tt);
 
  //  double contract_epsilon_tensor(const Field_G& Fmunu_1, const Field_G& Fmunu_2);
 
  void set_field_strength(const Field_G& U);
};
#endif