упорядочивание строк в матрице гессе

[baleog]

Добрый день!

Столкнулся с тем, что в доставшихся в "наследство" расчетах выполняется процедура т.н. "динамического разложения", которая по сути упорядочивает строки матрицы Гессе в непонятном для меня порядке.

По факту такое упорядочивание дает значительный прирост в скорости решения системы уравнений. Подскажите, пожалуйста, из-за чего происходит такой эффект? Почему изменение порядка строк в матрице Гессе ускоряет расчет?


Спасибо!

[мат-ламер]

Про динамическое разложение не знаю. Но может матрица разреженая ленточная. И переупорядочивание может уменьшить ширину ленты.

[baleog]

Да, проблему уже решил сам, но мат-ламер
прав - дело было именно в работе с разреженными матрицами! На этапе разложения индексные операции слишком дорого стоили

[baleog]

Возник еще один вопрос. При построении матрицы Гессе используется приближение

h = j' * diag(r) * j

где j - матрица якоби, ' - транспонирование, r - вектор весов.

Подскажите, почему гессиан можно заменить на такую вот аппроксимацию?