Я использую scipy.optimize.minimize, чтобы свести к минимуму простую функцию правдоподобия журнала. Матрица Гессиана, похоже, плохо себя ведет. import scipy.optimize as op
def lnlike(theta, n, bhat,
С минимизатором L-BFGS-B в scipy, можно ли получить приблизительный обратный Гессиан, который вычисляется внутри? Имея его в неявной факторизованной форме, так что можно вычислить произвольные обратны
Я хочу найти хребты для данного изображения. (Хребты не ребра!) Примером может служить, как на картинке ниже Я думаю, что матрица Hessian будет работать интуитивно. Таким образом, я жестко закодировал