From 42caa1c66beed263c9d512a9ee3a6d8277b6b2f6 Mon Sep 17 00:00:00 2001
From: Colm Talbot <colm.talbot@ligo.org>
Date: Thu, 28 Feb 2019 07:36:16 +1100
Subject: [PATCH] use more time samples in interpolation
---
bilby/gw/conversion.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/bilby/gw/conversion.py b/bilby/gw/conversion.py
index dd63232fd..e4ffcff23 100644
--- a/bilby/gw/conversion.py
+++ b/bilby/gw/conversion.py
@@ -1107,7 +1107,7 @@ def generate_time_sample_from_marginalized_likelihood(
time_prior_array = likelihood.priors['geocent_time'].prob(times)
time_post = (np.exp(time_log_like - max(time_log_like)) * time_prior_array)
- keep = (time_post > max(time_post) / 10)
+ keep = (time_post > max(time_post) / 1000)
time_post = time_post[keep]
times = times[keep]
sample['geocent_time'] = Interped(times, time_post).sample()
--
GitLab