feat: use hp.map2alm_lsq for lmax higher than 1.5 nside

galsci · Sep 24, 2022 · 07f0076 · 07f0076
1 parent 86f4a79
commit 07f0076
Showing 1 changed file with 10 additions and 3 deletions.
diff --git a/pysm3/models/template.py b/pysm3/models/template.py
@@ -184,9 +184,16 @@ def apply_smoothing_and_coord_transform(
         if input_alm:
             alm = input_map.copy()
         else:
-            alm = hp.map2alm(
-                input_map, lmax=lmax, use_pixel_weights=True if nside > 16 else False
-            )
+            if lmax <= 1.5 * nside:
+                alm = hp.map2alm(
+                    input_map,
+                    lmax=lmax,
+                    use_pixel_weights=True if nside > 16 else False,
+                )
+            else:
+                alm, error, n_iter = hp.map2alm_lsq(
+                    input_map, lmax=lmax, mmax=lmax, tol=1e-6, maxiter=100
+                )
         if fwhm is not None:
             hp.smoothalm(alm, fwhm=fwhm.to_value(u.rad), inplace=True, pol=True)
         if rot is not None: