function randomGamma(k, theta) {
if (k === 1) {
return d3r.randomExponential(1 / theta);
}
var U = d3r.randomUniform(),
N = d3r.randomNormal(),
d = (k < 1 ? k + 1 : k) - 1 / 3,
c = 1 / (3 * Math.sqrt(d)),
multiplier = k < 1 ? () => Math.pow(U(), 1 / k) : () => 1;
return function() {
do {
do {
var x = N(),
v = 1 + c * x;
} while (v <= 0);
v *= v * v;
var u = 1 - U();
} while (u >= 1 - 0.0331 * x * x * x * x && Math.log(u) >= 0.5 * x * x + d * (1 - v + Math.log(v)));
return d * v * multiplier() * theta;
};
}