Centre Blaise Pascal » Bench4GPU

}

"""

def KernelCodeCuda():

    BlobCUDA= """

#else

#define TRACKPOINTS

#endif

__device__ float nothing(float x)

{

  return(x);

   uint sizex=(uint)gridDim.x*blockDim.x;

   uint sizey=(uint)gridDim.y*blockDim.y;

   // Perform trajectory for each pixel, exit on hit

  float m,rs,ri,re,tho;

  q=-2;

  raie=Line;

  float rp0,rpp,rps;

  int nh;

  float zp,fp;

  // Perform trajectory for each pixel

  int q,raie;

  m=Mass;

  ri=InternalRadius;

  re=ExternalRadius;

  tho=Angle;

  q=-2;

  raie=Line;

  float zp=0,fp=0;

  // Autosize for image, 25% greater than external radius

  bmx=1.25e0f*re;

   // Perform trajectory for each pixel

  int q,raie;

  m=Mass;

  ri=InternalRadius;

  re=ExternalRadius;

  tho=Angle;

  raie=Line;

  float bmx,db,b,h;

  float zp=0,fp=0;

  // Autosize for image

  // Perform trajectory for each pixel

  m=Mass;

  rs=2.e0f*m;

  re=ExternalRadius;

  int nh;

  // Autosize for image

  bmx=1.25e0f*re;

                    numpy.float32(ExternalRadius),

                    numpy.float32(Angle),

                    numpy.int32(Line),

                    block=(Threads,Threads,1))

    elif Method=='EachCircle':

        EachCircleCU(zImageCU,fImageCU,