Scipy Hessian с перекрестным произведением Numpy ⇐ Python

[Anonymous]

Мне нужно вычислить гессиан функции, включающей перекрестное произведение, в Python. В качестве эксперимента рассматриваю простую функцию

Код: Выделить всё

def func(x):
return np.sum(np.cross(x[0:3], x[3:6]))

и вычислите гессиан с помощью

Код: Выделить всё

x0 = np.array([1,2,3,4,5,6])
H = scipy.differentiate.hessian(func, x0)

Однако это не удается из-за ошибки.

ValueError: несовместимые измерения для перекрестного произведения
(размер должен быть 2 или 3)

Я предполагаю, что проблема связана с этой строкой из документации scipy.dependentiate.hessian:

аргумент f должен быть векторизован, чтобы принять массив формы (m, ...). Первая ось представляет m входных данных f; остальные оси, обозначенные эллипсами, предназначены для вычисления функции по нескольким абсциссам за один вызов.

но у меня нет знаний Python, чтобы это исправить. Может ли кто-нибудь направить меня в правильном направлении?

Подробнее здесь: https://stackoverflow.com/questions/798 ... ss-product