1

テスト コードで std::accumulate を使用すると、予期しない結果が得られます。double の大きなベクトルを合計しようとしていますが、何らかの理由で値がオーバーフローしています: 

#include <iostream>
#include <vector>
#include <functional>
#include <numeric>

using namespace std;

double sum(double x, double y)
{
    // slows things down but shows the problem:
    //cout << x << " + " << y << endl;
    return (x+y);
}

double mean(const vector<double> & vec)
{
    double result = 0.0;

    // works:
    //vector<double>::const_iterator it;
    //for (it = vec.begin(); it != vec.end(); ++it){
    //      result += (*it);
    //}

    // broken:
    result = accumulate(vec.begin(), vec.end(), 0, sum);

    result /= vec.size();

    return result;
}


int main(int argc, char ** argv)
{

    const unsigned int num_pts = 100000;

    vector<double> vec(num_pts, 0.0);

    for (unsigned int i = 0; i < num_pts; ++i){
        vec[i] = (double)i;
    }

    cout << "mean = " << mean(vec) << endl;

    return 0;
}

合計内の cout からの部分的な出力:

2.14739e+09 + 65535
2.14745e+09 + 65536
-2.14748e+09 + 65537
-2.14742e+09 + 65538
-2.14735e+09 + 65539

正しい出力 (反復):

平均 = 49999.5

誤った出力 (蓄積を使用):

平均 = 7049.5

私はおそらく疲れた間違いを犯していますか？以前に蓄積を正常に使用しました...

ありがとう

4

1 に答える 1

9

doubleに aを渡す必要がありますaccumulate:

result = accumulate(vec.begin(), vec.end(), 0.0, sum);
                                            ^^^

それ以外の場合は、 を使用して累積が実行されint、結果が double に変換されます。

于 2013-06-05T21:35:27.140 に答える