以更快的方式使用healpy自旋变换

Question

我使用的代码需要使用map2alm和alm2map（主要是后者）在healpy中进行频繁的spin-2变换。然而，当我越来越高的时候，这种转变不可避免地会变慢。

我正在使用openmp的healpy，并且由于需要大量的变换，这仍然需要很长时间才能完成一项工作。有没有进一步的方法来尝试加速可能适用于这种情况的转换？例如，因为我只对部分天空感兴趣，有没有办法只使用alm2map重建指定像素而不是整个天空的值，这会更快吗？

感谢您的帮助

Answer 1

I don't think healpy can sped this up, but you might be interested in checking out the other Cl estimators, in case you're interested in the power spectra and not only the alm.

For the full sky (in HEALPix):

• PolSpice
• XPol

For smaller fields:

• Poker

For reasons described here, the estimation of the power spectra is not straightforward for a partially masked sky. HEALPy does not correct for this; the other packages I linked above do.