在cuda中實現隨機發生器

Question

我想在cuda中實現這個功能作爲設備/全局函數，以便獲得高斯分佈的隨機數。

double gasdev2() { 
double ran3n(long *seed); 
// double genrand64_real3(); 
static int iset=0; 
static double gcos; 
double tmp1,tmp2; 

if (iset==0) { 
tmp1=sqrt(-2*log(ran3n(&seed))); 
tmp2=pi2*ran3n(&seed); 
// tmp1=sqrt(-2*log(genrand64_real3())); 
// tmp2=pi2*genrand64_real3(); 
    gcos=tmp1*cos(tmp2); 
iset=1; 
return tmp1*sin(tmp2); 
//return 1; 
}else{ 
    iset=0; 
    return gcos; 
//return 1; 
} 
} 

這個功能將在這些函數調用和串行代碼基本上使用這些都是這樣

for(int i=0;i<NTO;i++){ 
    Frdx[j]=gasdev2()*ranm[j]*tconst; 
    Frdy[j]=gasdev2()*ranm[j]*tconst; 
    Frdz[j]=gasdev2()*ranm[j]*tconst; 
    } 

Answer 1

我建議不要自己實現它，但使用Thrust提供的隨機算法：

uint32_t seed = 1234; 
thrust::default_random_engine rng(seed); 
thrust::uniform_real_distribution<float> dist(0.0f, 1.0f); 
float random_value_1 = dist(rng); 
float random_value_2 = dist(rng); 

您可以在主機和設備代碼中使用它。 看看Thrust examples。