小编Qua*_*tes的帖子

在我的c ++应用程序中,我有一个范围(0,1)的双精度矢量,我必须尽可能准确地计算它的总和.感觉这个问题应该先解决,但我找不到任何东西.

显然迭代矢量上的每个项目并且如果矢量大小很大并且有些项目明显小于其他项目,则执行sum + = vect [i]会累积一个重大错误.

我目前的解决方案是这个功能:

double sumDoubles(vector<double> arg)// pass by copy
{
  sort(arg.rbegin(),arg.rend());  // sort in reverse order
  for(int i=1;i<=arg.size();i*=2)
    for(int j=0;j<arg.size()-i;j+=(2*i))
        arg[j]+=arg[j+i];
  return arg[0];
}

基本上它按升序对输入进行排序并计算成对总和:

A + B + C + d + E + F + G + H =((A + B)+(C + d))+((E + F)+(G + H))

就像构建二叉树一样,但要做到位.排序应确保在每个步骤中两个加数的大小相当.

上面的代码确实比具有累积和的单个循环执行得更好.但是我很好奇是否可以进一步提高精度而不会降低性能太多.

我的理解是,现代c ++编译器采用以下方式的快捷方式:

if(true)
{do stuff}

但是如下:

bool foo(){return true}
...
if(foo())
{do stuff}

Or:

class Functor
{

 public:
        bool operator() () { return true;}

}

...

Functor f;

if(f()){do stuff}