Add fabs to prevent negative Spl2 and Smi2
Description
This adds a small patch for IMRPhenomXP(HM) to catch a pathology resulting in NaNs that can occur when the in-plane spins are extremely small (~ numerical precision) leading to negative values of Smi2. Flagged up by @cecilio.garcia-quiros.
API Changes and Justification
Backwards Compatible Changes
-
This change introduces no API changes -
This change adds new API calls
Backwards Incompatible Changes
-
This change modifies an existing API -
This change removes an existing API
If any of the Backwards Incompatible check boxes are ticked please provide a justification why this change is necessary and why it needs to be done in a backwards incompatible way.
Review Status
A regression test has been performed with the corner plot . Mismatches agree to ~ numerical precision. Output from regression test attached here.
Other notes:
- Master branch corresponds to git commit 1016f628 (~ 5th Nov 2020) here:
- Fixed branch corresponds to git commit 73260889 (branched from current master ~ 13th Nov) here:
https://git.ligo.org/waveforms/reviews/lalsuite/-/tree/phenomxphm_patch
Edited by Maria Haney