Speaker
Description
PolyChord
was originally advertised encouraging users to experiment with their own clustering algorithms. Identifying clusters of nested sampling live points is critical for PolyChord
to perform nested sampling correctly. We have updated the Python
interface of PolyChordLite
to allow straightforward substitution of different clustering methods.
Recent reconstructions of the primordial matter power spectrum $\mathcal P_\mathcal R(k)$ with a flex-knot revealed that the K-Nearest-Neighbours algorithm used by PolyChord
cannot reliably detect the two posterior modes caused by cosmic variance and detector resolution. After exploring a number of different algorithms, we have found the X-means algorithm to be a reliable substitute for the power spectrum reconstruction.
This work prompted the development of additions fo the post-processing tool anesthetic
, allowing posterior modes corresponding to different live pint clusters to be analysed and plotted independently.