Новый алгоритм ускоряет расчёт матриц для анализа реликтового излучения

Авторы взяли не саму физику реликтового излучения, а один из самых тяжёлых вычислительных кусочков в его анализе — построение матриц, которые описывают, как маска неба смешивает разные угловые масштабы сигнала. Обычно такие матрицы считают рекурсивными схемами, то есть шаг за шагом пересчитывают много промежуточных значений. Здесь же они нашли, какие именно комбинации нужны на практике, и для них заранее посчитали небольшой набор опорных значений, а всё остальное аккуратно восстановили из них. Благодаря этому расчёт стал быстрее, но формулы при этом не пришлось упрощать или приближать.

Главная идея в том, что авторы использовали особую структуру коэффициентов Вигнера–3j — это математические объекты, которые часто появляются, когда связывают разные угловые моды на сфере. Важно не просто то, что метод быстрый, а то, что он остаётся точным: на реалистичных масках неба он даёт те же численные ответы, что и обычные схемы. То есть ускорение не покупается ценой потери качества, что для CMB-анализа особенно важно, потому что здесь даже небольшие ошибки в матрицах потом могут разойтись по всей цепочке расчётов.

Ещё один полезный момент: метод хорошо распараллеливается. Его можно гонять и на обычных процессорах, и на видеокартах, причём GPU добавляет ещё один большой рывок по скорости. Авторы отдельно показывают, что подход удобен не только для полной матрицы, но и для задач, где нужен лишь ограниченный набор мод вместо всей таблицы целиком — например, для полосно-ограниченных матриц связи и аналитических оценок ковариации, то есть для расчётов, где важна не вся картина, а только нужный кусок.