下面是这段代码是我的一个算法中用来求和以及求平方和的kernel函数:
__global__ static void CompSumAndSquare(int rate,int I_n,int size,int width,int wsize,int *image,float *sum,float *sumOfSquare)
{
int x,y,p,i,j,ix,k;
const int pos = threadIdx.x+blockIdx.x*blockDim.x;
const int max_dimension=blockDim.x*gridDim.x;
k=wsize/2;
for(i=0;i<rate;i++)
{
p=pos+i*max_dimension;
if(p<size){
x=p%I_n+wsize/2;
y=p/I_n;
ix=y*width+x;
sum[ix]=0;
sumOfSquare[ix]=0;
for(j=-k;j<=k;j++)