def _h_ars(x,params):
    """ln p(Omega) for ARS"""
    mO, sO2= params
    return -0.5*(x-mO)**2./sO2+numpy.log(x)
def _h_ars ( x , params ) : mO , sO2 = params return - 0.5 * ( x - mO ) ** 2. / sO2 + numpy . log ( x )
ln p(Omega) for ARS