3

问题

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

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

东道主

主机 1

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

主机 2

  • Intel(R) Core(TM) i5 CPU 650 @ 3.20GHz
  • g++ (GCC) 4.7.2 20120921 (红帽 4.7.2-2)
  • 3802480 实现后失败

主机 3

  • Intel(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;
}

我真的很感激任何帮助！

4

2 回答 2

1

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

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

于 2012-12-11T12:07:52.883 回答
0

x86-32 和 x86-64 之间的一大区别是 x87（80 位 FP）和 SSE（64 位 FP）之间的区别。即 32 位 Atom 有 16 个额外的位可以工作。

于 2012-12-07T14:21:34.320 回答