fix large-scale ell to k conversion (h/t Boris Bolliet)

simonsobs · Jan 29, 2021 · 6828bab · 6828bab
1 parent 8b285f4
commit 6828bab
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/hmvec/cosmology.py b/hmvec/cosmology.py
@@ -398,7 +398,7 @@ def limber_integral(ells,zs,ks,Pzks,gzs,Wz1s,Wz2s,hzs,chis):
          f = interp1d(ks,Pzks[0],bounds_error=True)
     Cells = np.zeros(ells.shape)
     for i,ell in enumerate(ells):
-        kevals = ell/chis
+        kevals = (ell+0.5)/chis
         if zs.size>1:
             # hack suggested in https://stackoverflow.com/questions/47087109/evaluate-the-output-from-scipy-2d-interpolation-along-a-curve
             # to get around scipy.interpolate limitations