alg_rnd_gauge.h

Go to the documentation of this file.

//----------------------------------------------------------------------

#ifndef __ROTATEGAUGE__CD
#define __ROTATEGAUGE__CD

#include <util/lattice.h>
#include <alg/alg_base.h>
#include <alg/common_arg.h>
#include <util/vector.h>

CPS_START_NAMESPACE




class AlgRandomGauge: public Alg
{
  private:

    char *cname;

    Float theta;

    // parint out matrix (for debugging)
    void printMatrix( const Matrix& x );
  
    void RandMatrix( Float theta , Matrix& x , int sub );

    void UpperRandMatrix( Float theta , Matrix& x ) { 
        RandMatrix( theta, x, 0 ); 
    }

    void LowerRandMatrix( Float theta , Matrix& x ) { 
        RandMatrix( theta, x, 1 ); 
    }

    void MixedRandMatrix( Float theta , Matrix& x ) { 
        RandMatrix( theta, x, 2 ); 
    }

    int NX, NY, NZ, NT;
    int index(int, int, int, int);
        
  public:

    AlgRandomGauge ( Lattice& latt , CommonArg *c_arg );


    ~AlgRandomGauge() {;}


    void set_theta( Float t ) { theta=t; } 
  
    void run();


    void free() { AlgLattice().FixGaugeFree(); }
};



class AlgRotateGauge : public Alg
{
private:

    char  *cname;
    int NX, NY, NZ, NT;
    int index(int, int, int, int, int);    
  
public:

  AlgRotateGauge( Lattice& latt, CommonArg *c_arg ):
    Alg(latt,c_arg),
    cname("RotateGauge")
  {;}
  
  ~AlgRotateGauge() {;}

  // Apply the gauge transformation.
  void run();
    

private:

  Matrix GetMat(const Matrix *m_offset, const int *x, int nu,
                const int *node_sites, const int *m_dir_offset);

};

CPS_END_NAMESPACE
#endif

---