Новый алгоритм FMM ускорил моделирование космоса в коде RAMSES
Международная группа исследователей разработала новый решатель уравнения Пуассона для популярного астрофизического кода RAMSES. Алгоритм основан на методе быстрых мультиполей (FMM) и обеспечивает линейную вычислительную сложность O(N), что критически важно для моделирования крупномасштабных структур Вселенной.
Код RAMSES широко используется для гидродинамических и гравитационных расчетов в космологии. Новый решатель, построенный на предыдущем унигрид-алгоритме Ли и Тейсье (2026), добавляет поддержку адаптивного сеточного уточнения (AMR) и адаптивного шага по времени (ATS). Это позволяет эффективно распределять вычислительные ресурсы между областями с разной детализацией.
В работе представлены несколько новшеств: использование множества FMM-деревьев для каждого уровня AMR, объединенное дерево для всех активных уровней для оптимизации поиска соседей, а также концепция near-field для обеспечения симметрии сил на границах уровней.
При тестировании на разнообразных задачах новый метод показал превосходное согласие с эталонным многосеточным (MG) решателем на уровне процентов. Однако FMM продемонстрировал лучшие свойства сохранения импульса на границах между крупной и мелкой сетками.
Несмотря на дополнительные вычислительные затраты из-за пространственно-временной адаптивности, FMM показал более высокую масштабируемость по сравнению с MG. Разрыв в производительности увеличивается с ростом размеров задачи, что делает алгоритм особенно привлекательным для симуляций с высоким разрешением.
Разработка может ускорить моделирование формирования галактик, скоплений и других космологических процессов, требующих огромных вычислительных мощностей. Алгоритм уже интегрирован в RAMSES и доступен для сообщества.
