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:
  int m_Niter;         
  double m_Enorm;      
 
  Bridge::VerboseLevel m_vl; 
 
  Field_G *m_At;
  Field_G *m_Dlt;
  Field_G *m_Ut;
  Field_G *m_vt;
  Field_G *m_wt;
 
 public:
 
  Projection_Maximum_SU_N(const Parameters& params)
  { init(params); }
 
  ~Projection_Maximum_SU_N() {tidyup(); }
 
  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:
 
  void init(const Parameters& params);
 
  void tidyup();
 
  void reset_field(const Field_G& Uref);
 
  void maxTr(Field_G& U, const Field_G& V);
 
  void maxTr_SU2(const int, const int, Field_G&, Field_G&, Field_G&);
 
  int mindex(const int i, const int j, const int Nc)
  { return i + j * Nc; }
 
  
#ifdef USE_FACTORY
 private:
 
  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_params::Register("Maximum_SU_N", create_object_with_params);
    return init;
  }
#endif
};
#endif