使用/为TBB构建分组数据结构

时间:2015-06-24 00:44:42

标签: c++ parallel-processing tbb

最近我一直在考虑使用TBB而不是boost.threads来加速开发。通常,parallel_for在大多数情况下都有效,但我的情况有点复杂。

有一个结构数组需要根据成员变量进行排序的计算。这是因为变量值与计算期间访问的数据有关,根据这种情况分组结构将允许串行设计中的高速缓存一致性。

#include <tbb/tbb.h>
#include <iostream>

struct thing
{
  float value_one;
  float value_two;
  unsigned int sort_id;
};

class functor
{
  thing* m_array;
public:
  functor(thing* _array) : m_array(_array) {;}
  void operator()(const tbb::blocked_range<unsigned int>& r) const
  {
    for(int i = r.begin(); i != r.end(); ++i)
    {
      //Doing a computation with array
      m_array[i].value_one = m_array[i].value_two * m_array[i].value_two;
    }
  }
};

int main(int argc, char const *argv[])
{
  unsigned int n = 10;
  thing* array = new thing[n];

  // Note the ordered id groups
  array[0].sort_id = 1;
  array[1].sort_id = 1;
  array[2].sort_id = 1;
  array[3].sort_id = 2;
  array[4].sort_id = 3;
  array[5].sort_id = 5;
  array[6].sort_id = 5;
  array[7].sort_id = 9;
  array[8].sort_id = 9;
  array[9].sort_id = 9;

  // Do something parallel with array here...
  // parallel_for(tbb::blocked_range<unsigned int>(0, n, 2), functor(array));

  delete[] array;

  return 0;
}

上面给出了一个简化的例子,但实际上我很可能有一个大约30-60万个元素的数组。

据我所知,parallel_for会将数组划分为分组范围。但是我希望每个范围都包含特定id的所有结构。我不介意范围是否包含多个id的结构,只要它们是顺序的并且包含这两个id的所有结构。

int count = 0;
thing** blocks = new thing*[7];
int* sizes = new int[7];

int current_id = 0;
for(unsigned int i = 0; i < n; ++i)
{
  if(array[i].sort_id != current_id)
  {
    current_id = array[i].sort_id;
    blocks[count] = &array[i];
    sizes[count] = 1;
    ++count;
  }
  else
  {
    sizes[count - 1] += 1;
  }
}

parallel_for(tbb::blocked_range<unsigned int>(0, count, 2), functor(blocks, sizes));

我是否应该以某种方式将数组划分为另一个数组所指向的较小块然后并行化(如上面的代码所示),如果是这样,那么这样做的有效方法是什么,或者给出最佳的示例?是否有替代parallel_for(例如task_group)更适合这个问题?

1 个答案:

答案 0 :(得分:1)

这个问题对我来说仍然不太清楚,因为你正在混合目标和可能的方法。

如果需要对数组进行排序,可以使用parallel_sort

如果需要为排序数组构建索引,其中sort_id作为键映射到给定元素组所在的索引,并使用concurrent_unordered_map来存储组(如果有大量组)并使用parallel_for来构建它。

如果组的数量少于大约数百个,您可以使用std::mapstd::unordered_map并使用parallel_reduce来构建部分地图并将它们合并到最终状态。

最后,如果您拥有适合组的数据结构,则可以根据需要在组中使用parallel_for

P.S。这样:

  

分组范围称为将被添加到堆栈进行计算的任务

对我来说听起来很奇怪。有一个用户仿函数(或C ++ 11 lambda),可以并行调用来处理不同的范围[begin;end)。如果您将仿函数称为“任务”,则可以,但不要将其与tbb::task混合使用。