Новый метод быстрых мультиполей ускорил расчеты гравитации в симуляциях RAMSES
Исследователи разработали новый масштабируемый решатель уравнения Пуассона для кода RAMSES, основанный на методе быстрых мультиполей (FMM). Алгоритм обладает линейной вычислительной сложностью и предназначен для ускорения гравитационных расчетов в астрофизических симуляциях. Работа опубликована на arXiv.
Код RAMSES широко используется для моделирования гидродинамики и гравитации при изучении формирования галактик и крупномасштабной структуры Вселенной. Ранее в нем применялся многосеточный метод (MG), итеративный алгоритм, требующий многократных V-циклов на иерархии сеток. FMM предлагает альтернативный подход.
В методе FMM строится иерархия сеток поверх существующей декартовой сетки. Вычисления производятся за один восходящий проход, в ходе которого накапливаются и сдвигаются мультипольные разложения, и один нисходящий проход для распространения локальных разложений. Это исключает необходимость итераций.
Численные тесты показали, что FMM достигает точности, сопоставимой с MG для гладких потенциалов. Особое преимущество метода проявляется в задачах с изолированными граничными условиями, поскольку он не требует приближенных граничных условий Дирихле, необходимых для MG-схем.
Хотя теоретически FMM требует примерно в 30 раз больше операций с плавающей точкой, чем MG, его более высокая арифметическая интенсивность обеспечивает сопоставимую производительность. При этом FMM демонстрирует лучшую масштабируемость на параллельных системах, что делает его перспективным для крупномасштабных симуляций.
Новый решатель может ускорить моделирование в астрофизике, особенно для задач с изолированными границами, таких как формирование отдельных галактик или скоплений. Работа выполнена в рамках развития методов вычислительной астрофизики.
