I developed the batman package for fast calculation of transit light curves in Python. The package supports calculation of light curves for any radially symmetric stellar limb darkening law, using a new integration algorithm for models that cannot be quickly calculated analytically.

In typical use, batman can calculate a million model light curves in well under 10 minutes for any limb darkening profile.

For more info, check out the docs.