Efficient generation of microlensing magnification maps with GPUs

ABSTRACT The impending discovery and monitoring of hundreds of new gravitationally lensed quasars and supernovae from upcoming ground and space-based large-area surveys such as LSST, Euclid, and Roman necessitates the development of improved numerical methods for studying gravitational microlensing. We present in this work the fastest microlensing map generation code currently publicly available. We utilize graphics processing units to take advantage of the inherent parallelizable nature of creating magnification maps, in addition to using (1) the fast multipole method to reduce the runtime dependence on the number of microlenses and (2) inverse polygon mapping to reduce the number of rays required. The code is available at https://github.com/weisluke/microlensing/.