我在cuda的__syncthreads（）後丟失了數據

-1

我試圖找到一個數組的最大值..我從CUDA Maximum Reduction Algorithm Not Working獲得了幫助。並做一些自己的修改。然而，我正在運行它的16個數據。我發現在內核代碼共享內存中只複製第4個數據。休息都會丟失。我把兩個cuPrintf..1st printf顯示的數據放在共享內存中。但第二cuPrintf只是__syncthreads後..並表示從0線程ID 4 onwords ..請幫助的#include 的#include 的#include 的#include 的#include 的#include 「cuPrintf.cu」的#include「 cuPrintf.cuh「我在cuda的__syncthreads（）後丟失了數據

__device__ float MaxOf2(float a, float b) 
{ 
    if(a > b) return a; 
    else   return b; 
} 

__global__ void findMax(int size,float *array_device , float *outPut) 
{ 
    extern __shared__ float sdata[]; 
    int tid = threadIdx.x; 
    int i = blockIdx.x*blockDim.x + threadIdx.x; 

    if(i< size) 
    { 
    sdata[tid] = array_device[i]; 
    cuPrintf(" array_d[%d]===%f, sdata[%d]===%f\n ",i,array_device[i],tid,sdata[tid]); 
    __threadfence(); 

    } 
    __syncthreads(); 

    if(tid<size) 
     cuPrintf(" array_d[%d]===%f, sdata[%d]===%f\n ",i,array_device[i],tid,sdata[tid]); 

    for (int s=blockDim.x/2; s>0; s=s>>1)//s=blockDim.x/2 
    { 
     if (tid < s) 
     { 
      sdata[tid]= MaxOf2(sdata[tid],sdata[tid+s]); 

     } 

     __syncthreads(); 
    } 
    if (tid == 0) outPut[blockIdx.x] = sdata[0]; 

} 

int main() 
{ 
    long double M = pow(2,20); 
    long double N = 2; 
    int noThreadsPerBlock = 512 ; 
    printf("\n Provide the array Size N.(array will be of size N * 2^20) :-"); 
    scanf("%Lf",&N); 
    long int size = 16; 
    int numOfBlock = (int)size /noThreadsPerBlock + 1; 
    printf("\n num of blocks==%ld",numOfBlock); 

    float *array_device , *outPut; 
    float array_host[]={221,100,2,340,47,36,500,1,33,4460,5,6,7,8,9,11}; 
    cudaMalloc((void **)&array_device, size*sizeof(float)); 
    cudaMalloc((void **)&outPut, size*sizeof(float)); 
    cudaError_t error0 = cudaGetLastError(); 
    printf("\n 0CUDA error: %s\n", cudaGetErrorString(error0)); 
    printf("size===%ld",size); 

    cudaMemcpy(array_device, array_host, size*sizeof(float), cudaMemcpyHostToDevice); 
    cudaError_t error1 = cudaGetLastError(); 
    printf("\n1CUDA error: %s\n", cudaGetErrorString(error1)); 


    while(size>1) 
    { 
     cudaPrintfInit(); 
     findMax<<< numOfBlock,noThreadsPerBlock>>>(size,array_device, outPut);cudaPrintfDisplay(stdout, true); 
     cudaPrintfEnd(); 
     cudaError_t error2 = cudaGetLastError(); 
     printf(" 2CUDA error: %s\n", cudaGetErrorString(error2)); 
     cudaMemcpy(array_device, outPut, size*sizeof(float), cudaMemcpyDeviceToDevice); 
     size = numOfBlock; 
     printf("\n ****size==%ld\n",size); 
     numOfBlock = (int)size /noThreadsPerBlock + 1; 
    } 

    cudaMemcpy(array_host, outPut, size*sizeof(float), cudaMemcpyDeviceToHost); 
    cudaError_t error3 = cudaGetLastError(); 
    printf("\n3CUDA error: %s\n", cudaGetErrorString(error3)); 
    for(int i=0;i<size;i++) 
      printf("\n index==%d ;data=%f ",i,array_host[i]); 
    return 0; 
}

來源

2013-08-28 Abhidip Bhattacharyya

首先，你還沒有在內核啓動時指定動態大小的共享內存。它應該如下所示：'findMax <<< numOfBlock，noThreadsPerBlock，sizeof（float）* noThreadsPerBlock >>>' 其次，第二個cuPrintf上if（tid Maku

是否真的有必要發佈很多註釋掉的代碼？您要求我們閱讀並更正您的代碼，至少花時間正確格式化，以便我們閱讀很容易。幫助我們幫助你..... – talonmies

我說如果條件只是爲了禁止這麼多的打印，我只有16個數據。我試圖把輸出..但geting一些問題來保存它，因爲它是「一些代碼在那裏......」。 –

我按照要求發佈我的評論作爲答案。

首先，你沒有在內核啓動時指定動態大小的共享內存。它應該是這個樣子： findMax<<< numOfBlock,noThreadsPerBlock,sizeof(float)*noThreadsPerBlock>>>

其次，什麼是落後於第二cuPrintf條件if(tid<size)概念？提供程序輸出也可能有所幫助。

來源

2013-08-30 12:50:11 Maku

感謝您的回答。它有助於我們將這些廢棄的問題從無人問候的清單中解脫出來。 – talonmies

我在cuda的__syncthreads（）後丟失了數據

回答

相關問題