def _t(self):
a = self._alpha_s
ad = self._alphad_s
zd = self._zd
pg = sp.polygamma(0, zd)
return (a * self.Hsd / self.betad) * self._gammazd * (self._yd ** (ad + 1) * self._gammazd *
sp.hyperReg_2F2(zd, self._xd) +
pg - self.betad * np.log(self._yd))
def _lnk_a_a(self):
return sp.polygamma(1, self._zk)/self.beta**2