OpenMP并行for循环运行一段时间后仅单线程执行问题求助

结论 #

首先可以明确：传递const引用不会导致OpenMP并行失效，const引用仅为只读参数的传递方式，不会触发隐式线程同步或互斥，你怀疑的这个点不成立。

可能的问题根源 #

• 最常见的原因是动态调度下的任务分配不均：你当前使用schedule(dynamic, 10)的调度策略，每次给线程分配10个连续迭代块。当循环执行到后期，剩余未分配的迭代总数不足24*10时，部分线程干完所有分配的任务后就会进入闲置状态；如果不同迭代的实际计算量差异极大，还会出现快线程早早干完所有任务、只有少数慢线程还在处理高耗时迭代的情况，完全符合你观测到的“轮流只有单个线程输出结果、单核心占用高”的现象。
• 内部函数隐含锁竞争：如果你的type_cluster、cluster_breakpoints、dynamic_programming三个函数内部隐式调用了带全局锁的接口（比如标准库的rand()、无线程安全版本的内存分配器、全局静态变量读写等），当某段时间这类接口调用频率极高时，所有线程会排队等锁，表现出来就是只有单线程在执行。
• 内存分配瓶颈：如果你的Cluster对象体积较大，循环内部频繁的vector扩容、对象拷贝会触发大量内存申请，glibc默认的malloc实现存在全局锁，高频率申请时也会出现近似串行执行的现象。

排查方案 #

• 先确认position_clusters的总规模，如果总迭代数本身不大，运行后期剩余迭代不足以喂满24个线程是正常现象，不属于bug。
• 将调度策略改为schedule(static)重新测试，如果修改后后期不再出现单线程执行的问题，即可确认是动态调度下任务计算量不均导致的，可通过调整动态调度的块大小（比如把10改成1）优化负载均衡。
• 给循环内部用到的所有vector提前调用reserve()预留容量，避免频繁扩容触发内存申请，测试问题是否缓解。
• 检查三个内部调用函数的实现，确认是否用到全局静态变量、带锁的系统接口，移除这类依赖即可解决锁竞争导致的串行化问题。

内容的提问来源于stack exchange，提问作者Paghillect