MulensModel calculates light curves of microlensing events. Both single and binary lens events are modeled and various higher-order effects can be included: extended source (with limb-darkening), annual microlensing parallax, and satellite microlensing parallax. The code is object-oriented and written in Python3, and requires AstroPy (ascl:1304.002).
https://ui.adsabs.harvard.edu/abs/2018ascl.soft03006P and https://ui.adsabs.harvard.edu/abs/2019A%26C....26...35P ; please see additional citation information here: https://github.com/rpoleski/MulensModel#acknowledgements