вывод уравнений для облаков точек в симптах или альтернативах

Question

Предположим, у вас есть облако точек (хранится в R, размерная матрица 3xN) в трехмерной ортогональной системе координат, и вы хотите получить аналитические выражения для сложных функций R (включая дифференцирование, суммирование и т. Д. c) .

В качестве простого примера:

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

Есть ли способ сделать это симпати (или что-нибудь еще ...), например, обобщив x,y,z до X={x_i}, Y={y_i}, Z={z_i}?

Я хотел бы получить общее выражение с дельтой Dira c, чтобы мне не нужно было обрабатывать все параметры по отдельности. Если это вариант, я бы предпочел обозначение Эйнштейна.

Спасибо в ожидании

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