This function samples
sim_gig(num_sim, lambda, psi, chi)
Number to generate
Index of modified Bessel function of third kind.
Second parameter of GIG. Should be positive.
Third parameter of GIG. Should be positive.
The density of
Hörmann, W., Leydold, J. Generating generalized inverse Gaussian random variates. Stat Comput 24, 547-557 (2014).
Leydold, J, Hörmann, W.. GIGrvg: Random Variate Generator for the GIG Distribution. R package version 0.8 (2023).