GaussObsModel.py 文件源码

def c_Func(nu, logdetB, m, kappa):
    ''' Evaluate cumulant function at given params.

    Returns
    --------
    c : scalar real value of cumulant function at provided args
    '''
    if logdetB.ndim >= 2:
        logdetB = np.log(np.linalg.det(logdetB))
    D = m.size
    dvec = np.arange(1, D + 1, dtype=np.float)
    return - 0.5 * D * LOGTWOPI \
           - 0.25 * D * (D - 1) * LOGPI \
           - 0.5 * D * LOGTWO * nu \
           - np.sum(gammaln(0.5 * (nu + 1 - dvec))) \
        + 0.5 * D * np.log(kappa) \
        + 0.5 * nu * logdetB