gcc - boost::uniform_on_sphere 在几百万次正确实现后突然失败，但仅限于某些主机

Question

问题

在正确生成二维随机向量一段时间后，boost::uniform_on_sphere 分布突然生成一个值为 -nan 的向量。我已经在三台机器上测试了包含的程序 - 在其中两台机器上观察到错误，但在第三台机器上没有观察到错误。有人知道会发生什么吗？

编辑:如果使用相同的类型，它会发生在所有主机上。

主办方

主机1

• AMD Opteron(tm) 处理器 6174
• g++ (GCC) 4.4.6 20120305(红帽 4.4.6-4)
• 实现 3802480 次后失败

主机2

• 英特尔(R) 酷睿(TM) i5 CPU 650 @ 3.20GHz
• g++ (GCC) 4.7.2 20120921(红帽 4.7.2-2)
• 实现 3802480 次后失败

主机3

• 英特尔(R) Atom(TM) CPU D2700 @ 2.13GHz
• g++(Ubuntu/Linaro 4.6.3-1ubuntu5)4.6.3
• 在 3802480 次实现后失败，输出略有不同

实现

在主机 1 和 2 上:

3802470: -4.8961880803e-01 , -8.7193667889e-01
3802471: 9.9225074053e-01 , -1.2425158173e-01
3802472: 6.5411877632e-01 , -7.5639182329e-01
3802473: -9.8332953453e-01 , -1.8183253706e-01
3802474: 7.1217632294e-01 , -7.0200067759e-01
3802475: -9.9968332052e-01 , 2.5166392326e-02
3802476: 9.9412262440e-01 , 1.0826008022e-01
3802477: -6.2786966562e-01 , 7.7831840515e-01
3802478: 5.7143938541e-01 , 8.2064425945e-01
3802479: 5.8261138201e-01 , 8.1275087595e-01
3802480: -nan , -nan
3802481: 9.2151606083e-01 , 3.8834002614e-01
3802482: 8.6448800564e-01 , -5.0265353918e-01
3802483: -9.1891586781e-01 , 3.9445358515e-01
3802484: -9.1544634104e-01 , 4.0244001150e-01

在主机 3 上，使用 float 而不是 float_t:

3802470: -4.8961877823e-01 , -8.7193661928e-01
3802471: 9.9225074053e-01 , -1.2425158918e-01
3802472: 6.5411871672e-01 , -7.5639182329e-01
3802473: -9.8332953453e-01 , -1.8183253706e-01  <- exactly the same as above
3802474: 7.1217626333e-01 , -7.0200061798e-01
3802475: -9.9968332052e-01 , 2.5166388601e-02
3802476: 9.9412262440e-01 , 1.0826008022e-01
3802477: -6.2786966562e-01 , 7.7831846476e-01
3802478: 5.7143932581e-01 , 8.2064431906e-01   <- slightly different
3802479: 5.8261138201e-01 , 8.1275087595e-01   <- exactly the same
3802480: -nan , -nan
3802481: 9.2151612043e-01 , 3.8834002614e-01
3802482: 8.6448800564e-01 , -5.0265347958e-01
3802483: -9.1891586781e-01 , 3.9445355535e-01
3802484: -9.1544634104e-01 , 4.0244001150e-01

计划

这是简单地使用 g++ bug.cpp 编译的。打开 -O3 优化并没有改变结果。

#include <boost/circular_buffer.hpp>
#include <boost/random/variate_generator.hpp>
#include <boost/random/uniform_on_sphere.hpp>
#include <boost/random.hpp>
#include <iostream>
#include <fstream>
using namespace std;

int main(int argc, const char *argv[])
{
    typedef boost::mt19937                                              GeneratorType;
    typedef boost::uniform_on_sphere<float_t>                           DistributionType;
    typedef boost::variate_generator<GeneratorType, DistributionType >  VariateType;
    typedef boost::circular_buffer<DistributionType::result_type>       BufferType;
    GeneratorType       gen;
    DistributionType    dist(2);
    VariateType         variate(gen,dist);
    const int           BUFSIZE = 10;

    gen.seed(11);
    BufferType buf(BUFSIZE);
    long n(0);
    while (1){
        cout << "n: " << n << "\r" << flush;
        DistributionType::result_type tmp = variate();
        buf.push_back(tmp);
        if (isnan(tmp[0])) {
            cout << "n: " << n << "       " << endl;
            cout << tmp[0] << " , " << tmp[1] << endl;
            ofstream fout("debug.out");
            for (int i=0; i<BUFSIZE; i++)
                fout << buf[i][0] << " " << buf[i][1] << endl;
            fout.close();
            ofstream gen_out("gen.out");
            gen_out << gen;
            gen_out.close();
            exit(1);
        }
        n++;
    }

    return 0;
}

我非常感谢任何帮助!

Answer 1

这似乎是由于 boost/random/uniform_on_sphere.hpp 中的错误而发生的。在 N 维中，生成 N 个正态分布数作为 N 向量的分量，然后对其进行归一化。在 2d 中，获得两个零的概率显然不可忽略，由于归一化，导致每个分量的计算为 0/0=NaN。

解决方法是针对小尺寸手动编程此分布。