OpenMP并行编程简易教程.docVIP

OpenMP并行编程简易教程(一)近来在学校做实习，需要使用OpenMP来优化一些SAT问题的算法，便接触到了OpenMP，它是基于多处理器的,分享同一内存的,并行的编程模式。对于多核的CPU计算机，仅需要通过增加简单的语句，便可以得到立竿见影的性能提升。本文教程由自己的学习笔记而来，通过4个例子简单的阐述OpenMP的基本编程语句。例子1 helloworld.c01 #include 02 int main( argc, argv)03 int argc;04 char **argv;05 {06 #pragma omp parallel07 printf( Hello world!\n );08 return 0;09 }程序helloworld.c的本质就是在屏幕上打印字符串”Hello world!”。”#pragma omp parallel”是一条OpenMP标准的语句，它的含义是让它后面的语句按照多线程来执行。需要注意的是每个线程都去做相同的事情。[root@localhost zf]# gcc –fopenmp helloworld.c[root@localhost zf]# ./a.outHello world!编译，执行程序。屏幕上打印出了一遍”Hello world”。-fopenmp 是gcc编译支持OpenMP程序的参数，gcc4.2以上的版本默认支持OpenMP。由于系统环境变量中的NUM_OMP_THREADS的值默认为1，所以程序在执行的时候只使用了一个线程来执行。[root@localhost zf]# NUM_OMP_THREADS=5[root@localhost zf]# export NUM_OMP_THREADS[root@localhost zf]# ./a.outHello world!Hello world!Hello world!Hello world!Hello world!给表示线程数量的环境变量NUM_OMP_THREADS赋值为5并导出，再执行程序，得到5遍的”Hello world!”，说明程序将打印语句用5个线程分别执行了一遍。当然，我们不希望受到运行系统环境变量的限制，所以我们也可以通过替换06行代码为” #pragma omp parallel num_threads(10)”，编译之后再执行程序，我们得到10遍的”Hello world!”，这时无论环境变量NUM_OMP_THREADS的值为多少，我们依然只能得到10遍的”Hello world!”。例子2 RankNum.c01 #include 02 #include 03 int main( int argc, char **argv)04 {05 int rank, size;06 #pragma omp parallel private(rank)07 {08 rank= omp_get_thread_num();09 size= omp_get_num_threads();10 printf( Hello world! Im %d of %d\n, rank, size );11 }12 return 0;13 }导入omp.h文件以支持OpenMP的函数。函数omp_get_threads_num() 返回当前的线程标号。函数omp_get_num_threads() 返回当前的线程数量。（不再程序中定义线程数，使用环境变量中的数据，当前是10）[root@localhost zf]# gcc –fopenmp RankNum.c[root@localhost zf]# ./a.outHello world! I’m 8 of 10Hello world! I’m 0 of 10Hello world! I’m 4 of 10Hello world! I’m 6 of 10Hello world! I’m 3 of 10Hello world! I’m 2 of 10Hello world! I’m 7 of 10Hello world! I’m 1 of 10Hello world! I’m 9 of 10Hello world! I’m 5 of 10编译，执行，得到10条输出，10个线程分别执行了打印语句。线程执行的顺序随机，与编号无关。未完待续…… OpenMP并行程序设计（一） OpenMP是一个支持共享存储并行设计的库，特别适宜多核CPU上的并行程序设计。今天在双核CPU机器上试了一下OpenMP并行程序设计，发现效率方面超出想象，因此写出来分享