diff options
-rw-r--r-- | null.py | 4 |
1 files changed, 3 insertions, 1 deletions
@@ -89,8 +89,10 @@ def _null_geodesic_kernel(times, t0, events, pu_t_0, gtX, gtX_dX, gtX_dt): gtt_0 = gtt(t0, 0.0).flatten()[0] gXX_0 = gXX(t0, 0.0).flatten()[0] + gtX_0 = gtX(t0, 0.0).flatten()[0] - pu_X_0 = np.sqrt(-gtt_0 / gXX_0) * pu_t_0 + term_t = gtX_0 * pu_t_0 + pu_X_0 = (-term_t + np.sqrt(term_t * term_t - gtt_0 * gXX_0 * pu_t_0 * pu_t_0)) / gXX_0 args = { # workaround for older versions of scipy not supporting args for fun |