This adds a mixture-model power law distribution, useful for hyperparameter problems.
>>> a = bilby.core.prior.DoublePowerLaw(alpha0=-1, alpha1=1, xi=0.3, minimum=10, maximum=100) >>> plt.hist(a.sample(10000), bins=30) >>> plt.show()