在C++中創建統一的隨機數的最快方法是什麼？

-1

我需要在for循環中生成統一的隨機數。 for循環用於1000000個數字。還有另一個for循環，它在2000裏面。所以我生成了2 * 10^9的統一隨機數。我用下面的方法：在C++中創建統一的隨機數的最快方法是什麼？

#include <random> 
    double zeta; 
    unsigned seed = std::chrono::system_clock::now().time_since_epoch().count(); 
     auto uniform_rand = bind(uniform_real_distribution<double>(0,1), mt19937(seed)); 
    for(int j=0; j<1000000; j++) 
       for(int i=0; i<2000; i++) 
     zeta=-eta/2.0+uniform_rand()*eta; // in the range (-eta/2, +eta/2) 
theta[i]+=zeta; 
    end 
end

來源

2017-10-17 Oliver Range

'std :: bind'往往被lambdas取代，無論是在可讀性和性能方面。 – chris

這是一個非常多的隨機數字。如果你需要這麼多，爲什麼不把它分解成線程來利用多核？ – tadman

對不起，我不明白你是什麼意思@tadman –

它幾乎與你的一樣。我只是沒有看到任何需要一個活頁夾或lambda。

編輯：我也將發生器從mt19937更改爲minstd_rand，這使得代碼快了88倍。使用mt19937代碼的性能與問題中的代碼大致相同。

#include <random> 

int main() { 
    /*std::mt19937 gen(std::random_device{}());*/ 
    // Changing mt19937 to minstd_rand makes the code run 88 times faster! 
    std::minstd_rand gen(std::random_device{}()); 
    std::uniform_real_distribution<double> dist(0, 1); 
    for(unsigned int i = 0; i < 1000000; ++i) { 
     for(unsigned int j = 0; j < 2000; ++j) { 
      double anotherRandomNumber = dist(gen); 
      // Do whatever you want with generated random number. 
     } 
    } 
}

來源

2017-10-17 17:29:29 navyblue

它比我寫的真的快嗎？ –

它是否在不同的運行中創建不同的隨機數字？ –

@OliverRange是的，它每次使用'std :: random_device'來播放隨機生成器。但有時它不會工作，請參閱https://stackoverflow.com/questions/18880654/why-do-i-get-the-same-sequence-for-every-run-with-stdrandom-device-with-mingw –

使用這種算法從here：

uint64_t s[2] = { 0x41, 0x29837592 }; 

static inline uint64_t rotl(const uint64_t x, int k) { 
    return (x << k) | (x >> (64 - k)); 
} 

uint64_t next(void) { 
    const uint64_t s0 = s[0]; 
    uint64_t s1 = s[1]; 
    const uint64_t result = s0 + s1; 

    s1 ^= s0; 
    s[0] = rotl(s0, 55)^s1^(s1 << 14); // a, b 
    s[1] = rotl(s1, 36); // c 

    return result; 
} 

double uniform() { 
    return next()*(1.0/18446744073709551616.0); 
}

這比我的機器上的例子更快4倍。

注意：您需要種子s，或許與std::random_device。

來源

2017-10-17 17:55:41 geza

能否請您提供一個用於循環的答案，就像我的情況一樣？ –

@OliverRange：將這段代碼粘貼到你的代碼中，用'uniform（）'替換'uniform_rand（）'。 – geza

我不知道在哪裏發佈它。主要和之前的循環後？ –

在C++中創建統一的隨機數的最快方法是什麼？

回答

相關問題