fieldStrength.h

Go to the documentation of this file.

#ifndef FIELDSTRENGTH_INCLUDED
#define FIELDSTRENGTH_INCLUDED

#include "staple_lex.h"

#include "Field/shiftField_lex.h"

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


class FieldStrength
{
 public:
  static const std::string class_name;

 protected:
  Bridge::VerboseLevel m_vl;

 private:
  ShiftField_lex m_shift;
  Staple_lex m_staple;

 public:
  FieldStrength()
    : m_vl(CommonParameters::Vlevel()) {}

  virtual ~FieldStrength() {}

 private:
  // non-copyable
  FieldStrength(const FieldStrength&);
  FieldStrength& operator=(const FieldStrength&);

 public:
  void construct_Fmunu_1x1(Field_G& Fmunu,
                           const int mu, const int nu, const Field_G& U);

  void construct_Fmunu_1x2(Field_G& Fmunu,
                           const int mu, const int nu, const Field_G& U);

  void construct_Fmunu_1x1_traceless(Field_G& Fmunu,
                                     const int mu, const int nu, const Field_G& U);

  void construct_Fmunu_1x2_traceless(Field_G& Fmunu,
                                     const int mu, const int nu, const Field_G& U);

  void construct_Fmunu_plaq_traceless(Field_G& Fmunu,
                                      const int mu, const int nu, const Field_G& U);

  double contract(const Field_G& Fmunu_1, const Field_G& Fmunu_2);

  void contract_at_t(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2);

  void contract_at_t(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2,
                     const std::vector<int>& momentum_sink,
                     const std::vector<int>& source_position);

  void contract_at_x(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2);

  void contract_at_x(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2,
                     const std::vector<int>& momentum_sink,
                     const std::vector<int>& source_position);

  void contract_at_y(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2);

  void contract_at_y(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2,
                     const std::vector<int>& momentum_sink,
                     const std::vector<int>& source_position);

  void contract_at_z(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2);

  void contract_at_z(std::vector<double>& corr_global,
                     const Field_G& Fmunu_1, const Field_G& Fmunu_2,
                     const std::vector<int>& momentum_sink,
                     const std::vector<int>& source_position);

  void global_corr_x(std::vector<double>& corr_global,
                     const std::vector<double>& corr_local);

  void global_corr_y(std::vector<double>& corr_global,
                     const std::vector<double>& corr_local);

  void global_corr_z(std::vector<double>& corr_global,
                     const std::vector<double>& corr_local);

  void global_corr_t(std::vector<double>& corr_global,
                     const std::vector<double>& corr_local);
};

#endif