-
Notifications
You must be signed in to change notification settings - Fork 19
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 support for single detector maps #690
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looks fine. One small improvement noted if you wish.
job_ops.mapmaker.apply(obs_data) | ||
# See if user wants separate detector maps | ||
if hasattr(otherargs, "detmaps") and otherargs.detmaps: | ||
my_dets = data.all_local_detectors() |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This might cause problems if some detectors are cut in all observations. You could replace this line with
my_dets = data.all_local_detectors(flagmask=defaults.det_mask_invalid)
to get the global detectors which are not cut by major issues.
if data.comm.comm_world is None: | ||
all_dets = my_dets | ||
else: | ||
from matplotlib.cbook import flatten |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Interesting, I had never seen this method! Definitely useful for the typical use case of calling gather
/ allgather
on python lists, resulting in a list of lists. We could consider doing a namespace import as toast.mpi.flatten
and use it more generally.
This PR requires hpc4cmb/toast#743 to be merged first. |
89ddb21
to
79a1bc6
Compare
This PR is superseded by #877 |
No description provided.