Nf subtraction bug (#201)

LSSTDESC · Jul 17, 2024 · 76c3d41 · 76c3d41
1 parent 39a9894
commit 76c3d41
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/pymaster/workspaces.py b/pymaster/workspaces.py
@@ -734,8 +734,9 @@ def compute_coupled_cell(f1, f2):
 
     cls = np.array([[hp.alm2cl(a1, a2, lmax=lmax)
                      for a2 in alm2] for a1 in alm1])
-    for i in range(len(alm1)):
-        cls[i, i, :] -= Nf
+    if Nf != 0:
+        for i in range(len(alm1)):
+            cls[i, i, :] -= Nf
     cls = cls.reshape([ncl, lmax+1])
     return cls