12#if (defined __cplusplus)
76 double *dSdpi,
double *pp,
unsigned int *iu,
unsigned int *
id,
Complex gamval[20],
Complex_f gamval_f[20],
const unsigned short gamin[16],
77 Complex *sigval,
Complex_f *sigval_f,
unsigned short *sigin,
const Complex jqq,
const float beta,
const float akappa,
78 const float c_sw,
const int stepl,
const float dt,
double *ancg,
int *itot,
const float proby);
112 double *dSdpi,
double *pp,
unsigned int *iu,
unsigned int *
id,
Complex gamval[20],
Complex_f gamval_f[20],
const unsigned short gamin[16],
113 Complex *sigval,
Complex_f *sigval_f,
unsigned short *sigin,
const Complex jqq,
const float beta,
const float akappa,
114 const float c_sw,
const int stepl,
const float dt,
double *ancg,
int *itot,
const float proby);
148 double *dSdpi,
double *pp,
unsigned int *iu,
unsigned int *
id,
Complex gamval[20],
Complex_f gamval_f[20],
const unsigned short gamin[16],
149 Complex *sigval,
Complex_f *sigval_f,
unsigned short *sigin,
const Complex jqq,
const float beta,
const float akappa,
150 const float c_sw,
const int stepl,
const float dt,
double *ancg,
int *itot,
const float proby);
167#if (defined __cplusplus)
void cuGauge_Update(const double d, double *pp, Complex *ut[2], dim3 dimGrid, dim3 dimBlock)
CUDA wrapper for the gauge update during the integration step of the HMC.
int OMF2(Complex *ut[2], Complex_f *ut_f[2], Complex *X0, Complex *X1, Complex *Phi, double *dk[2], float *dk_f[2], double *dSdpi, double *pp, unsigned int *iu, unsigned int *id, Complex gamval[20], Complex_f gamval_f[20], const unsigned short gamin[16], Complex *sigval, Complex_f *sigval_f, unsigned short *sigin, const Complex jqq, const float beta, const float akappa, const float c_sw, const int stepl, const float dt, double *ancg, int *itot, const float proby)
OMF second order five step integrator.
int Leapfrog(Complex *ut[2], Complex_f *ut_f[2], Complex *X0, Complex *X1, Complex *Phi, double *dk[2], float *dk_f[2], double *dSdpi, double *pp, unsigned int *iu, unsigned int *id, Complex gamval[20], Complex_f gamval_f[20], const unsigned short gamin[16], Complex *sigval, Complex_f *sigval_f, unsigned short *sigin, const Complex jqq, const float beta, const float akappa, const float c_sw, const int stepl, const float dt, double *ancg, int *itot, const 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 OMF4(Complex *ut[2], Complex_f *ut_f[2], Complex *X0, Complex *X1, Complex *Phi, double *dk[2], float *dk_f[2], double *dSdpi, double *pp, unsigned int *iu, unsigned int *id, Complex gamval[20], Complex_f gamval_f[20], const unsigned short gamin[16], Complex *sigval, Complex_f *sigval_f, unsigned short *sigin, const Complex jqq, const float beta, const float akappa, const float c_sw, const int stepl, const float dt, double *ancg, int *itot, const float proby)
OMF fourth order eleven step integrator.
int Gauge_Update(const double d, double *pp, Complex *ut[2], Complex_f *ut_f[2])
Gauge update for the integration step of the HMC.
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.
#define Complex
Double precision complex number.
#define Complex_f
Single precision complex number.
dim3 dimGrid
Default grid size. First component is normally nt. Second and third depend whatever is needed to get ...
dim3 dimBlock
Default block size. Usually 128.