C OpenMP指令

2019-09-03 00:07:28 阅读：141 来源： 互联网

标签：c openmp

我有一个循环,我正在尝试并行化,并在其中我填充容器,说一个STL地图.然后考虑下面的简单伪代码,其中T1和T2是一些任意类型,而f和g是整数参数的一些函数,分别返回T1,T2类型：

#pragma omp parallel for schedule(static) private(i) shared(c)
for(i = 0; i < N; ++i) {
   c.insert(std::make_pair<T1,T2>(f(i),g(i))
}

这看起来相当直观,似乎应该是平凡的并行化,但它并没有像我预期的那样加速.相反,它导致我的代码中的运行时错误,因为容器中填充了意外值,可能是由于竞争条件.我甚至尝试过设置障碍而不是什么,但一切都没有用.允许它工作的唯一方法是使用如下的关键指令：

#pragma omp parallel for schedule(static) private(i) shared(c)
for(i = 0; i < N; ++i) {
#pragma omp critical
   {
      c.insert(std::make_pair<T1,T2>(f(i),g(i))
   }
}

但是这种渲染在上面的例子中使用omp的全部意义毫无用处,因为一次只有一个线程正在执行循环的大部分(容器插入语句).我在这里错过了什么？如果不改变编写代码的方式,有人可以解释一下吗？

解决方法:

除非f()和g()是极其昂贵的函数调用,否则这个特殊的例子并不适合并行化.

> STL容器不是线程安全的.这就是你获得比赛条件的原因.因此,访问它们需要同步 – 这使您的插入过程本身是顺序的.
>正如另一个答案所提到的那样,并行性有很多开销.因此,除非f()和g()非常昂贵,否则你的循环没有做足够的工作来抵消并行性的开销.

现在假设f()和g()是非常昂贵的调用,那么你的循环可以像这样并行化：

#pragma omp parallel for schedule(static) private(i) shared(c)
    for(i = 0; i < N; ++i) {
        std::pair<T1,T2> p = std::make_pair<T1,T2>(f(i),g(i));

#pragma omp critical
       {
          c.insert(p);
       }
    }

标签：c,openmp
来源： https://codeday.me/bug/20190902/1794081.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

C OpenMP指令