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.