CUDA Указатели, Визитные Ядра

Question

Если я использую указатели на функции C, как:

void processCalcNorm(float* a, float* b, float* c, float* d, float* e, float* f) 
    { 
      *a = *a + *b; 
      *c = *c + *d; 
      *e = *e + *f; 
    } 

for(id = 0; id < 1000; id++) 
    { 
      processCalcNorm(&xcord[id],&lvelox[id],&ycord[id],&lveloy[id],&zcord[id],&lveloz[id]); 
    } 

Как следует ядро ​​будет вызвана, когда она выполняется?

Answer 1

Нечто подобное должно работать (написано в браузере, не тестировалось):

__global__ void processCalcNorm_kernel(float* a, float* b, float* c, float* d, float* e, float* f, int len) 
    { 
     int idx = threadIdx.x + blockDim.x*blockIdx.x; 
     if (idx < len){ 
      a[idx] = a[idx] + b[idx]; 
      c[idx] = c[idx] + d[idx]; 
      e[idx] = e[idx] + f[idx];} 
    } 

#define DATA_LEN 1000 
#define nTPB 256 
... 
processCalcNorm_kernel<<<(DATA_LEN+nTPB-1)/nTPB, nTPB>>>(d_xcord,d_lvelox,d_ycord,d_lveloy,d_zcord,d_lveloz,DATA_LEN); 

Где d_... переменные являются копией устройства из одноименного переменного хоста, соответственно, установленная с cudaMalloc и cudaMemcpy звонков, как это (с использованием xcord в качестве примера):

float *d_xcord; 
cudaMalloc((void **)&d_xcord, DATA_LEN*sizeof(float)); 
cudaMemcpy(d_xcord, xcord, DATA_LEN*sizeof(float), cudaMemcpyHostToDevice); 

(и создать аналогичную последовательность для других переменных)

Обратите внимание, что цикл for в вашем исходном C-коде больше не требуется, поскольку графический процессор эффективно обрабатывает каждую итерацию вашего цикла for-loop одним вызовом ядра.