def log_prob(self, x): from numpy.linalg import det mu = self.mu S = self.sigma D = len(mu) q = self.__q(x) return -0.5 * (D * log(2 * pi) + log(abs(det(S)))) - 0.5 * q ** 2