C++如何优化2D粒子碰撞模拟性能 实现5000+粒子流畅运行

你当前程序帧率下降的核心原因是碰撞检测采用了全量两两比对的O(n²)暴力算法：500个粒子每帧需要执行25万次碰撞判断，当粒子数涨到5000时，判断次数直接飙升到2500万次，运算量翻100倍，必然出现明显卡顿。以下是可落地的优化路径，改完后流畅跑5000粒子完全没有压力。

1 最高优先级优化：用均匀空间网格把碰撞检测复杂度降到O(n) #

这是2D粒子仿真领域最通用、收益最高的优化手段，改完性能直接提升1~2个数量级：

• 网格尺寸设置为粒子直径（你当前粒子size为3，直径6，网格就设为6×6像素），保证任何发生碰撞的两个粒子，一定位于同一个网格或者相邻网格中。
• 每帧执行碰撞检测前，先清空网格缓存，遍历所有粒子，根据粒子的(x,y)坐标计算它所属的网格ID，把粒子的索引/指针存入对应网格的粒子列表中。
• 碰撞检测阶段，不再遍历全部粒子，对每个粒子，仅检查它所在网格+周围8个相邻网格内的粒子即可，距离超过1个网格的粒子完全不可能发生碰撞，直接跳过。
  按这个逻辑，每个粒子平均只需要和3~10个邻近粒子做碰撞判断，5000粒子每帧总判断量仅几万次，相比原来的2500万次开销可以忽略。

2 砍掉现有代码中的无效开销，零成本提性能 #

你当前代码里有不少浪费性能的写法，改掉之后不需要改架构就能提不少帧率：

• 修复双重循环的重复计算：当前两层循环都从0遍历到数组末尾，i=1,j=2和i=2,j=1是同一对粒子，会重复做两次碰撞判断，直接把内层循环改成从i+1开始遍历，直接砍掉一半运算量。另外全局的temp变量没必要，碰撞交换速度时用局部变量即可，全局变量访存效率更低。
• 替换错误且低效的碰撞判定：你当前写的四个方向轴对齐碰撞判断既不符合圆形粒子的碰撞逻辑，判断分支也多。直接用平方距离判定碰撞，省掉开根号的开销：

float dx = i.x - j.x;
float dy = i.y - j.y;
float dist_sq = dx*dx + dy*dy;
float r_sum = i.size + j.size;
if (dist_sq <= r_sum * r_sum) {
    // 此处写碰撞响应逻辑
}

• 移除重复的随机数初始化：不要在gas::p函数里反复调用srand(time(0))，time(0)精度为秒，同一秒内初始化的粒子随机值会完全重复，且函数调用本身有开销，随机种子在main函数开头初始化一次就够。
• 修复启动延迟bug：你当前glutTimerFunc第一次调用的延迟设为1000ms，会导致程序打开后第一秒空白，改成和后续帧一致的1000/60即可。

3 面向后续功能扩展的进阶优化 #

如果后续还要加粒子相互作用、渲染特效等功能，可以提前做这些优化：

• 替换立即模式渲染：你当前用的glBegin/glEnd逐顶点提交的老OpenGL接口，粒子数多了渲染本身就会成为瓶颈，换成顶点缓冲区（VBO）批量提交所有粒子的顶点数据，5000个点的渲染开销几乎可以忽略。
• 多线程并行计算：你代码已经引入了thread头文件，可以把网格按块拆分，分配给不同CPU核心并行计算粒子运动、碰撞检测，在多核CPU上性能还能再提3~5倍。
• 修正碰撞响应逻辑：你当前直接交换两粒子速度的写法，仅适用于质量完全相同的粒子正碰场景，斜碰、质量不同时物理效果完全错误，后续加功能时可以补全基于动量守恒的二维弹性碰撞公式，同时避免同一帧内同一对粒子被重复处理导致的速度异常。
• 稳定时间步：当前固定t=1的时间步在粒子速度较快时容易出现穿模，可以把每帧的物理更新拆成2~3个小步长迭代，不需要加复杂的连续碰撞检测就能解决大部分穿模问题。

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