auto rv = gammaVar;
static assert(isRandomVariable!(typeof(rv)));
import mir.random.engine;
auto x = rv(rne);

import mir.random.engine;
Random* gen = threadLocalPtr!Random;
auto rv = GammaVariable!double(1, 1);
auto x = rv(gen);