Comparison of null geodesics computed by the semi-analytical integrator geokerr to those computed by our fully numerical code. In each case, the geodesic represented by the black dots was computed using geokerr while the geodesic represented by the white dotted line is the output of RAPTOR (dots were chosen for the former geodesic because its variable step-size can otherwise create interpolation issues). The gray spheroid represents the black hole’s outer event horizon. MKS coordinates were used with ϵ = 0.001. The parameters for these geodesics are listed in Table 1.