Чтобы решить проблему Пуассона в 2D, у меня есть очень разреженная матрица с полосами, состоящая из трех полос (https://people.eecs.berkeley.edu/~demmel/cs267/lecture24/DiscretePoisson2.gif).
Используя линейный решатель ScaLAPACK, я могу получить из n ^ 4 матрицы до ~ n ^ 3 с использованием полосового хранилища. Однако для большой области с сотнями нулей между полосами более подходящим будет разреженный решатель.
Поэтому я ищу библиотеку разреженных линейных решателей , совместимую с MPI . Я был бы признателен за несколько альтернативных библиотек, так как для моего приложения можно использовать только предустановленное программное обеспечение. Заранее спасибо за любую помощь:)