OpenMP基础：编程模型、编译指令、并行区域、线程管理
好，咱们开始聊OpenMP。说实话，我第一次接触OpenMP是在一个图像处理项目里，当时老板催得紧，串行代码跑一张4K图要十几秒，我急得直挠头。后来一个老同事说：「你试试OpenMP，加一行指令就行。」我当时半信半疑，结果一用，真香。
OpenMP说白了就是一套给C/C++和Fortran用的并行编程接口。它不是新语言，而是编译器指令的集合。你想想看，你原本的串行代码，加几行#pragma omp，就能自动跑多线程了。嗯，就是这么简单粗暴。
2.1 OpenMP编程模型
OpenMP基于Fork-Join模型。什么意思呢？程序启动时只有一个主线程（Master Thread），遇到并行区域时，主线程会派生出（Fork）一组线程来干活。干完活，这些线程再合并（Join）回主线程。
核心要点：
- 串行部分：只有主线程在执行
- 并行区域：多个线程同时干活
- 隐式同步：并行区域结束时，所有线程自动同步
我在项目中遇到过一个问题：有人以为并行区域里的变量都是共享的，结果数据乱得一塌糊涂。其实OpenMP对变量的处理分两种——共享（shared）和私有（private）。默认情况下，定义在并行区域外的变量是共享的，区域内的局部变量是私有的。这个坑，我踩过。
2.2 编译指令：#pragma omp
OpenMP的编译指令格式很固定：
#pragma omp 指令名 [子句[, 子句]...]
常用的指令有这些：
| 指令 | 作用 |
|---|---|
parallel |
定义一个并行区域 |
for |
将for循环分配给多个线程执行 |
sections |
将不同代码块分配给不同线程 |
single |
指定某段代码只由一个线程执行 |
master |
指定某段代码只由主线程执行 |
举个例子，最简单的并行区域：
#include <stdio.h>
#include <omp.h>
int main() {
#pragma omp parallel
{
int tid = omp_get_thread_num();
printf("Hello from thread %d\n", tid);
}
return 0;
}
编译时记得加-fopenmp（GCC）或/openmp（MSVC）。我第一次编译忘了加，结果程序还是串行跑的，我还纳闷怎么没加速……
2.3 并行区域：parallel指令
#pragma omp parallel是OpenMP的基石。它后面跟的代码块，会被多个线程同时执行。每个线程都跑同样的代码，但处理的数据可以不同。
你可以用子句来控制并行区域的行为：
num_threads(N)：指定线程数if(condition)：条件并行，满足条件才开多线程private(var)：让每个线程拥有变量的私有副本shared(var)：显式声明变量为共享default(none)：要求显式声明所有变量的属性
我的习惯：写并行区域时，我总爱加default(none)。虽然多打几个字，但能逼自己把每个变量的属性想清楚。这习惯帮我避免了好几次数据竞争。
看个带子句的例子：
#pragma omp parallel num_threads(4) default(none) shared(a, b) private(i)
{
int tid = omp_get_thread_num();
i = tid;
a[i] = b[i] * 2;
}
2.4 线程管理
OpenMP的线程管理，说白了就是控制「谁干活、干多少、怎么分」。常用的API函数有：
| 函数 | 作用 |
|---|---|
omp_get_thread_num() |
获取当前线程ID（0到N-1） |
omp_get_num_threads() |
获取并行区域中的线程总数 |
omp_set_num_threads(N) |
设置后续并行区域的线程数 |
omp_get_max_threads() |
获取系统允许的最大线程数 |
omp_in_parallel() |
判断当前是否在并行区域内 |
线程数怎么设？我一般遵循这个原则：
- CPU密集型任务：线程数 = 物理核心数
- 混合型任务（有I/O）：线程数 = 物理核心数 × 2
- 千万别设超过硬件线程数太多，否则上下文切换会吃掉性能
我曾经踩过的坑：有次我设了num_threads(128)，以为线程越多越快。结果我的CPU只有8核，128个线程互相抢资源，性能反而比串行还差。嗯，从那以后我学乖了，先查omp_get_max_threads()再设线程数。
还有一个细节：线程ID从0开始，主线程的ID是0。你可以在并行区域里用if(tid == 0)来让主线程干一些特殊的事，比如初始化、汇总结果。
2.5 环境变量控制
除了在代码里设线程数，你还可以通过环境变量来控制OpenMP的行为：
OMP_NUM_THREADS：设置默认线程数OMP_SCHEDULE：设置循环调度策略OMP_DYNAMIC：是否允许动态调整线程数OMP_NESTED：是否允许嵌套并行
比如在Linux终端里：
export OMP_NUM_THREADS=8
./my_program
这样就不用改代码了。我个人习惯在开发阶段用代码硬编码线程数做测试，上线前再改成环境变量控制，灵活一些。
2.6 小结
这一章我们聊了OpenMP的Fork-Join模型、编译指令怎么写、并行区域怎么用、线程怎么管。你想想看，其实核心就三件事：
- 用
#pragma omp parallel划出并行区域 - 用子句控制变量属性和线程数
- 用API函数获取线程信息
下一章我们会深入#pragma omp for，看看怎么把循环并行化。到时候我会分享一个我优化图像滤波的实战案例，保证干货满满。