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.

Sign up for GitHub

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

Synchrotron with local-spectra modulation #152

Merged
merged 18 commits into from
Feb 18, 2023
Merged

Synchrotron with local-spectra modulation #152

Changes from 1 commit
Commits
Show all changes
18 commits
Select commit Hold shift + click to select a range
9db4fcb
update synchrotron notebook, also produces artifacts
zonca Feb 15, 2023
ec9e6c7
add modulation to beta
zonca Feb 15, 2023
e4ab74c
renamed input file
zonca Feb 15, 2023
d7b6567
Revert "add modulation to beta"
zonca Feb 16, 2023
b4576fc
use modulation from previous beta release for beta dust
zonca Feb 16, 2023
048f55d
fix: ell max for beta products was too low
zonca Feb 16, 2023
1cc803b
remove numpy type from filename
zonca Feb 17, 2023
0f24f2f
change filenames and use modulation from prev release
zonca Feb 17, 2023
a1a3579
Clean execution, ell_max now 2.5 nside capped at 16384
zonca Feb 17, 2023
c2825ee
add units to modulation alms
zonca Feb 17, 2023
8e5047f
implemented separate maps for modulation of beta
zonca Feb 17, 2023
11bdeb4
galplane fix for synchrotron
zonca Feb 17, 2023
732b70d
fix filename
zonca Feb 17, 2023
8488df5
galplane fix filename was wrong
zonca Feb 17, 2023
0a702f9
pass the amplitude modulation in presets for s6
zonca Feb 17, 2023
bf20554
save memory by deleting unneeded objects
zonca Feb 17, 2023
1c93669
fix lmax and names
zonca Feb 17, 2023
c691dfd
other slurm script
zonca Feb 17, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
save memory by deleting unneeded objects
  • Loading branch information
@zonca
zonca committed Feb 17, 2023
commit bf20554b006945b53a2608fe94b47ad10bc5d7e0
43 changes: 23 additions & 20 deletions docs/preprocess-templates/utils_synch_generate_map.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -115,6 +115,15 @@
"map_log_pol_tens_large_scale = hp.alm2map(alm_log_pol_tens_large_scale.astype(np.complex128), nside=output_nside)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"del alm_log_pol_tens_large_scale"
]
},
{
"cell_type": "markdown",
"metadata": {},
Expand All @@ -128,7 +137,7 @@
"metadata": {},
"outputs": [],
"source": [
"modulate_alm = { k:hp.read_alm(output_dir_raw/f\"synch_{k}_modulation_alms_lmax768.fits.gz\").astype(np.complex128) for k in [\"temperature\",\"polarization\"] }"
"modulate_alm = { k:hp.read_alm(output_dir_raw/f\"synch_{k}_modulation_alms_lmax192.fits.gz\").astype(np.complex128) for k in [\"temperature\",\"polarization\"] }"
]
},
{
Expand Down Expand Up @@ -166,6 +175,7 @@
"\n",
"alm_log_pol_tens_small_scale = [hp.almxfl(each, np.ones(output_lmax+1)) for each in alm_log_pol_tens_small_scale]\n",
"map_log_pol_tens_small_scale = hp.alm2map(alm_log_pol_tens_small_scale, nside=output_nside)\n",
"del alm_log_pol_tens_small_scale\n",
"map_log_pol_tens_small_scale[0] *= hp.alm2map(modulate_alm[\"temperature\"], output_nside)\n",
"map_log_pol_tens_small_scale[1:] *= hp.alm2map(modulate_alm[\"polarization\"], output_nside)\n",
"assert np.isnan(map_log_pol_tens_small_scale).sum() == 0"
Expand All @@ -191,6 +201,15 @@
"map_log_pol_tens = map_log_pol_tens_large_scale + map_log_pol_tens_small_scale"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"del map_log_pol_tens_large_scale, map_log_pol_tens_small_scale"
]
},
{
"cell_type": "code",
"execution_count": null,
Expand Down Expand Up @@ -241,16 +260,9 @@
"metadata": {},
"outputs": [],
"source": [
"hp.mollview((output_map - log_pol_tens_to_map(map_log_pol_tens))[0])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"hp.mollview(output_map[0])"
"if output_nside < 4096:\n",
" hp.mollview((output_map - log_pol_tens_to_map(map_log_pol_tens))[0])\n",
" hp.mollview(output_map[0])"
]
},
{
Expand Down Expand Up @@ -278,15 +290,6 @@
"add_metadata([output_dir / f\"synch_template_nside{output_nside}.fits\"], coord=\"G\", unit=\"uK_RJ\", ref_freq=\"23 GHz\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"ls -lh $output_dir"
]
},
{
"cell_type": "code",
"execution_count": null,
Expand Down