关于OpenMP的归约操作reduction

最新推荐文章于 2024-04-30 11:04:25 发布

转载最新推荐文章于 2024-04-30 11:04:25 发布 · 1.1k 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/cckk/p/10012934.html

本文探讨了在使用OpenMP进行并行编程时常见的陷阱，特别是关于归约(reduction)指令的正确使用。通过对比不同函数调用方式，说明了直接操作变量与通过指针或引用传递参数对归约效果的影响。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这里提一个重要的点

像这样

double sum = 0;
void ff()
{
    sum += 0.5;
}
//main()
#pragma omp parallel for reduction(+:sum)
    for (int i = 1; i <= 40000; i++)
    {
        ff();
    }
    cout << sum << endl;

对归约对象的操作“隐藏”在其他函数内部，reduction(+:sum)是没有效果的！

但：

double sum = 0;
void ff(double *x)
{
    *x+= 0.5;
}
//main    
#pragma omp parallel for reduction(+:sum)
    for (int i = 1; i <= 40000; i++)
    {
        ff(&sum);
    }
    cout << sum << endl;

或：

double sum = 0;
void ff(double &x)
{
    x+= 0.5;
}
//main    
#pragma omp parallel for reduction(+:sum)
    for (int i = 1; i <= 40000; i++)
    {
        ff(sum);
    }
    cout << sum << endl;