如何使用减少并行化这个 for 循环？

Question

我试图通过使用 Openmp 使这个 for 循环并行化，我认识到这个循环中有减少，所以我添加了“#pragma omp parallel for reduction(+,ftab)”，但它没有工作，它给了我这个错误： 错误：未找到“ftab”的用户定义缩减。

   #pragma omp parallel for reduction(+:ftab)
    for (i = 1; i <= 65536; i++) ftab[i] += ftab[i-1];

Answer 1

你要做的操作是prefix sum。它可以并行完成。一种简单的方法是将thrust::inclusive_scan 与 OpenMP 或 TBB 后端一起使用。

thrust::inclusive_scan(thrust::omp::par, ftab, ftab + 65536, fab);

或

thrust::inclusive_scan(thrust::tbb::par, ftab, ftab + 65536, fab);

您也可以按照the Wikipedia page 中的说明自行实现。