docs: def get_entropy(self):
          """Computes (analyticaly) the entropy of the Gaussian distribution.
          """
          dim = self.mu.shape[0]
          entropy = 0.5 * (dim * (numpy.log(2. * numpy.pi) + 1.) + numpy.log(self.det))
          return entropy