Xi, Tan and Liu proposed two simple algorithms to generate gamma random numbers based on
the ratio-of-uniforms method and logarithmic transformations of gamma random variable.
This implementation is Algorithm 3.1 in the reference. It does not restrict the range of the
It, however, assumes the scale parameter is always 1.
Note that for when the shape parameter < 0.01, it returns numbers in logarithmic scale.
Take exp to convert.
"B. Xi, K. M. Tan, and C. Liu, "Two New Ratio-of-Uniforms Gamma Random Number Generators,"
Technical Report #10-03, Purdue University."