2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




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

 
 
 
 Posted automatically
Сообщение07.02.2017, 14:55 
 i  Тема перемещена из форума «Околонаучный софт» в форум «Computer Science»
Причина переноса: тематика.

 
 
 [ Сообщений: 2 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group