Использование healpy spin трансформаций быстрее

Question

Я использую код, который должен часто преобразовывать спин-2 в healpy, используя как map2alm, так и alm2map (преимущественно последний). Однако, когда я перехожу на все более и более высокие стороны, это преобразование неизбежно становится медленнее.

Я использую healpy с openmp, и из-за огромного объема необходимых преобразований это все еще приводит к длительному завершению работы. Есть ли еще какой-нибудь способ попытаться ускорить преобразование, которое может быть применимо к этой ситуации? Например, поскольку меня интересует только часть неба, есть ли способ восстановить значения только в указанных пикселях, а не во всем небе при использовании alm2map, и будет ли это быстрее?

Спасибо за любую помощь

Answer 1

Я не думаю healpy может ускорить это, но вам может быть интересно проверить другие оценки Cl, если вас интересуют спектры мощности, а не только alm,

Для полного неба (в HEALPix):

• PolSpice
• Xpol

Для небольших полей:

• покер

По причинам, описанным здесь, оценка спектров мощности не является прямой для частично замаскированного неба. HEALPy не исправляет это; другие пакеты, которые я связал выше, делают.