OpeMP并行程序设计时需要注意的事项

OpeMP并行程序设计时需要注意的事项

OpeMP程序设计时需考虑一下两点：
1）当循环次数较少时，如果分成过多的线程来执行的话，可能会使得总的运行时间高于较少线程或一个线程的执行情况，并且会增加能耗；

2）如果设置的线程数量远大于CPU的核数的话，那么存在着大量的任务切换和调度的开销，也会降低整体的效率。

OpenMP只能并行化for循环，它不会并行while和do-while循环，而且只能并行循环次数在for循环外面就确定了的for循环。

    （2）、循环变量只能是整型和指针类型（不能是浮点型）

（3）、循环语句只能是单入口单出口的。循环内部不能改变index，而且里面不能有goto、break、return。但是可以使用continue，因为它并不会减少循环次数。另外exit语句也是可以用的，因为它的能力太大，他一来，程序就结束了。

3）single命令

Single用来指定某块程序由一个线程执行，除非使用了nowait子句，否则在此线程执行期间其他线程都处于等待状态，直到遇到single构造快处隐含的barrier。

在并行区域内使用共享变量时，如果存在写操作，必须对共享变量加以保护；否则，不要轻易使用共享变量，尽量将共享变量的访问转化为私有变量的访问或者对变量进行加锁操作。

循环迭代变量在循环构造区域里是私有的。声明在循环构造区域内的自动变量是私有。

注意，如果在并行区域内不加锁保护就直接对共享变量就进行写操作，存在数据竞争问题，会导致不可预测的异常结果。共享数据作为private,firstprivate ,lastprivate,threadprivate和reduction子句的参数进入并行区域后，就变成线程私有，不需要加锁保护。

任务调度

OpenMP中，任务调度主要使用在并行的for循环中。当循环由每次迭代的计算量不等时，如果简单的将各个线程分配相同的次数的迭代，则会造成各个线程的计算负载不均衡，使得有些线程先执行完，有些后执行完，造成某些CPU核空闲，影响程序性能。

线程同步

1）barrier 是线程同步命令，但是开销非常大，如果大量使用会引起程序的性能急速下降，因此在程序中应尽量减少barrier命令的使用。

2）critical命令，用在一段代码临界区之前，临界区在同一时间内只能有一个线程执行它，其他线程要执行临界区的线程则需要排队来执行它。

需要注意的是，critical命令后的结构快不能从块中跳出来，也不能从外面跳进去。

3）atomic命令

atomic命令也是对临界区的一种操作，与critical命令不同的是它对应于操作系统中的原子操作。

 

atomic功能最强大的在于能够实现对数组进行更新，看如下代码：

#include <stdio.h> 
#include <stdlib.h> 
#include <omp.h> 
void main() 
{
	float *x, *y, *work1, *work2;
	int *index;
	int n, i;
	n = 10;
	x = (float*)malloc(n*sizeof(float));
	y = (float*)malloc(n*sizeof(float));
	work1 = (float*)malloc(n*sizeof(float));
	work2 = (float*)malloc(n*sizeof(float));
	index = (int*)malloc(10 * sizeof(float));
	for (i = 0; i < n; i++) {
		index[i] = (n - i) - 1;
		x[i] = 0.0;
		y[i] = 0.0;
		work1[i] = i;
		work2[i] = i*i;
	}
#pragma omp parallel for  shared(x,y,index,n) 
	for (i = 0; i < n; i++) {
#pragma omp atomic 
		x[index[i]] += work1[i];
		y[i] += work2[i];
	}
	for (i = 0; i < n; i++)
		printf("%d %g %g\n", i, x[i], y[i]);
}