C++ gamma_distribution返回无穷大

c++ gamma_distribution returning infinity

本文关键字：返回 无穷大 distribution gamma C++      更新时间：2023-10-16

我最近写了一些代码，需要从伽马分布中大量绘制。我使用标准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 << " millionn";
        // 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 << "endn";
    std::cin.get();
    return 0;
}

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

目前，我已经通过捕获无限值并使它们成为大数字来修补这个问题。但是，如果有人对为什么/如何发生这种情况有任何见解，我们将不胜感激！

我用以下方法确定了随机种子：

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

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

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

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

相关文章：