CosmoCoffee

[Zhen Hou]

Dear Antony,

I was trying to run lenspix with polarization turned on for high nside values, e.g., nside=4096 and 8192. However, I find all the elements of the output Q and U map are '-NaN' if I set lmax larger than 7000. In fact, for lmax lower than 7000, I only tested lmax=6500 that the output looked fine. The T map is always good for large lmax.

I have checked the input Cls which contained ell from 2 to 15100 with delta_ell=1 for [2,10000] and delta_ell=100 for [10000, 15100]. I don't think there is problem in the input Cls.

I am using lenspix version April 2014 and Healpix 3.11, both compiled by intel fortran 14.0.2.

Here are relevant lines in my params.ini

nside = 4096
lmax = 7000
cls_file = planck_lensing_wp_highL_bestFit_20130627_lenspotentialCls.dat
lens_method = 1
interp_factor = 0.375
interp_method = 1
mpi_division_method = 3
want_pol = T
want_map = T
want_alm = F
rand_seed = 320

I am not sure if any one met this problem before. Could you help to look into it? Thank you!

Cheers,
Zhen

[Pavel Motloch]

Hi Zhen,

my guess is your input Cls are flawed as they do not satisfy a physicality condition C_TT*C_EE > C_TE^2. Then Lenspix generates NaN a_lm, which then propagate into the maps.

This can happen if in CAMB you choose scalar lmax > 8000 (which is the lmax of the provided high precision interpolation file) due to complicated CAMB interpolation scheme.

Can hack this in Lenspix code which generates alms, *_Sim or so.

Pavel