astaple_lex.h

Go to the documentation of this file.

 
#ifndef ACCEL_ASTAPLE_LEX_INCLUDED
#define ACCEL_ASTAPLE_LEX_INCLUDED
 
// core library
#include "lib/IO/bridgeIO.h"
 
// alt-code
#include "lib_alt_Accel/Field/shiftAField_lex.h"
 
 
template<typename AFIELD>
class AStaple_lex
{
 public:
  typedef typename AFIELD::real_t real_t;
  static const std::string class_name;
 
 private:
  int m_Ndf, m_Nvol, m_Ndim;
 
  Bridge::VerboseLevel m_vl;
 
  ShiftAField_lex<AFIELD> *m_shift;
 
  AFIELD m_v1, m_v2, m_v3; 
  AFIELD m_Umu;            
 
  void init();
 
  void tidyup();
 
 public:
 
  AStaple_lex(){ init(); }
 
  ~AStaple_lex(){ tidyup(); }
 
  // void set_parameters(const Parameters& params);
 
  void plaquette(real_t& plaq, const AFIELD&);
 
  void plaq_s(real_t& plaq_s, const AFIELD&);
 
  void plaq_t(real_t& plaq_t, const AFIELD&);
 
  void staple(AFIELD&, const AFIELD&, const int mu);
 
  void upper(AFIELD&, const AFIELD&, const int mu, const int nu);
 
  void lower(AFIELD&, const AFIELD&, const int mu, const int nu);
 
  void shift_forward(AFIELD& v, const int ex1,
                     const AFIELD& w, const int ex2, const int mu)
  { m_shift->forward(v, ex1, w, ex2, mu); }
 
  void shift_backward(AFIELD& v, const int ex1,
                     const AFIELD& w, const int ex2, const int mu)
  { m_shift->backward(v, ex1, w, ex2, mu); }
 
};
#endif