def log_prob(self, x): a, b = self['alpha'], self['beta'] return a * log(b) - gammaln(clip(a, 1e-308, 1e308)) + \ (a - 1) * log(clip(x, 1e-308, 1e308)) - b * x