YNOGK: Calculating null geodesics in the Kerr spacetime

Discussion topics for individual codes
Post Reply
User avatar
owlice
Guardian of the Codes
Posts: 1135
Joined: Wed Aug 04, 2004 4:18 pm
Location: Washington, DC

YNOGK: Calculating null geodesics in the Kerr spacetime

Post by owlice » Mon May 13, 2013 4:20 am

YNOGK: Calculating null geodesics in the Kerr spacetime

Abstract: YNOGK, written in Fortran, calculates the null geodesics in the Kerr spacetime. It uses Weierstrass' and Jacobi's elliptic functions to express all coordinates and affine parameters as analytical and numerical functions of a parameter $p$, which is an integral value along the geodesic. The information about the turning points do not need to be specified in advance by the user, allowing applications such as imaging, the calculation of line profiles or the observer-emitter problem to become root finding problems. Elliptic integrations are computed by Carlson's elliptic integral method, which allows fast computation.

Credit: Yang, Xiaolin; Wang, Jiancheng

Site: http://www1.ynao.ac.cn/~yangxl/yxl.html
https://ui.adsabs.harvard.edu/abs/2013ApJS..207....6Y

Bibcode: 2013ascl.soft05008Y

ID: ascl:1305.008
Last edited by Ada Coda on Thu Aug 20, 2020 10:41 pm, edited 1 time in total.
Reason: Updated code entry.
Code-seeking owl at your service

Post Reply