L-BFGS алгоритм

[allchemist]

Здравствуйте.
Пытаюсь реализовать квазиньютоновский алгоритм L-BFGS, но он упорно ведет себя не так, как надо.
Буду очень благодарен, если напишете схему алгоритма, для проверки.
В частности, интересует механизм обновления гессиана.
PS. Книгу Гилла и Мюррея читал.