Работаю с разреженными матрицами (МКЕ) и приходится использовать разложение Холецкого, где нужно переупорядочить матрицу методом минимальной степени.
Использую библиотеку NL
http://www.uic.nnov.ru/~zny/nl/index.html
Если кто работал с библиотекой или просто знает в чем дело подскажите, почему если не использовать упорядочение - решение по Холецкому системы UTDUxB=X есть, стоит переупорядочить - выводит огромные числа.
система уравнений 24х24 - матрица жесткости одного 8-ми узлового кирпичика. Граничные условия обеспечены.
Вектор перестановок выглядит как 23 22 21 ... 0
переставляет строки и столбцы вроде правильно
Или подскажите библиотеку на Си, где реализованы эти алгоритмы работы с спарс матрицами
16.12.2008, 01:41 