search

Libsharp / libsharp

Library for fast spherical harmonic transforms, see http://arxiv.org/abs/1303.4945
GNU General Public License v2.0
24 stars 17 forks source link

Equivalent of HEALPix iter #19

Closed zonca closed 5 years ago

zonca commented 5 years ago

map2alm in healpy/HEALPix accepts a iter argument, is there anything equivalent in libsharp/MPI?

mreineck commented 5 years ago

Libsharp doesn't have any equivalent at the moment. As far as I understand, the iterative solver only makes sense for Healpix. Also, the library would have to allocate space for a temporary set of a_lm or a map.

Implementing this on the user side is not hard fortunately. For an exmple, see map2alm_iter in Healpix_cxx/alm_healpix_tools.cc. On the high-level language side, this can be done with about 10 lines. For libsharp it would be much harder.