A new approach and code for spinning black holes in modified gravity

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.