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(const Parameters& params)
    : m_vl(CommonParameters::Vlevel())
  {
    m_filename_output = "stdout";
    m_flag_field_set  = 0;
 
    set_parameters(params);
  }
 
  ~EnergyMomentumTensor() {}
 
 public:
  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; }
 
  void get_parameters(Parameters& params) const;
 
  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