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

Minor fixes to noise estimation work #730

Merged
merged 1 commit into from
Jan 25, 2024

Conversation

tskisner
Copy link
Member

Some small fixes, discovered with unit tests and batch jobs with 8 or 16 processes per group:

  • When duplicating an observation, also duplicate per-detector flags.

  • When redistributing an observation, reset per-detector flags before setting.

  • Make the common mode removal an option (default True) prior to noise estimation. This is useful if the data has already had the common mode removed prior to this operator.

  • When high-pass filtering data that is distributed by time slices there may be some chunks that are fully flagged. This should not be an error. Instead, the flagged timestream is set to zero.

  • When communicating overlaps, ensure that the span used (lagmax

    • half_average) is less than the number of samples.

with 8 or 16 processes per group:

* When duplicating an observation, also duplicate per-detector
  flags.

* When redistributing an observation, reset per-detector flags
  before setting.

* Make the common mode removal an option (default True) prior
  to noise estimation.  This is useful if the data has already
  had the common mode removed prior to this operator.

* When high-pass filtering data that is distributed by time slices
  there may be some chunks that are fully flagged.  This should
  not be an error.  Instead, the flagged timestream is set to zero.

* When communicating overlaps, ensure that the span used (lagmax
  + half_average) is less than the number of samples.
Copy link
Member

@keskitalo keskitalo left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good!

@tskisner tskisner merged commit 112d388 into noise_estimation_fix Jan 25, 2024
@tskisner tskisner deleted the noise_estimation_fix_fix branch January 25, 2024 22:47
keskitalo added a commit that referenced this pull request Jan 26, 2024
* Fix parallelization problems in noise estimation

* Small fixes for self.view=None

* Some small fixes, discovered with unit tests and batch jobs (#730)

with 8 or 16 processes per group:

* When duplicating an observation, also duplicate per-detector
  flags.

* When redistributing an observation, reset per-detector flags
  before setting.

* Make the common mode removal an option (default True) prior
  to noise estimation.  This is useful if the data has already
  had the common mode removed prior to this operator.

* When high-pass filtering data that is distributed by time slices
  there may be some chunks that are fully flagged.  This should
  not be an error.  Instead, the flagged timestream is set to zero.

* When communicating overlaps, ensure that the span used (lagmax
  + half_average) is less than the number of samples.

* Ignore flagged detectors in poly2d fit

---------

Co-authored-by: Theodore Kisner <tskisner@users.noreply.github.com>
Co-authored-by: Ted Kisner <tskisner.public@gmail.com>
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