TBB之parallel_reduce

循环做的是reduce(最终输出的是单一的数)，例如求和：

float SerialSumFoo( float a[], size_t n ) {    float sum = 0;    for( size_t i=0; i!=n; ++i )        sum += Foo(a[i]);    return sum;}

你能使用模板类parallel_reduce并行这个循环：

float ParallelSumFoo( const float a[], size_t n ) {    SumFoo sf(a);    parallel_reduce( blocked_range<size_t>(0,n), sf );    return sf.my_sum;}

这个类SumFoo说明了reduce的细节，例如，如何对子集求和，如何把他们结合起来，这是类SumFoo的定义：

class SumFoo {    float* my_a;public:    float my_sum;    void operator()( const blocked_range<size_t>& r ) {        float *a = my_a;        float sum = my_sum;        size_t end = r.end();        for( size_t i=r.begin(); i!=end; ++i )            sum += Foo(a[i]);        my_sum = sum;    }    SumFoo( SumFoo& x, split ) : my_a(x.my_a), my_sum(0) {}    void join( const SumFoo& y ) {my_sum+=y.my_sum;}    SumFoo(float a[] ) :        my_a(a), my_sum(0)        {}};

注意不同于ApplyFoo，首先，operator()不是一个const，这是因为它一定改变SumFoo::my_sum。第二，SumFoo有一个splitting构造函数和一个join方法一定被表示为parallel_reduce工作，splitting构造函数接受指向原始对象的 引用和一个标识符split，这个标识符由TBB库定义，使用split区分splitting构造函数和拷贝构造函数。

当一个工作线程可用时，任务调度器决定让parallel_reduce调度splitting构造函数为工作线程创建一个子任务，当子任务完成时，parallel_reduce使用方法join累计子任务的结果，下图显示在一个工作线程可用的情况下split-join次序：

如果一个工作线程不可用，the second half of the iteration使用相同的body对象做的reduce，在the first half结束后the reduction of the second half才会开始。

因为相同的body可能被使用累加多个子范围，operator()不能过早抛弃累加值是重要的，下面的的代码是错误的定义SumFoo::operator()。

class SumFoo {...public:    float my_sum;    void operator()( const blocked_range<size_t>& r ) {    ...    float sum = 0; // WRONG – should be "sum = my_sum".    ...    for( ... )        sum += Foo(a[i]);    my_sum = sum;    }    ...};

上面的代码是错误的，它只返回最后子域的部分和，而不是全部子域。