performance improvements for generating pareto/NBD draws
added implementation of Gamma/Gompertz/NBD model (by Bemmaor/Glady)
implemented custom slice sampler (to improve performance)