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 "bridgeIO.h"
using Bridge::vout;


//- parameters class
class Parameters_Projection_Maximum_SU_N : virtual public Parameters
{
 public:
  Parameters_Projection_Maximum_SU_N();
};
//- end

class Projection_Maximum_SU_N : public Projection
{
 private:
  int    m_Niter;  
  double m_Enorm;  

 public:
  Projection_Maximum_SU_N()
  {
    //- defaults
    m_Niter = 100;
    m_Enorm = 1.0e-12;
  }

  ~Projection_Maximum_SU_N() {}

  void set_parameters(const Parameters& params);
  void set_parameters(const int Niter, const double Enorm);

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

  void force_recursive(Field_G& Xi, Field_G& iTheta,
                       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, Field_G& V);

  //- maximization by SU(2) subgroup.
  void maxTr_SU2(int, int, Field_G&, Field_G&, Field_G&);

  //- matrix index for convenience.
  int mindex(int i, int j, int Nc)
  {
    return i + j * Nc;
  }
};
#endif