catalog2alm parallel

Author: damonge

pymaster/utils.py

@@ -849,11 +849,12 @@ def _map2alm_healpy(map, spin, sht_info, alm_info):
           'healpy': _alm2map_healpy}
 
 
-def _alm2catalog_ducc0(alms, positions, spin, lmax):
+def _alm2catalog_ducc0(alms, positions, spin, lmax, nthreads=0):
     alms = np.atleast_2d(alms)
     values = ducc0.sht.synthesis_general(alm=alms, spin=spin,
                                          lmax=lmax, loc=positions.T,
-                                         epsilon=1E-5)
+                                         epsilon=1E-5,
+                                         nthreads=nthreads)
     return values