python - 在 sympy 或替代方案中导出点云方程

Question

假设您R在 3D 正交坐标系中有一个点云（存储在 , 3xN 维矩阵中），并且您想导出 R 的复杂函数（包括微分、求和等）的解析表达式。

举个简单的例子：

A_il(R)=\sum_jk df(R)/dR_ij  df(R)/dR_kl 

有没有办法在 sympy （或其他任何东西......）中做到这一点，例如泛化x,y,z到X={x_i}, Y={y_i}, Z={z_i}？

我想用狄拉克增量获得一般表达式，这样我就不需要单独处理所有选项。如果这是一个选项，我更喜欢爱因斯坦符号。

感谢期待

Answer 1

在这里回复您的评论。也许这可以满足您的要求：

In [12]: X = IndexedBase('X')                                                                                                                  

In [13]: i, j = symbols('i, j')                                                                                                                

In [14]: X[i].diff(X[j])                                                                                                                       
Out[14]: 
δ   
 i,j