++++++++++++++++++
- Second release of the code for resubmission to AJ.
- Major code re-write. Redesigned user interface, easier to use,
faster, and more flexible.
- Maps are now instantiated via a single `Map` object that supports
both spherical harmonic coefficients and limb darkening coefficients,
as well as arbitrary combinations of both.
- The `Map` class now accepts a `nwav` keyword specifying the number of
wavelength bins (default 1). This allows users to easily and efficiently
compute wavelength-dependent light curves from spectral surface maps.
- The `Primary` and `Secondary` classes replaced the old `Star` and
`Planet` classes for increased generality. These are now actual subclasses
of `Map`, making them even easier to use.
- Several modifications to the computation of the light curves were made to
increase speed. The largest speed improvement is in the computation of
gradients. The code is no longer limited by the `STARRY_NGRAD` parameter,
and arbitrary number of gradients can now be taken efficiently.
- Derivatives with respect to the map coefficients are no longer labeled
as `Y_{0,0}`, `Y_{1,-1}`, ... Instead, they are labeled by a single key
`y`, whose value is a vector of derivatives corresponding to each element
in the spherical harmonic vector `y`. The same applies to the limb
darkening coefficients, whose derivatives are stored in they key `u`.
- Minor changes to the names of some keywords and class properties, and
to some call sequences.
- Minor bug fixes here and there.