def compute_phi_cavity(self):
phi = 0
for d in range(self.Dout):
muhatd_val = self.muhat[d].get_value()
Suhatd_val = self.Suhat[d].get_value()
# s, a = npalg.slogdet(Sud_val)
a = np.log(npalg.det(Suhatd_val))
phi += 0.5 * np.dot(muhatd_val.T, np.dot(npalg.inv(Suhatd_val), muhatd_val))[0, 0]
phi += 0.5 * a
return phi
评论列表
文章目录