(一) OpenMP介绍

OpenMP是一种用于共享内存并行系统的多线程程序设计方案，支持的编程语言包括C、C++。OpenMP提供对并行算法的高层抽象描述，特别适合在多核CPU机器上的并行程序设计。进程是CPU资源分配的最小单位，线程是CPU调度的最小单位。线程是建立在进程的基础上的一次程序运行单位，进程是应用程序的执行实例。

OpenMP采用fork-join的执行模式，开始的时候只存在一个主线程，当需要进行并行计算的时候，派生出若干个分支线程来执行并行任务。当并行代码执行完成之后，分支线程会合，并把控制流程交给单独的主线程。

编译制导指令以#pragma omp开始，后边跟具体的功能指令，格式如:#pragma omp 指令[子句 子句]

(1) 用户定义声明，使用#pragma omp declare reduction。

#pragma omp declare reduction(omp_insert:std::vector:omp_out.insert(omp_out.end(),omp_in.begin(),omp_in.end()))
std::vector vec;
#pragma omp parallel for reduction(omp_insert:vec)
for(int i=0;i<1000;i++) vec.push_back(i);

(2) omp中使用parallel制导指令标识代码中的并行段，形式为：

#pragma omp parallel
{
   每个线程都会执行大括号里的代码
}

omp_set_num_threads(4); //设置线程数，一般设置的线程数不超过CPU核心数，这里开4个线程执行并行代码段

#pragma omp parallel
for(int i=0;i<2;i++)
   printf("i = %d, I am Thread %d\n", i, omp_get_thread_num());
i = 0, I am Thread 0 
i = 0, I am Thread 1 
i = 1, I am Thread 0 
i = 1, I am Thread 1 
i = 0, I am Thread 2 
i = 1, I am Thread 2 
i = 0, I am Thread 3 
i = 1, I am Thread 3

#pragma omp parallel for
for (int i = 0; i < 6; i++)
printf("i = %d, I am Thread %d\n", i, omp_get_thread_num());
i = 4, I am Thread 2
i = 2, I am Thread 1
i = 0, I am Thread 0
i = 1, I am Thread 0
i = 3, I am Thread 1
i = 5, I am Thread 3
线程0执行i=0和1,线程1执行i=2和3，线程2执行i=4，线程3执行i=5,线程0就是主线程