simonsobs / pixell

A rectangular pixel map manipulation and harmonic analysis library derived from Sigurd Naess' enlib.
Other
42 stars 32 forks source link

Faster almxfl? #273

Open msyriac opened 1 month ago

msyriac commented 1 month ago

optweight.alm_c_utils.lmul is allegedly 80% faster than pixell.curvedsky.almxfl https://github.com/AdriJD/optweight/blob/0b286eba73d712b4b3a54ebcdd0371de495432d0/cython/alm_c_utils.pyx#L157

@amaurea @adrijd keeping this open until we understand this / incorporate Adri's version into pixell

Re lmul, yes:

image

Perhaps worth opening a pixell issue for this, I agree we should streamline all of these dependencies/duplicated functionality across our core codes.

Originally posted by @zatkins2 in https://github.com/simonsobs/mnms/issues/9#issuecomment-2397513461

amaurea commented 1 month ago

What shape arrays is this for? Do you have a test case?