Разреженная несимметричная ленточная матрица

[maxym]

Здраствуйте!
Посоветуйте выбрать алгоритм для разреженной матрицы.
Для метода Ньютона-Рафсона нужно решить СЛАУ с матрицей Якоби(1000*1000). Матрица структурно-симметричная( но не абсолютно симметричная)
Она ленточная. Ширина ленты - около 15 ( т. е. порядка 1% от длины диагонали), на ленте будет много нулей.
Так или иначе, ее нужно обращать.
Смотрел LU-декомпозицию, но она не учитывает разреженности.
QR-декомпозиция - сложная в реализации, но на выходе, имхо, получим тоже что и LU.
В инете куча алгоритмов (Алгоритм Катхилла — Макки, Холецкого например), но они только для симметричных матриц.
Смотрел Писсанецки, читал Брамеллера.
Там говорится, что нужно найти приемлемый метод упорядоченного исключения Гаусса именно для моей матрицы. Но как в самом Брамеллере, так и в гугле я ничего по этому поводу не нашел.
Возможно ленту этой матрицы можно как то сузить, а потом обратить(или разложить и потом обратить например обратным ходом).

[Karan]

Тема перемещена из форума «Околонаучный софт» в форум «Computer Science»
Причина переноса: тематика.