From 47b41befa538cd3b51e4b4a1cd3d623e7fa258fd Mon Sep 17 00:00:00 2001
From: Gregory Ashton <gregory.ashton@ligo.org>
Date: Tue, 20 Apr 2021 12:28:33 +0100
Subject: [PATCH] Add errstate ignore to the gw.conversion module
---
bilby/gw/conversion.py | 21 ++++++++++++---------
1 file changed, 12 insertions(+), 9 deletions(-)
diff --git a/bilby/gw/conversion.py b/bilby/gw/conversion.py
index 1ded6605e..5c4577d6b 100644
--- a/bilby/gw/conversion.py
+++ b/bilby/gw/conversion.py
@@ -254,16 +254,18 @@ def convert_to_lal_binary_black_hole_parameters(parameters):
for angle in ['tilt_1', 'tilt_2', 'theta_jn']:
cos_angle = str('cos_' + angle)
if cos_angle in converted_parameters.keys():
- converted_parameters[angle] =\
- np.arccos(converted_parameters[cos_angle])
+ with np.errstate(invalid="ignore"):
+ converted_parameters[angle] =\
+ np.arccos(converted_parameters[cos_angle])
if "delta_phase" in original_keys:
- converted_parameters["phase"] = np.mod(
- converted_parameters["delta_phase"]
- - np.sign(np.cos(converted_parameters["theta_jn"]))
- * converted_parameters["psi"],
- 2 * np.pi
- )
+ with np.errstate(invalid="ignore"):
+ converted_parameters["phase"] = np.mod(
+ converted_parameters["delta_phase"]
+ - np.sign(np.cos(converted_parameters["theta_jn"]))
+ * converted_parameters["psi"],
+ 2 * np.pi
+ )
added_keys = [key for key in converted_parameters.keys()
if key not in original_keys]
@@ -512,7 +514,8 @@ def chirp_mass_and_mass_ratio_to_total_mass(chirp_mass, mass_ratio):
Mass of the lighter object
"""
- return chirp_mass * (1 + mass_ratio) ** 1.2 / mass_ratio ** 0.6
+ with np.errstate(invalid="ignore"):
+ return chirp_mass * (1 + mass_ratio) ** 1.2 / mass_ratio ** 0.6
def component_masses_to_chirp_mass(mass_1, mass_2):
--
GitLab