Unified MCMC interface
MCMC implementation of toa_phoa_snr
No need to set GSL error handler here

These functions do not raise GSL errors.
bayestar_realize_coincs requires PSD input

and has more sensible default arguments
bayestar_realize_coincs gets a progress bar