energyMomentumTensor.h

Go to the documentation of this file.

#ifndef ENERGY_MOMENTUM_TENSOR_INCLUDED
#define ENERGY_MOMENTUM_TENSOR_INCLUDED

#include "fieldStrength.h"

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


class EnergyMomentumTensor
{
 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;
  int m_flag_field_set;

 public:
  EnergyMomentumTensor()
    : m_vl(CommonParameters::Vlevel())
  {
    m_filename_output = "stdout";
    m_flag_field_set  = 0;
  }

  ~EnergyMomentumTensor() {}

 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_EMT(const double t_flow);

  double measure_EMT_at_t(const double t_flow);

  double measure_EMT_at_t_FT(const double t_flow);

  double measure_EMT_at_x(const double t_flow);

  double measure_EMT_at_x_FT(const double t_flow);

  double measure_EMT_at_y(const double t_flow);

  double measure_EMT_at_y_FT(const double t_flow);

  double measure_EMT_at_z(const double t_flow);

  double measure_EMT_at_z_FT(const double t_flow);

  void set_field_strength(const Field_G& U);

 private:
  int factor(const int mu, const int nu);

  int index_munu2i(int mu, int nu);
};
#endif