def _logp_tau(self, value):
chol_tau = self.chol_tau
k = chol_tau.shape[0]
mu = self.mu
delta = value.reshape((-1, k)) - mu
delta_trans = tt.dot(chol_tau.T, delta.T)
After Change
diag = tt.nlinalg.diag(chol_tau)
delta_trans = tt.dot(chol_tau.T, delta.T)
result = n * k * tt.log(2 * np.pi)
result -= 2.0 * n * tt.sum(tt.log(diag))
result += (delta_trans ** 2).sum()
result = -0.5 * result
return bound(result, ok)
class MvStudentT(Continuous):
R