3d SVD

请问一下,taichi.h 里面的svd分解为什么跟Eigen里面的svd分解不一样呢? 同一个算法,前者能顺利运行,后者算一会就炸。
Besides, taichi.h 有3阶矩阵的svd分解吗?

有的:taichi/sifakis_svd.h at master · taichi-dev/taichi · GitHub
你的 Eigen 的计算是放在 Taichi kernel 里面进行的嘛?

对。 我在C++版的2D case里面,把taichi.h自带的svd分解那两行代码替换成为了Eigen的svd分解,但是计算的时候,body触碰到边界的时候很快就爆炸了。 使用taichi.h自带的svd分解,代码运行没有这个问题。