diff options
author | Anton Khirnov <anton@khirnov.net> | 2023-03-16 14:00:20 +0100 |
---|---|---|
committer | Anton Khirnov <anton@khirnov.net> | 2023-03-16 14:00:20 +0100 |
commit | add65ec667d42fd7a63ac64d328c3e6efe78b618 (patch) | |
tree | 20a7047937ab6c38a398ade9403ee7afbaba581f | |
parent | cfe507106e8c7f182c48690acd2b286c303564b0 (diff) |
null: fix the null-vector constraint with non-zero shift
-rw-r--r-- | null.py | 3 |
1 files changed, 2 insertions, 1 deletions
@@ -62,7 +62,8 @@ def _null_geodesic_axis_rhs(Lambda, f, gu = np.linalg.inv(g) # enforce the momentum to be a null vector - mom_t = np.sqrt(-gXX / gtt) * mom_X + termX = gtX * mom_X + mom_t = (-termX + np.sqrt(termX * termX - gtt * gXX * mom_X * mom_X)) / gtt G = np.zeros((4, 4, 4)) G[0, 0, 0] = 0.5 * gtt_dt |