CAR debugged

LSSTDESC · Dec 27, 2023 · 0022a3f · 0022a3f
1 parent 33ecb7e
commit 0022a3f
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/pymaster/utils.py b/pymaster/utils.py
@@ -64,7 +64,8 @@ def from_rectpix(cls, n_theta, theta_min, d_theta,
 
         # CC weights according to ducc0
         i_theta_0 = int(theta_min/d_theta+0.5)
-        nth_full = int(np.pi/d_theta+0.5)
+        # +1 below because both poles get a ring
+        nth_full = int(np.pi/d_theta+0.5)+1
         w = ducc0.sht.experimental.get_gridweights('CC', nth_full)
         w *= d_phi/(2*np.pi)
         weight_th = w[i_theta_0:i_theta_0+n_theta]
@@ -269,7 +270,7 @@ def reform_map(self, maps):
         if self.flip_ph:
             maps = maps[..., :, ::-1]
         # Flatten last two dimensions
-        maps = maps.reshape(maps.shape[::-2]+(self.npix,))
+        maps = maps.reshape(maps.shape[:-2]+(self.npix,))
         return maps
 
     def _map_compatible(self, mp):