We constrain extensions of general relativity that can be described by a reflection coefficient as described in https://arxiv.org/abs/2001.00041.
Python 3.7, SciPy 1.2
python/main.py includes all the relevant routines. These are:
Create interpolated digitized EHT image. The csv files characterizing the digitized EHT image are stored in digitize_eht/.
Determine properties of image such as ring diameter, width, deviation from circularity, fractional central brightness, etc.
Get intensities of reflected light rays. Includes kerr_geodesic which calculates the Kerr geodesic of a light ray using odeint.
Add reflected image to the original EHT image.
Check if a modified gravity image is consistent with the EHT observtions.
Creates exclusion plots.