Question

我最近编写了一些代码，需要从gamma分布中进行大量绘制。我使用标准的gamma_distribution方法实现了这一点，但我偶尔发现（一旦在蓝色的月亮中），它返回的值为1.＃INF，没有明显的原因。

以下是我发现此问题的最小示例。在我的测试中，我发现问题通常发生在第十亿次迭代左右。

#include "stdafx.h"
#include <iostream>
#include <random>

std::random_device rd;
std::default_random_engine generator(rd());

int _tmain(int argc, _TCHAR* argv[])
{
    // create gamma distribution and random variable x
    std::gamma_distribution<double> rgamma(2.0,1.0);
    double x;

    // loop through a large number of iterations
    for (unsigned long int i=0; i<int(4e9); i++) {

        // print update to console every million iterations
        if ((i+1)%int(1e6)==0)
            std::cout << "iteration: " << (i+1)/1e6 << " million\n";

        // draw new value of x from gamma distribution
        x = rgamma(generator);

        // if x==infinity then break
        if ((1.0/x)==0) {
            std::cout << "Error at iteration " << (i+1) << ": x=" << x << "\n";
            std::cin.get();
            exit(1);
        }
    }
    // print message if reach end of loop
    std::cout << "end\n";
    std::cin.get();
    return 0;
}

我很想知道其他人是否可以复制这个问题。我不知道这是否相关，但上面的程序是在Visual Studio 2010中编写为Win32应用程序，并在具有8核Intel处理器的Windows机器上运行。

暂时我只是通过捕获无限值并使它们成为大数来修补此问题。但是，如果有人对此为何/如何发生有任何了解，我们将不胜感激！

Answer 1

我用随机种子确定了随机种子：

//std::random_device rd;
std::default_random_engine generator(-246744094);

并且可以在8.32亿次迭代中始终如一地重现此错误。我得到了不同的编译器（VS2013 32 / 64bit和英特尔C ++ 2016 64位），以及不同的机器（我的桌面i7和集群节点Xeon E5），结果相同。我的rgamma.max（）也同意Bob的观点。

然后我用MinGW编译它（g ++ main.cpp -o main.exe -std = c ++ 0x -O3） - 它运行时没有产生任何无限。

SO - 我假设：它是通过Visual Studio链接的库中的某种错误，无论您使用的是MS还是英特尔编译器。不确定是否可以报告/修复 - 或者是否可以更快地找到一些可靠的gamma源代码。

c ++ gamma_distribution返回无穷大

1 个答案: