|
baleog |
|
|
|
Добрый день!
Столкнулся с тем, что в доставшихся в "наследство" расчетах выполняется процедура т.н. "динамического разложения", которая по сути упорядочивает строки матрицы Гессе в непонятном для меня порядке.
По факту такое упорядочивание дает значительный прирост в скорости решения системы уравнений. Подскажите, пожалуйста, из-за чего происходит такой эффект? Почему изменение порядка строк в матрице Гессе ускоряет расчет?
Спасибо!
|
|
|
|
 |
|
мат-ламер |
|
|
|
Про динамическое разложение не знаю. Но может матрица разреженая ленточная. И переупорядочивание может уменьшить ширину ленты.
|
|
|
|
|
|
baleog |
|
|
|
Да, проблему уже решил сам, но мат-ламер
прав - дело было именно в работе с разреженными матрицами! На этапе разложения индексные операции слишком дорого стоили
|
|
|
|
|
|
baleog |
|
|
|
Возник еще один вопрос. При построении матрицы Гессе используется приближение
h = j' * diag(r) * j
где j - матрица якоби, ' - транспонирование, r - вектор весов.
Подскажите, почему гессиан можно заменить на такую вот аппроксимацию?
|
|
|
|
|
