clover.h

Go to the documentation of this file.

#pragma once
#include <su2hmc.h>

typedef struct{
   float *diag;
   Complex_f *offd;
}Bilinear_a;

typedef struct{
   float diag[2];
   Complex_f offd;
}Bilinear;

void ByGenLeft(Complex_f a[nc],const unsigned short gen);
void ByGenRight(Complex_f a[nc],const unsigned short gen);

void Clover_SU2plaq(Complex_f *ut[nc], Complex_f Leaves[nc], unsigned int *iu,  int i, int mu, int nu);
void Half_Leaves(Complex_f *hLeaves[2],Complex_f *ut[2], unsigned int *iu,unsigned int *id,\
      const unsigned short mu,const unsigned short nu);
void Leaf(Complex_f Leaves[nc],Complex_f *ut[nc], unsigned int *iu, unsigned int *id, unsigned int i,\
      const unsigned short mu, const unsigned short nu,const unsigned short leaf);
void Clover(Complex_f *clover[2], Complex_f *ut[2], unsigned int *iu, unsigned int *id);
void ByClover(Complex *phi, Complex *r, Complex *clover[2], Complex *sigval, const float akappa, unsigned short *sigin, bool dag);
void ByClover_f(Complex_f *phi, Complex_f *r, Complex_f *clover[2], Complex_f *sigval, const float akappa, unsigned short *sigin, bool dag);
void HbyClover(Complex *phi, Complex *r, Complex *clover[2],Complex *sigval, const float akappa, unsigned short *sigin,bool dag);
void HbyClover_f(Complex_f *phi, Complex_f *r, Complex_f *clover[2],Complex_f *sigval, const float akappa, unsigned short *sigin,bool dag);

void CalcXmunu(Bilinear_a Xmunu, Complex_f *X1, Complex_f *X2, const Complex_f *sigval,\
               const unsigned short *sigin,const unsigned short mu, const unsigned short nu);
void Clov_Force(double *dSdpi, Complex_f *ut[2], Complex_f *X1, Complex_f *X2, const Complex_f *sigval, const unsigned short *sigin,\
                  unsigned int *iu, unsigned int *id, const float akappa);
int Init_clover(Complex **sigval, Complex_f **sigval_f,unsigned short **sigin, float c_sw);
void Clover_free(Complex_f *clover[nc]);
 
#ifdef USE_GPU
#ifdef __cplusplus
extern "C"
{
#endif

int cuClover(Complex_f *clover[nc],Complex_f *ut[nc], unsigned int *iu, unsigned int *id);
void cuByClover(Complex *phi, Complex *r, Complex *clover[nc],Complex *sigval, const float akappa, unsigned short *sigin,bool dag);
void cuHbyClover(Complex *phi, Complex *r, Complex *clover[nc],Complex *sigval, const float akappa, unsigned short *sigin,bool dag);
void cuByClover_f(Complex_f *phi, Complex_f *r, Complex_f *clover[nc],Complex_f *sigval, const float akappa,unsigned short *sigin,bool dag);
void cuHbyClover_f(Complex_f *phi, Complex_f *r, Complex_f *clover[nc],Complex_f *sigval, const float akappa,unsigned short *sigin,bool dag);
void cuCalcXmunu(Bilinear_a Xmunu, Complex_f *X1, Complex_f *X2, const Complex_f *sigval,\
      const unsigned short *sigin,const unsigned short mu, const unsigned short nu);
int cuClov_Force(double *dSdpi, Complex_f *ut[nc], Complex_f *X1, Complex_f *X2, const Complex_f *sigval,\
      const unsigned short *sigin, const unsigned int *iu, const unsigned int *id, const float akappa);
#ifdef __cplusplus
}
#endif
#endif