18 #ifndef BZ_RANDOM_GAMMA
24 template<typename T =
double, typename IRNG =
defaultIRNG,
30 F(T numeratorDF, T denominatorDF)
32 setDF(numeratorDF, denominatorDF);
33 mindenom = 0.085 * tiny(T());
36 void setDF(T _dfn, T _dfd)
38 BZPRECONDITION(_dfn > 0.0);
39 BZPRECONDITION(_dfd > 0.0);
43 ngamma.setMean(dfn/2.0);
44 dgamma.setMean(dfd/2.0);
49 T xnum = 2.0 * ngamma.random() / dfn;
50 T xden = 2.0 * ngamma.random() / dfd;