projection_Maximum_SU_N.h

Go to the documentation of this file.

 
#ifndef PROJECTION_MAXIMUM_SU_N_INCLUDED
#define PROJECTION_MAXIMUM_SU_N_INCLUDED
 
#include "projection.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class Projection_Maximum_SU_N : public Projection
{
 public:
  static const std::string class_name;
 
 private:
  Bridge::VerboseLevel m_vl;
 
  int m_Niter;     
  double m_Enorm;  
 
 public:
  Projection_Maximum_SU_N()
    : m_vl(CommonParameters::Vlevel())
  {
    //- defaults
    m_Niter = 100;
    m_Enorm = 1.0e-12;
  }
 
  Projection_Maximum_SU_N(const Parameters& params)
    : m_vl(CommonParameters::Vlevel())
  {
    //- defaults
    m_Niter = 100;
    m_Enorm = 1.0e-12;
 
    set_parameters(params);
  }
 
  ~Projection_Maximum_SU_N() {}
 
  void set_parameters(const Parameters& params);
  void set_parameters(const int Niter, const double Enorm);
 
  void get_parameters(Parameters& params) 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);
 
 private:
  //- maximization of ReTr[U^\dag V].
  void maxTr(Field_G& U, const Field_G& V);
 
  //- maximization by SU(2) subgroup.
  void maxTr_SU2(const int, const int, Field_G&, Field_G&, Field_G&);
 
  //- matrix index for convenience.
  int mindex(const int i, const int j, const int Nc)
  {
    return i + j * Nc;
  }
 
#ifdef USE_FACTORY
 private:
  static Projection *create_object()
  {
    return new Projection_Maximum_SU_N();
  }
 
  static Projection *create_object_with_params(const Parameters& params)
  {
    return new Projection_Maximum_SU_N(params);
  }
 
 public:
  static bool register_factory()
  {
    bool init = true;
    init &= Projection::Factory::Register("Maximum_SU_N", create_object);
    init &= Projection::Factory_params::Register("Maximum_SU_N", create_object_with_params);
    return init;
  }
#endif
};
#endif