Я написал код, который решает большие системы PDE, используя некоторые методы дискретизации, которые в основном включают решение больших, разреженных систем Ax = b много раз на каждом шаге по времени.
В настоящее время я использую решатель PARDISO (избиблиотека intel MKL), которая является прямой LU-факторизацией A для решения системы.Я хотел бы сравнить этот метод с использованием итерационных решателей (которые при использовании предварительных преобразователей могли бы работать лучше, так как я мог бы использовать один и тот же предварительный обработчик на многих временных шагах, если моя матрица Якоби не сильно меняется).
Тогда мой вопрос: какую библиотеку вы предлагаете для разреженных итерационных решателей в фортране?Я нашел один (SLATEC), который был написан в 1993 году, поэтому мне интересно, есть ли что-то более производительное, которое было написано совсем недавно?
Спасибо :))