11#if (defined __cplusplus)
68 Complex *Phi,
double *dk4m,
double *dk4p,
float *dk4m_f,
float *dk4p_f,
double *dSdpi,
double *pp,
70 float beta,
float akappa,
int stepl,
float dt,
double *ancg,
int *itot,
float proby);
103 Complex *Phi,
double *dk4m,
double *dk4p,
float *dk4m_f,
float *dk4p_f,
double *dSdpi,
double *pp,
105 float beta,
float akappa,
int stepl,
float dt,
double *ancg,
int *itot,
float proby);
107 Complex *Phi,
double *dk4m,
double *dk4p,
float *dk4m_f,
float *dk4p_f,
double *dSdpi,
double *pp,
109 float beta,
float akappa,
int stepl,
float dt,
double *ancg,
int *itot,
float proby);
112 void cuGauge_Update(
const double d,
double *pp,
Complex *u11t,
Complex *u12t, dim3 dimGrid, dim3 dimBlock);
115#if (defined __cplusplus)
121__global__
void cuGauge_Update(
const double d,
double *pp,
Complex *u11t,
Complex *u12t,
int mu);
int OMF2(Complex *u11t, Complex *u12t, Complex_f *u11t_f, Complex_f *u12t_f, Complex *X0, Complex *X1, Complex *Phi, double *dk4m, double *dk4p, float *dk4m_f, float *dk4p_f, double *dSdpi, double *pp, int *iu, int *id, Complex *gamval, Complex_f *gamval_f, int *gamin, Complex jqq, float beta, float akappa, int stepl, float dt, double *ancg, int *itot, float proby)
OMF integrator. Each trajectory step takes the form of p->p+dt/2,u->u+dt,p->p+dt/2 In practice this i...
int OMF4(Complex *u11t, Complex *u12t, Complex_f *u11t_f, Complex_f *u12t_f, Complex *X0, Complex *X1, Complex *Phi, double *dk4m, double *dk4p, float *dk4m_f, float *dk4p_f, double *dSdpi, double *pp, int *iu, int *id, Complex *gamval, Complex_f *gamval_f, int *gamin, Complex jqq, float beta, float akappa, int stepl, float dt, double *ancg, int *itot, float proby)
int Gauge_Update(const double d, double *pp, Complex *u11t, Complex *u12t, Complex_f *u11t_f, Complex_f *u12t_f)
Gauge update for the integration step of the HMC.
int Leapfrog(Complex *u11t, Complex *u12t, Complex_f *u11t_f, Complex_f *u12t_f, Complex *X0, Complex *X1, Complex *Phi, double *dk4m, double *dk4p, float *dk4m_f, float *dk4p_f, double *dSdpi, double *pp, int *iu, int *id, Complex *gamval, Complex_f *gamval_f, int *gamin, Complex jqq, float beta, float akappa, int stepl, float dt, double *ancg, int *itot, float proby)
Leapfrog integrator. Each trajectory step takes the form of p->p+dt/2,u->u+dt,p->p+dt/2 In practice t...
int Momentum_Update(const double d, const double *dSdpi, double *pp)
Wrapper for the momentum update during the integration step of the HMC.
Header for random number configuration.
Defines the constants of the code and other parameters for loop dimensions. Each subroutine includes ...
#define Complex
Double precision complex number.
#define Complex_f
Single precision complex number.