Задача нахождения собственных значений трехмерной гессенской матрицы - PullRequest
0 голосов
/ 09 января 2019

Я применяю эту функцию (https://stackoverflow.com/a/31207520/8948867) к трехмерному массиву (5x5x5), чтобы получить гексановую матрицу 3x3 этого объема. Форма полученной гессианской матрицы составляет 3x3x5x5x5.

Моя цель состоит в том, чтобы вычислить 3 собственных значения двумерной гессенской матрицы, но, поскольку моя гессенская матрица имеет более двух измерений (3x3x5x5x5), форма вычисленных собственных значений составляет 3x3x5x5.

Я не знаю, заключается ли проблема в гессианской матрице или в функции для оценки собственных значений, или я, вероятно, пропускаю некоторые промежуточные шаги.

Это код, который я использую (гессиан - это функция, показанная в ссылке выше):

hss = hessian (arr)
from numpy import linalg as LA
e_val, e_vec = LA.eigh (hss)
...