An improved version of this code base can be found in gh-levy-simulation.
The implementation of the methodology introduced in Point process simulation of Generalised inverse Gaussian processes and estimation of the Jaeger Integral by Simon Godsill and Yaman Kındap is provided.