С успехом использовал пакет SPARSKIT Дж.Саада
http://www-users.cs.umn.edu/~saad/softw ... rskit.html
Написан на фортране 77 и работает, как швейцарские часы.
Реализует в основном метод сопряженных градиентов и его модификации
с предобуславливателями. Не требует дополнительных
библиотек типа blas и т.д.
Очень понятный интерфейс. Главное, все работает.
Самое забавное, что я перевел исходник этого пакета на язык С
с помощью известной утилиты f2c.exe и тоже все работет замечательно.
Many thank to Saad! Good thing.
В списке отсутствует еще один пакет, который я тоже
успешно использовал: UMFPACK Тима Девиса
http://www.cise.ufl.edu/~davis/
Пакет сразу написан на С и реализует прямые методы, которые
работают как кувалда! Тоже очень понятный интерфейс, но
требует дополнительную библиотеку blas.
Но это небольшое затруднение. Я ее всю собрал и очень доволен.
Этот пакет рекомендую для задач, с не очень большим числом неизвестных.
Так, например, для решения несжимаемых ур. Навье-Стокса я ее
использовал на шаге решения ур. Пуассона. С RAM 1 Гиг максимальная
трехмрная сетка была где-то от 70 до 100 тыс ячеек. Больше не влезало из-за
очень больших L и U матриц (хотя и разреженных).
Для LES рачетов это мало. Нужно до миллиона. Поэтому использовал
итерационные методы, а именно SPARSKIT.
Вобщем, рекомендую.