def __call__(self, x): from scipy.special import gammaln return self.a * x - \ self.n * gammaln(numpy.clip(x, 1e-308, 1e308)) + \ self.b * log(x), \ self.a - self.n * psi(x) + self.b / x