Basic BS working with a small glitch: the pdf functions must return 0-dimensional arrays.
BS works properly now, estimates log Z from log posterior values instead of Z from posterior values.
Added a very simple example notebook.