Skip to content

Commit 081496e

Browse files
minor in coordchange (#56)
* minor in coordchange * [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci * added also pe in output * [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci --------- Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>
1 parent 8bfeef7 commit 081496e

1 file changed

Lines changed: 10 additions & 3 deletions

File tree

PyART/analytic/coordschange.py

Lines changed: 10 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -271,6 +271,10 @@ def eob_ID_to_ADM(eob_Wave, verbose=False, PN_order=2, rotate_on_x_axis=True):
271271
'x2' : x coordinate of puncture 2
272272
'D' : coordinate separation between punctures
273273
'x_offset' : offset to be added to x coordinates
274+
'qe' : EOB coordinates (x,y)
275+
'pe' : EOB momenta (px, py), mu-normalized
276+
'qe_chk' : EOB coordinates (x,y), EOB->ADM->EOB
277+
'pe_chk' : EOB momenta (px, py), mu-normalized, EOB->ADM->EOB
274278
"""
275279
# Get info from EOB dynamics
276280
q = eob_Wave.pars["q"]
@@ -308,6 +312,8 @@ def eob_ID_to_ADM(eob_Wave, verbose=False, PN_order=2, rotate_on_x_axis=True):
308312
x2 = -d_ADM * q / (q + 1)
309313
x_offset = -b_par + d_ADM / (q + 1)
310314

315+
qe_check, pe_check = Adm2Eob(qa, pa, nu, PN_order=PN_order)
316+
311317
# wrap output
312318
out = {
313319
"q_cart": qa,
@@ -318,12 +324,13 @@ def eob_ID_to_ADM(eob_Wave, verbose=False, PN_order=2, rotate_on_x_axis=True):
318324
"x2": x2,
319325
"D": d_ADM,
320326
"x_offset": x_offset,
327+
"qe": qe,
328+
"pe": pe,
329+
"qe_chk": qe_check, # EOB->ADM->EOB
330+
"pe_chk": pe_check, # EOB->ADM->EOB
321331
}
322332

323333
if verbose:
324-
# for testing
325-
qe_check, pe_check = Adm2Eob(qa, pa, nu, PN_order=PN_order)
326-
327334
dashes = "-" * 50
328335
print("{}\nPunctures\n{}".format(dashes, dashes))
329336
print("b_par : {:.15f}".format(b_par))

0 commit comments

Comments
 (0)