CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

lenspix NaN in output Q and U map when lmax > 7000
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Zhen Hou



Joined: 02 Mar 2009
Posts: 6
Affiliation: KICP, University of Chicago

PostPosted: October 16 2014  Reply with quote

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
Back to top
View user's profile [ Hidden ]
Pavel Motloch



Joined: 14 Oct 2016
Posts: 6
Affiliation: U Chicago

PostPosted: March 16 2017  Reply with quote

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
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.