'Compute the differential entropy of the dirichlet distribution.
Parameters
%(_dirichlet_doc_default_callparams)s
Returns
h : scalar
Entropy of the Dirichlet distribution'

def entropy(self, alpha):

    alpha = _dirichlet_check_parameters(alpha)
    alpha0 = np.sum(alpha)
    lnB = _lnB(alpha)
    K = alpha.shape[0]
    out = ((lnB + ((alpha0 - K) * scipy.special.psi(alpha0))) - np.sum(((alpha - 1) * scipy.special.psi(alpha))))
    return _squeeze_output(out)
