关于linux下的openmp编程基础[转]

OpenMP是专门针对共享地址空间的平行计算机提供的并行计算库，在Intel C++和Visual C++ 8.0里通过#pragma支持。用OpenMP，可以不必去写诸如CreateThread之类的线程管理代码，多线程程序写起来比较简洁。而且OpenMP提供了很丰富的指令，对于同步共享变量、合理分配负载等任务，都提供了有效的支持，很值得一学。不过因为这个东西用起来很轻松，稍一疏忽就能酿成愚蠢错误，所以肯定是要经常用才能驾驭好的。一开始用它，犯错误是难免的。下载安装gcc4.2后，编译的时候加上-fopenmp参数即可。
   gcc   -fopenmp   test_omp.c
   g++   -fopenmp   test_omp.cc
   gfortran   -fopenmp   test_omp.f90(95)
   尽管很多文献都说gcc4.2以上版本才支持openmp，笔者所用的gcc编译器版本为4.1.1，实验表明也支持openmp,安装fedora 7以上版本，gcc默认就支持openmp.
[xuyc@sv168 openmp]$ gcc -v
使用内建 specs。
目标： i386-redhat-linux
配置为：../configure --prefix=/usr --mandir=/usr/share/man --infodir=/usr/share/info --enable-shared --enable-threads=posix --enable-checking=release --with-system-zlib --enable-__cxa_atexit --disable-libunwind-exceptions --enable-libgcj-multifile --enable-languages=c,c++,objc,obj-c++,java,fortran,ada --enable-java-awt=gtk --disable-dssi --enable-plugin --with-java-home=/usr/lib/jvm/java-1.4.2-gcj-1.4.2.0/jre --with-cpu=generic --host=i386-redhat-linux
线程模型：posix
gcc 版本4.1.1 20070105 (Red Hat 4.1.1-52)

例子1   helloworld.c
01 #include
//#include "omp.h" 如果用intel的icc 编译器，就要写上这行
02 int main( argc, argv)
03 int argc;
04 char **argv;
05 {
06 #pragma omp parallel
07 printf( "Hello world!\n" );
08 return 0;
09 }
程序helloworld.c的本质就是在屏幕上打印字符串”Hello world!”。
“#pragma omp parallel”是一条OpenMP标准的语句，它的含义是让它后面的语句按照多线程来执行。需要注意的是每个线程都去做相同的事情。
[xuyc@sv168 openmp]$ gcc –fopenmp –o   helloworld   helloworld.c
[xuyc@sv168 openmp]$ OMP_NUM_THREADS=10
[xuyc@sv168 openmp]$ export OMP_NUM_THREADS
[xuyc@sv168 openmp]$ ./helloworld  
Hello world!
Hello world!
Hello world!
Hello world!
OMP_NUM_THREADS=2 ./helloc
如果不定义OMP_NUM_THREADS，默认会等于CPU数量，在8核心的机器上，会打印出8行"Hello World".

[xuyc@sv168 openmp]$cat /proc/cpuinfo   查看CPU核的数量，可以看到有四个核
[xuyc@sv168 openmp]$ cat /proc/cpuinfo
processor    : 0
vendor_id    : GenuineIntel
cpu family    : 6
model           : 15
model name    : Intel(R) Core(TM)2 Quad CPU Q6600   @ 2.40GHz
stepping        : 11
cpu MHz       : 2400.138
cache size    : 4096 KB
fdiv_bug        : no
hlt_bug       : no
f00f_bug        : no
coma_bug        : no
fpu          : yes
fpu_exception : yes
cpuid level     : 10
wp              : yes
flags           : fpu tsc msr pae mce cx8 apic mtrr mca cmov pat pse36 clflush d                                         ts acpi mmx fxsr sse sse2 ss ht tm pbe nx lm constant_tsc pni monitor ds_cpl vmx                                        est tm2 cx16 xtpr lahf_lm
bogomips        : 6002.73

processor    : 1
vendor_id    : GenuineIntel
cpu family    : 6
model           : 15
model name    : Intel(R) Core(TM)2 Quad CPU Q6600   @ 2.40GHz
stepping        : 11
cpu MHz       : 2400.138
cache size    : 4096 KB
fdiv_bug        : no
hlt_bug       : no
f00f_bug        : no
coma_bug        : no
fpu          : yes
fpu_exception : yes
cpuid level     : 10
wp              : yes
flags           : fpu tsc msr pae mce cx8 apic mtrr mca cmov pat pse36 clflush d                                         ts acpi mmx fxsr sse sse2 ss ht tm pbe nx lm constant_tsc up pni monitor ds_cpl                                        vmx est tm2 cx16 xtpr lahf_lm
bogomips        : 6002.73

processor    : 3
vendor_id    : GenuineIntel
cpu family    : 6
model           : 15
model name    : Intel(R) Core(TM)2 Quad CPU Q6600   @ 2.40GHz
stepping        : 11
cpu MHz       : 2400.138
cache size    : 4096 KB
fdiv_bug        : no
hlt_bug       : no
f00f_bug        : no
coma_bug        : no
fpu          : yes
fpu_exception : yes
cpuid level     : 10
wp              : yes
flags           : fpu tsc msr pae mce cx8 apic mtrr mca cmov pat pse36 clflush d                                         ts acpi mmx fxsr sse sse2 ss ht tm pbe nx lm constant_tsc up pni monitor ds_cpl                                        vmx est tm2 cx16 xtpr lahf_lm
bogomips        : 6002.73

   示例xx.c
   1   #include
   2   #include
   3   int main( int argc, char **argv)
   4   {
   5   int rank, size;
   6   #pragma omp parallel private(rank)
   7   {
   8   rank= omp_get_thread_num();
   9   size= omp_get_num_threads();
10   printf( "Hello world! I'm %d of %d\n", rank, size );
11   }
12   return 0;
13   }
[xuyc@sv168 openmp]$ gcc -fopenmp -o xx xx.c
[xuyc@sv168 openmp]$ ./xx
Hello world! I'm 0 of 4
Hello world! I'm 1 of 4
Hello world! I'm 2 of 4
Hello world! I'm 3 of 4

Microsoft Visual Studio .net 2005完全支持OpenMP编程[5]。Visual Studio. Net 2005 Professional安装之后，即可编写OpenMP程序，无须另外安装其它软件。当前的Visual Studio. Net 2005完全支持OpenMP 2.0标准。通过新的编译器选项/openmp来支持OpenMP程序的编译和连接，编译器会自动地将用户的代码和OpenMP在Windows下实现的库vcomp.dll连接在一起。程序在运行的时候会自动地寻找vcomp.dll。下面用Visual Studio.Net 2005生成一个新OpenMP项目OpenMP1。
类别：多核 查看评论