C++实现太空物体引力碰撞时间计算时小数值被解析为0致程序卡死

• 最致命的逻辑错误来自gamma的初始化语句：你写的double gamma = 166857 / 2500000000000000，等号右侧是两个纯整数类型的字面量（编译器会根据数值大小自动匹配对应的整数类型存储，不会识别为浮点数），C会先执行整数截断除法，再把结果转成double。分子远小于分母的情况下，整数除法的结果直接是0，也就是说gamma从程序启动开始就完全等于0，根本不是你预期的万有引力常数换算值。你之前排查觉得是C把极小值解析为0，实际上这个问题和极小浮点数解析完全无关——gamma从根上就算成了0，后续所有加速度计算的结果都是0，速度、位移永远没有变化，r的数值永远停在初始值1.5，循环自然永远跳不出去。这也是同一段逻辑在Python能正常运行的原因——Python3的整数除法用/时是真除法，不会做整数截断，算出来的gamma是正常的小数值。
• 全局变量的类型选择存在精度隐患：你把a1、a2、v1、v2、dt、elapsedTime全部定义成了单精度float，这类类型的有效数字只有6~7位，长时间模拟过程中累计的浮点误差会远大于双精度类型，很容易出现r的变化量低于float的精度下限、看起来数值完全不动的假卡死现象。
• 物理公式存在逻辑错误：匀加速运动的位移公式为s = v0*dt + 0.5*a*dt²，你写的distance函数漏掉了加速度项前面的0.5系数，就算程序能正常跑起来，算出来的碰撞时长也会和真实物理值有偏差。

• 修正gamma的初始化逻辑，强制触发浮点除法：只要给参与除法的字面量加上浮点数后缀就行，改成double gamma = 166857.0 / 2500000000000000.0，这样运算全程是双精度浮点数，不会出现整数截断为0的问题。
• 把所有不必要的float变量替换成double类型，把代码里float a1 = 0, a2 = 0, v1 = 0, v2 = 0, dt = 0.1, elapsedTime;整行替换为double a1 = 0, a2 = 0, v1 = 0, v2 = 0, dt = 0.1, elapsedTime = 0;，避免单精度的累计精度误差，同时补上elapsedTime的显式初始值。
• 修正distance函数的位移计算逻辑，把原函数里的return (a * dt * dt + v * dt);改成return (0.5 * a * dt * dt + v * dt);，保证运动学计算的物理正确性。

补充说明：以上修复完成后程序就不会再卡死，但因为你用的是固定时间步长的欧拉法模拟，当两物体距离非常近时，单步位移可能大于两物体的剩余间距，会出现r直接从阈值上方跳到阈值下方的情况。如果需要更高精度的碰撞时间结果，可以适当缩小dt的取值，或者在循环里判断当前步的总位移如果大于剩余的r值，就按比例折算最后一段的运行时间。

内容的提问来源于stack exchange，提问作者Anders Haarberg Eriksen