布料仿真系统并行化的算法？

Question

可以使用以下算法进行简单的布料模拟：

def tick(dt):
    for p1 in particles:
        for bound in p.bounds:
            p2 = bound.particle
            p2.vel += p1.pos + bound.stable_pos - p2.pos
        p1.pos += p1.vel * dt

作为测试，我在JavaScript中尝试了implementing it。不幸的是，这种情况很糟糕。随着颗粒数量的增加，性能下降得非常快，并且极限非常低。有没有办法并行化这个算法？你能把它形容为一个简单的伪代码吗？