Abstract
We discuss and implement a spectral method approach to computing stationary and axisymmetric black hole solutions and their properties in modified theories of gravity. The resulting code is written in the Julia language and is transparent and easily adapted to new settings. We test the code on both general relativity and on Einstein-scalar-Gauss–Bonnet gravity. It is accurate and fast, converging on a spinning solution in these theories with tiny errors (${\sim} \mathcal{O}\left(10^{-13}\right)$ in most cases) in a matter of seconds.
Licence information
This item is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Attribution 3.0 United States
Copyright statements
© 2023 The Author(s). Published by IOP Publishing Ltd