Мені потрібно вирішити ту саму розріджену лінійну систему (від 300x300 до 1000x1000) з багатьма правими сторонами (від 300 до 1000). На додаток до цієї першої проблеми, я також хотів би вирішити різні системи, але з тими ж ненульовими елементами (просто різними значеннями), тобто багато рідкісних систем з постійною схемою розрідженості. Мої матриці невизначені.
Виконання факторизації та ініціалізації не важливо, але ефективність етапу вирішення є. В даний час я розглядаю PaStiX або Umfpack, і, мабуть, я пограю з Petsc (який підтримує обидва рішення) Чи є бібліотеки, здатні скористатися моїми конкретними потребами (векторизація, багатопотокова передача) чи варто покладатися на загальні вирішувачі, і може, їх трохи змінити для моїх потреб?
Що робити, якщо розріджена матриця більша, до ?