在做隐式有限元时,有个很重要的物理量是Jaccobian矩阵.例如应力张量对应变张量的偏导.若我有一个应变场 strain=ti.Matrix.field(33floatshape=(XYZ)) 和应力场 stress=field(33floatshape=(XYZ)),有一个kernel能够使用每个位置点的strain的3x3矩阵求出该点的stress 3x3矩阵.现在想求\partial stress/\partial strain,每个位置点上一个3x3x3x3的矩阵.想请教两个问题,Taichi目前没办法定义高阶张量,如何定义该Jaccobian矩阵场,使用3x3x3x3xXxYxZ的标量场来计算?若是想使用自动求导方案,不知道有没有Taichi的范式可以参考一下?