def logdet(self):
Sigma = np.diag(1/self.J) + np.dot(self.U.T, self.U/self.N[:,None])
cf = sl.cho_factor(Sigma)
ld = np.sum(np.log(self.N)) + np.sum(np.log(self.J))
ld += np.sum(2*np.log(np.diag(cf[0])))
return ld
评论列表
文章目录