===================
- The ``bayestar-inject`` script now assumes that the source distribution is
specified per unit comoving volume per unit proper time, rather than per unit
comoving volume per unit observer time. This is in agreement with the
conventional definition for LIGO/Virgo astrophysical rates.
- The ``bayestar-inject`` and ``ligo-skymap-from-samples`` scripts now accept
an optional integer value for the ``-j`` flag to set the number of
subprocesses.
- ``ligo-skymap-from-samples`` will use all posterior samples if the value of
the ``--maxpts`` argument is greater than or equal to the number of posterior
samples.
- If the ``billiard`` package is present, then use it instead of the
``multiprocessing`` standard library module to parallelize
``ligo-skymap-from-samples`` so that the script's Python entry point can
be called from daemon processes (for example, inside Celery tasks).
- Switch from WMAP9 to Planck15 cosmological parameters.
- ``ligo.skymap.kde.Clustered2DSkyKDE.as_healpix()`` has an optional
``top_nside`` to allow for better initial grid, before refinement.
``ligo-skymap-from-samples`` has an additional ``--top-nside`` argument,
accordingly.