implement missing sampling functions for non-Gaussian distributions

left over from recent refactoring