alg_hmd.h

Go to the documentation of this file.

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


#ifndef INCLUDED_ALG_HMD_H
#define INCLUDED_ALG_HMD_H

CPS_END_NAMESPACE
#include <util/lattice.h>
#include <util/smalloc.h>
#include <util/pmalloc.h>
#include <alg/alg_base.h>
#include <alg/common_arg.h>
#include <alg/hmd_arg.h>
CPS_START_NAMESPACE


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

//------------------------------------------------------------------
class AlgHmd : public Alg
{
 private:
    char *cname;
   
 protected:

    HmdArg *hmd_arg;
  
    int g_size;       

    int Ncb;

    Matrix* mom;

    int alloc_sum;

 public:

  AlgHmd(Lattice& latt, CommonArg *c_arg, HmdArg *arg);


  virtual ~AlgHmd();

  virtual Float run() = 0;

};


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

//------------------------------------------------------------------
class AlgHmcPhi : public AlgHmd
{
 private:
    char *cname;


 protected:

    int n_frm_masses;     

    int n_bsn_masses;     

    int f_size;       

    CgArg **frm_cg_arg;

    CgArg **bsn_cg_arg;

    Vector** phi;

    Vector** bsn;

    Matrix* gauge_field_init;

    Vector** frm1;
    Vector** frm2;

    Vector** cg_sol_cur;

    Vector** cg_sol_prev;

    Float *h_f_init;    

    Float *h_f_final;   

    Float *delta_h_f;   

    Float *h_b_init;    

    Float *h_b_final;   

    Float *delta_h_b;   

 public:

  AlgHmcPhi(Lattice& latt, CommonArg *c_arg, HmdArg *arg);

  virtual ~AlgHmcPhi();

  Float run(void);
};

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

//------------------------------------------------------------------
class AlgHmcQPQ : public AlgHmd
{
 private:
    char *cname;


 protected:

    int n_frm_masses;     

    int n_bsn_masses;     

    int f_size;       

    CgArg **frm_cg_arg;

    CgArg **bsn_cg_arg;

    Vector** phi;

    Vector** bsn;

    Matrix* gauge_field_init;

    Vector** frm1;
    Vector** frm2;

    Vector** cg_sol_prev;
    Vector*** cg_sol;
    Vector*** vm;

    Float *h_f_init;    

    Float *h_f_final;   

    Float *delta_h_f;   

    Float *h_b_init;    

    Float *h_b_final;   

    Float *delta_h_b;   

 public:

  AlgHmcQPQ(Lattice& latt, CommonArg *c_arg, HmdArg *arg);

  virtual ~AlgHmcQPQ();

  Float run(void);
};

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

//------------------------------------------------------------------
class AlgHmcRHMC : public AlgHmd
{
 private:
    char *cname;

    void init();

    void massRenormalise(Float *mass, Float *trueMass, int degree, 
                         Float *shift, MassRenormaliseDir direction);

 protected:

    int n_frm_masses;     

    int n_bsn_masses;     

    int f_sites;       
    int f_vec_count;       
    int f_count;       
    int f_size;       

    CgArg **frm_cg_arg;

    CgArg **bsn_cg_arg;

    Vector** phi;

    Vector** bsn;

    Matrix* gauge_field_init;

    unsigned int **rng4d_init;
    unsigned int **rng5d_init;

    Matrix* gauge_field_final;

    Vector* frm;

    Vector** frmn;

    Vector** frmn_d;

    Float *h_f_init;    

    Float *h_f_final;   

    Float *delta_h_f;   

    Float *h_b_init;    

    Float *h_b_final;   

    Float *delta_h_b;   

    int total_size;

    Float *all_res;

    EigArg *eig_arg;

    Float **alpha;

 public:

  AlgHmcRHMC(Lattice& latt, CommonArg *c_arg, HmdArg *arg);

  AlgHmcRHMC(Lattice& latt, CommonArg *c_arg, HmdArg *arg, EigArg *e_arg);

  virtual ~AlgHmcRHMC();

  Float run(void);

  void generateApprox(HmdArg*);

  void dynamicalApprox();

};

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

//------------------------------------------------------------------
class AlgHmdR : public AlgHmd
{
 private:
    char *cname;

 protected:
    int n_frm_masses;     

    Float *flavor_time_step;

    int f_size;       

    CgArg **frm_cg_arg;

    Vector** phi;

    Vector* frm1;
    Vector* frm2;
    
 public:

  AlgHmdR(Lattice& latt, CommonArg *c_arg, HmdArg *arg);

  virtual ~AlgHmdR();

  Float run(void);
};

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

//------------------------------------------------------------------
class AlgHmdR2 : public AlgHmd
{
 private:
    char *cname;

 protected:
    int n_frm_masses;     

    Float *flavor_time_step;

    Float *force_coeff;

    int f_sites;       
    int f_vec_count;       
    int f_count;       
    int f_size;       

    CgArg **frm_cg_arg;

    Vector** phi;

    Vector** frmn;

    Vector** frmn_d;

    Float *shift;

    int light;

    int heavy;

 public:

  AlgHmdR2(Lattice& latt, CommonArg *c_arg, HmdArg *arg);

  virtual ~AlgHmdR2();

  Float run(void);
};

#endif

CPS_END_NAMESPACE

---