alg_smear2.h

Go to the documentation of this file.

//------------------------------------------------------------------
//------------------------------------------------------------------
#ifndef _ALG_SMEAR_CD_
#define _ALG_SMEAR_CD_
#include <config.h>
#include <alg/alg_base.h>
#include <alg/ape_smear_arg.h>
//#include <alg/kinetic_smear_arg.h>
//#include <alg/hyp_smear_arg.h>


CPS_START_NAMESPACE


class AlgSmear2:public Alg
{
private:
  char *cname;

  int bool_su3_proj;
  Float tolerance  ;
  int orthog       ;
  Matrix* lat_back ;
    
public:
  
  AlgSmear2( Lattice&   lat,
            CommonArg* ca ,
            int su3_proj );
  
  virtual ~AlgSmear2();
  
  void run();
  

  void set_orthog( int i ) { orthog=i; }
  

  int  get_orthog() const  { return orthog; }

  void set_tol( Float x ) { tolerance = x; }
  
  Float get_tol() const { return tolerance; }

protected:
  

  virtual void smear_link2( Matrix& link, 
                           int*    pos,
                           int      mu )=0;
  
  
};

void three_staple2(Lattice& latt, Matrix& link, int *pos, int u, int orth);

//void five_staple  (Lattice& latt, Matrix& link, int *pos, int u, int orth);

//void seven_staple (Lattice& latt, Matrix& link, int *pos, int u, int orth);

//void lepage_staple(Lattice& latt, Matrix& link, int *pos, int u, int orth);
  
int su3_proj( Matrix& x , Float tolerance );



class AlgApeSmear2:public AlgSmear2
{
private:
  
  Float c;
  char *cname;
  
public:

  AlgApeSmear2(Lattice&     lat,
              CommonArg*   ca ,
              ApeSmearArg* asa );
  
  ~AlgApeSmear2()
  {;}
  
  void run();

  void set_coef( Float x ) { c = x ; }

protected:
  
  void smear_link2( Matrix& mat,
                   int*    pos,
                   int      mu );
};


CPS_END_NAMESPACE
#endif /* _ALG_SMEAR_CD_ */

---