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:
  Bridge::VerboseLevel m_vl;
 
  unsigned long int m_flop;
  double m_time;
 
 public:
  Projection_Stout_SU3()
    : m_vl(CommonParameters::Vlevel())
  {
    init();
  }
 
  Projection_Stout_SU3(const Parameters& params)
    : m_vl(CommonParameters::Vlevel())
  {
    init();
    set_parameters(params);
  }
 
  ~Projection_Stout_SU3() {}
 
  void set_parameters(const Parameters& param);
 
  void get_parameters(Parameters& param) const;
 
  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();
  }
 
  static Projection *create_object_with_params(const Parameters& params)
  {
    return new Projection_Stout_SU3(params);
  }
 
 public:
  static bool register_factory()
  {
    bool init = true;
    init &= Projection::Factory::Register("Stout_SU3", create_object);
    init &= Projection::Factory_params::Register("Stout_SU3", create_object_with_params);
    return init;
  }
#endif
};
#endif