-
Notifications
You must be signed in to change notification settings - Fork 3
ESIGMA Backward integration #2
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
base: master
Are you sure you want to change the base?
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.
Copilot reviewed 2 out of 2 changed files in this pull request and generated no comments.
Comments suppressed due to low confidence (1)
esigmapy/generator.py:213
- The error message still references 'orbital_var_names' even though the valid names are now stored in 'all_orbital_var_names'. Please update the error message to reflect the correct variable.
f"{name} is not a valid orbital variable name. Available orbital variable names are: {orbital_var_names}."
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.
Copilot reviewed 2 out of 2 changed files in this pull request and generated 1 comment.
Comments suppressed due to low confidence (1)
esigmapy/generator.py:221
- The backward integration branch computes orbital parameters using ls.SimInspiralESIGMADynamicsBackwardInTime and then immediately calls ls.SimInspiralENIGMADynamics with the updated parameters. Please add a clarifying comment to document the rationale behind this two-step integration process.
elif f_ref > f_lower:
| mean_anomaly = l.data.data[-1] | ||
| f_start = f_lower | ||
| elif f_ref < f_lower: | ||
| itime = time.perf_counter() |
Copilot
AI
Apr 13, 2025
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.
The variable 'itime' is assigned multiple times across different branches. If these timing measurements are important, consider consolidating them or adding comments to clarify their distinct purposes.
| itime = time.perf_counter() | |
| timing_data["f_ref_less_start"] = time.perf_counter() # Start timing when f_ref < f_lower |
2af8da9 to
161b008
Compare
When generating ESIGMA waveforms, this PR will enable us to specify the initial conditions at a reference frequency, which can be larger than the starting GW frequency. Currently we only support the cases where the reference frequency is equal to or lower than the starting GW frequency.