===================
- Add ``bayestar-mcmc`` tool for pure Markov Chain Monte Carlo parameter
estimation, without sky map postprocessing but with options for holding
parameters at fixed values.
- Fix a corner case in the initialization of the ``distance.marginal_ppf``
solver that could cause NaN return values.
- Silence ``numpy.genfromtxt`` Unicode deprecation warning in
``ligo-skymap-plot-stats`` and update the minimum version of Numpy to 1.14.
See the related `Numpy changelog entry
<https://docs.scipy.org/doc/numpy/release.html#encoding-argument-for-text-io-functions>`_.
- Silence deprecation warning in ``ligo-skymap-plot-stats`` due to Matplotlib
renaming the ``hist`` method's keyword argument from ``normed`` to
``density``.
- The ``bayestar-realize-coincs`` tool now copies over spins from the input
``sim_inspiral`` table to the output ``sngl_inspiral`` table.
- Switch the FFT implementation from LAL (which calls `FFTW
<http://www.fftw.org>`_) to `scipy.fftpack
<https://docs.scipy.org/doc/scipy/reference/tutorial/fftpack.html>`_, which
is faster for small transform sizes (e.g. <= 1024).
- Add ``--seed`` option to ``bayestar-localize-coincs``,
``bayestar-localize-lvalert``, ``bayestar-mcmc``, and
``bayestar-realize-coincs``.
- Some reasonable sub-sample trigger interpolation schemes can return peak
times that are almost a full sample away from the maximum sample if the SNR
time series has a pronounced skew in one direction in the vicinity of the
maximum. Such an example occurs for the ``catmull-rom`` interpolation method
for the new unit tests in ``ligo.skymap.bayestar.tests.test_interpolation``.
Because of this, relax the tolerance of BAYESTAR's sanity check on
single-detector trigger times and SNR series timestamps to a full sample.
- Rewrite ``ligo-skymap-plot-stats`` to reduce code duplication.
- Add ``--measurement-error gaussian-noise`` option to
``bayestar-realize-coincs`` to simulate a matched filter in Gaussian noise.
- Remove deprecated module ``ligo.skymap.postprocess.detector_frame``.