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