We discuss methods for sampling Generalized Inverse Gaussian (GIG) distributed random variates both in the fixed and in the varying parameter case.
In particular we present a new algorithm for the varying parameter case. It is based on the acceptance-rejection method and has, different to all algorithms presented in the literature, a uniformly bounded rejection constant.