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

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

  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