假设我有以下代码
import numpy as np
f = lambda x,y: (np.sum(x) + np.sum(y))**2
x = np.array([1,2,3])
y = np.array([4,5,6])
df_dx
df_dy
df2_dx2
df2_dxdy
...
有一种快速计算这种函数的所有衍生物（单一和混合）的方法吗？模块应该在数组级执行经典的有限差分技术，即将元素h= tol元素添加到数组变量（取决于导数），计算函数并除以h。
（我的实际情况要复杂得多，因为它涉及来自DLL的数组值函数我无法修改......变量的数量是任意的，请不要专注于这个特定的玩具示例）< / p>