Добрый день!
Есть такая задача - есть матрица частных производных

(большая разреженная). И необходимо многократно пересчитывать разложение Холесского

Проблема в том, что иногда матрица

меняется не только численно, но и к ней в конец добавляются/удаляются строки. Хочется не пересчитывать полностью при этом уже посчитанную на итерации матрицу

, а пересчитать только некоторые элементы.
Т.е., например

размерностью 10000х5000, заполненность 0.001. Размерность

соответственно 5000х5000. Посчитали

. После этого надо добавить, например, 10 строчек к

(остальная часть

не меняется) и пересчитать

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