При моделировании кинетики химических реакций, расчете электронных схем и ряде других важных приложений возникает проблема численного решения задачи Коши для жестких систем обыкновенных дифференциальных уравнений. Современные методы решения жестких систем, как правило, используют обращение матрицы Якоби. В случае достаточно большой размерности дифференциальной задачи эффективность алгоритма интегрирования фактически полностью определяется временем обращения этой матрицы. Поэтому одним из способов повышения эффективности является распараллеливание данной стадии алгоритма.
Здесь реализован L-устойчивый (4,2)-метод [1] третьего порядка точности на многопроцессорных кластерах. Представлены результаты расчетов больших жестких систем для разного числа процессоров с указанием количественных характеристик многопроцессорности (время счета, коэффициент ускорения и т.п.). Распараллеливание вычислений обеспечивалось с помощью языка Fortran-DVM, разработанного в ИПМ им. Келдыша РАН.
[1] Новиков Е.А., Шитов Ю.А., Шокин Ю.И. Одношаговые безитерационные методы решения жестких систем / ДАН СССР, т.301, N 6, 1988. - C. 1310-1314.
Примечание. Тезисы докладов публикуются в авторской редакции
Ваши комментарии Обратная связь |
[Головная страница] [Конференции] |
© 1996-2000, Институт вычислительных технологий СО РАН, Новосибирск
© 1996-2000, Сибирское отделение Российской академии наук, Новосибирск