3

問題

しばらくの間、2 次元でランダムなベクトルを正しく生成した後、boost::uniform_on_sphere分布は突然値を持つベクトルを生成します-nan。含まれているプログラムを 3 台のマシンでテストしましたが、そのうちの 2 台でエラーが発生しましたが、3 台目では発生しませんでした。誰が何が起こっているのか考えていますか?

編集：同じタイプが使用されている場合、すべてのホストで発生します。

ホスト

ホスト 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 では、2 つのゼロを取得する確率は明らかに無視できないNaNため、正規化により各コンポーネントの計算は 0/0= になります。

回避策は、小さな次元に対してこの分布を手動でプログラムすることです。

于 2012-12-11T12:07:52.883 に答える
0

x86-32 と x86-64 の大きな違いの 1 つは、x87 (80 ビット FP) と SSE (64 ビット FP) の違いです。つまり、32 ビットの Atom には 16 ビットの余分なビットがあります。

于 2012-12-07T14:21:34.320 に答える