projection_Stout_SU3.h

Go to the documentation of this file.

#ifndef PROJECTION_STOUT_SU3_INCLUDED
#define PROJECTION_STOUT_SU3_INCLUDED

#include "projection.h"

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

// The following implementation only valid for Nc = 3 case.
#define  NC    3



class Projection_Stout_SU3 : public Projection
{
 public:
  static const std::string class_name;

 private:
  unsigned long int m_flop;
  double m_time;

 public:

  Projection_Stout_SU3()
  {
    init();
  }

  ~Projection_Stout_SU3() {}

  void set_parameters(const Parameters& param);

  void project(Field_G& U,
               const double alpha,
               const Field_G& C, const Field_G& Uorg);

  void force_recursive(Field_G& Xi, Field_G& iTheta,
                       const double alpha, const Field_G& Sigmap,
                       const Field_G& C, const Field_G& U);

  void print_stat();

 private:
  void exp_iQ(Field_G& e_iQ, const Field_G& iQ);
  void exp_iQ_bf(Field_G& e_iQ, const Field_G& iQ);

  void set_uw(double& u, double& w,
              const Mat_SU_N& iQ2, const Mat_SU_N& iQ3);

  void set_fj(dcomplex& f0, dcomplex& f1, dcomplex& f2,
              const double& u, const double& w);

  double func_xi0(const double w);
  double func_xi1(const double w);

  void init();

#ifdef USE_FACTORY
 private:
  static Projection *create_object()
  {
    return new Projection_Stout_SU3();
  }

 public:
  static bool register_factory()
  {
    return Projection::Factory::Register("Stout_SU3", create_object);
  }
#endif
};
#endif