distributions.py 文件源码

def gaussian_logpdf(mu,sigma,x):
    '''
    Non-positive standar deviations will be clipped

    Parameters
    ----------

    Returns
    -------
    '''
    mu,sigma,x = map(np.float128,(mu,sigma,x))
    x = (x-mu)/sigma
    return -0.5*x*x - slog(sigma) - logsqrt2pi