我使用带有反向欧拉积分的质量弹簧模型来实现基于 gpu 的布料模拟。使用带滤波器的共轭梯度求解器求解线性系统。一切都在 GPU 中完成。我认为作为一种隐式集成，它应该像许多论文指出的那样稳定，但就像隐式方法一样不稳定。大多数时候，当时间步长达到某个值（取决于刚度）时，CG 拒绝收敛并进入无限循环。
尽管我一遍又一遍地检查代码并阅读了很多论文，但我仍然找不到原因。布料正确移动，动画比使用显式积分的动画更令人信服，所以我猜力的计算是正确的。我错过了什么导致不稳定吗？
我已经被这个奇怪的问题解决了好几天......
有人可以帮忙吗？任何建议我都会非常感激
提前非常感谢！