製作CURAND從均勻分佈生成不同的隨機數

Question

我試圖使用CURAND庫來生成隨機數，這是完全相互獨立的從0到100因此我給時間作爲種子給每個線程並指定「 id = threadIdx.x + blockDim.x * blockIdx.x「作爲序列和偏移量。 然後獲得隨機數作爲浮動後，我乘以100，並採取其整數值。現在

，我面臨的問題是，它得到相同的隨機數的線程[0,0]和[0,1]，無論多少次，我跑這是11我無法代碼瞭解我做錯了什麼。請幫忙。

我貼我下面的代碼：

#include <stdlib.h> 
#include <stdio.h> 
#include <math.h> 
#include<curand_kernel.h> 
#include "util/cuPrintf.cu" 
#include<time.h> 

#define NE WA*HA //Total number of random numbers 
#define WA 2 // Matrix A width 
#define HA 2 // Matrix A height 
#define SAMPLE 100 //Sample number 
#define BLOCK_SIZE 2 //Block size 

__global__ void setup_kernel (curandState * state, unsigned long seed) 
{ 
int id = threadIdx.x + blockIdx.x + blockDim.x; 
curand_init (seed, id , id, &state[id]); 
} 

__global__ void generate(curandState* globalState, float* randomMatrix) 
{ 
int ind = threadIdx.x + blockIdx.x * blockDim.x; 
if(ind < NE){ 
    curandState localState = globalState[ind]; 
    float stopId = curand_uniform(&localState) * SAMPLE; 
    cuPrintf("Float random value is : %f",stopId); 
    int stop = stopId ; 
    cuPrintf("Random number %d\n",stop); 
    for(int i = 0; i < SAMPLE; i++){ 
      if(i == stop){ 
        float random = curand_normal(&localState); 
        cuPrintf("Random Value %f\t",random); 
        randomMatrix[ind] = random; 
        break; 
      } 
    } 
    globalState[ind] = localState; 
} 
} 

///////////////////////////////////////////////////////// 
// Program main 
///////////////////////////////////////////////////////// 

int main(int argc, char** argv) 
{ 

// 1. allocate host memory for matrix A 
unsigned int size_A = WA * HA; 
unsigned int mem_size_A = sizeof(float) * size_A; 
float* h_A = (float*) malloc(mem_size_A); 
time_t t; 

// 2. allocate device memory 
float* d_A; 
cudaMalloc((void**) &d_A, mem_size_A); 

// 3. create random states  
curandState* devStates; 
cudaMalloc (&devStates, size_A*sizeof(curandState)); 

// 4. setup seeds 
int n_blocks = size_A/BLOCK_SIZE; 
time(&t); 
printf("\nTime is : %u\n",(unsigned long) t); 
setup_kernel <<< n_blocks, BLOCK_SIZE >>> (devStates, (unsigned long) t); 
// 4. generate random numbers 
cudaPrintfInit(); 
generate <<< n_blocks, BLOCK_SIZE >>> (devStates,d_A); 
cudaPrintfDisplay(stdout, true); 
cudaPrintfEnd(); 
// 5. copy result from device to host 
cudaMemcpy(h_A, d_A, mem_size_A, cudaMemcpyDeviceToHost); 


// 6. print out the results 
printf("\n\nMatrix A (Results)\n"); 
for(int i = 0; i < size_A; i++) 
{ 
    printf("%f ", h_A[i]); 
    if(((i + 1) % WA) == 0) 
     printf("\n"); 
} 
printf("\n"); 

// 7. clean up memory 
free(h_A); 
cudaFree(d_A); 

} 

輸出，我得到的是：

時間是：1347857063 [0,0]：浮法隨機值：11.675105 [0,0 ]：隨機數11 [0,0]：隨機值0.358356 [0，1]：浮點隨機值爲：11.675105 [0,1]：隨機數11 [0,1]：隨機值0.358356 [ 1，0]：浮點隨機值爲：63.840496 [1,0]：隨機數63 [1,0]：隨機值0.696459 [1,1]：浮動隨機值：44.712799 [1,1]：隨機數44 [1,1]：隨機值0.735049

Answer 1

有幾個事情錯在這裏，我在這裏解決的第一批，讓你開始：

一般分

• 請檢查所有CUDA API調用的返回值，請參閱here獲取更多信息。
• 請運行cuda-memcheck檢查顯而易見的事情，如超出邊界訪問。

具體點

• 當分配的RNG狀態空間，你應該有空間，每個線程一個狀態（你現在有每個矩陣元素不是一個）。
• 您在setup_kernel（）中的線程ID計算錯誤，應該是threadIdx.x + blockIdx.x * blockDim.x（*而不是+）。
• 您使用線程ID的序列號以及偏移量，您應該只設置在cuRAND手冊中描述的偏移爲零：

提供最高質量的並行僞隨機數生成，每個 實驗應分配一個獨特的種子。在一個實驗中，每個計算線程應該分配一個唯一的序列號碼 。

最後你每塊運行兩個線程，這是非常低效的。查看「CUDA C編程指南」中的「最大化利用率」部分以獲取更多信息，但您應該希望啓動每塊32個線程（例如128,256）和大量塊（例如數萬）。如果問題很小，那麼可以考慮一次運行多個問題（或者在單個內核啓動時進行批處理，或者作爲不同流中的內核進行批處理以獲得併發執行）。