alg_hmc.h

Go to the documentation of this file.

#include<config.h>
CPS_START_NAMESPACE
//------------------------------------------------------------------
//------------------------------------------------------------------


#ifndef INCLUDED_ALG_HMC_H
#define INCLUDED_ALG_HMC_H

CPS_END_NAMESPACE
#include <util/lattice.h>
#include <util/smalloc.h>
#include <util/pmalloc.h>
#include <alg/alg_base.h>
#include <alg/hmc_arg.h>
#include <alg/alg_meas.h>
#include <alg/cg_stats.h>
#include <util/checksum.h>
#include <alg/alg_int.h>
CPS_START_NAMESPACE

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

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

class AlgHmc
{
 private:
    char *cname;

    int saveInitialState();
    void saveFinalState();
    void restoreInitialState();
    void restoreFinalState();
    int reproduceTest(int attempt);
    void reverseTest();
    void shiftStates(const int x,const int y,const int z, const int t);

 protected:

    int g_size;
    Matrix* gauge_field_init;

    Matrix* gauge_field_final;

//    unsigned int **rng4d_init;
//    unsigned int **rng5d_init;
      LRGState lrg_state;

    int config_no;

    Float h_init;   
    Float h_final;  
    Float delta_h;  
    Float h_delta;  

    Matrix* mom;

    HmcArg *hmc_arg;
    CommonArg *common_arg;

    AlgIntAB *integrator;

    unsigned int checksum[2];

 public:

  AlgHmc(AlgIntAB &integrator, CommonArg &c_arg, HmcArg &arg);

  virtual ~AlgHmc();

  Float run(void);
};

#endif

CPS_END_NAMESPACE

---