Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

add demo_QE notebook (and use NamedTuple for exps) #6

Open
wants to merge 2 commits into
base: master
Choose a base branch
from

Conversation

cmbant
Copy link

@cmbant cmbant commented Nov 30, 2021

If you want a basic QE demo, something like this seems to work (and shows as expected that skewness of QU is much lower than T recon)

@cmbant
Copy link
Author

cmbant commented Nov 30, 2021

Might want to tidy up for consistent use of gradient responses etc

@carronj
Copy link
Owner

carronj commented Nov 30, 2021

Thanks. I believe you adding the 'name' to get_config will break some codes. But can include something similar

@cmbant
Copy link
Author

cmbant commented Nov 30, 2021

Yes, sorry didn't test with anything else. Just wanted some way to play with expt config in notebook without having to modify the list of recognised definitions in the source.

@cmbant
Copy link
Author

cmbant commented Nov 30, 2021

Btw, weird stuff happens if you change LDres, HDres, eg.. kappa looks crazy if I set HDres = 9

@carronj
Copy link
Owner

carronj commented Nov 30, 2021

It makes little sense to LDres (~the data resolution) higher than HDres (~the sky resolution), but otherwise (e.g. both at 9) should work?

@cmbant
Copy link
Author

cmbant commented Nov 30, 2021

No, that looks crazy too. Maybe I misunderstood what they are - from the get_maps_lib doc string I thought HDres was used to change the physical patch size, and LDres the number of pixels in the patch.

@carronj
Copy link
Owner

carronj commented Dec 3, 2021

Thanks. with HD 9 it looks like the response at L close to the high-L limit 2 * lmax_cmb is so tiny and has(I guess) numerical noise that gets blown up by inverting this to normalize phi. These super-high Ls are useless modes anyways, so simplest is to truncate the ellmax of lib_qlm to something reasonable (like lmax_cmb)

@cmbant
Copy link
Author

cmbant commented Dec 3, 2021

I was trying to just change the patch size to something smaller, with resolution, lmax etc all the same. It seems there should be three numbers, the data res, sky res, and pixel number, what is the intended way to do this?

@carronj
Copy link
Owner

carronj commented Dec 3, 2021

What you wrote is correct, e.g. (9 10 ) gives the same as (10, 11) but on a box with sides smaller by a factor of two. (to avoid the kappa real-space you had, just better not use the lensing L ~ 6000 in remapping to real space)
(The sky resolution is fixed to 0.74 amin using the sim libs in this way, but one can define any othes adapting the sim lib definitions )

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants