使用for循环调用@ti.kernel修饰的类函数结果不一致

综合讨论区
,
#1

大家好，我在使用taichi进行mpm模拟的过程中发现以下两种情况的运行结果不同，请问是什么原因造成的？

图片
图片1910×1043 295 KB

#2

是很大的区别吗？有时候在并行计算两次结果是会不同的，带有一定随机性的。
不过一般来讲这两种运行方式应该不会有什么区别，我自己试了我的代码，是没啥区别的

#4

区别还是蛮大的，case1会出现数值发散的行为，而case2就不会。不过我发现如果我不用@ti.kernel修饰step函数，而是修饰里面的每一个子函数的时候，这两种情况的运行结果是一致的