Random variate generation for the generalized inverse Gaussian distribution

We provide a uniformly efficient and simple random variate generator for the entire parameter range of the generalized inverse Gaussian distribution. A general algorithm is provided as well that works for all densities that are proportional to a log-concave function φ, even if the normalization constant is not known. It requires only black box access to φ and its derivative.