def lnprob(param, x, y, dy):
    # mean, amplitude, stddev = param
    if param[-1] < 0:
        return -np.inf
    yp = gaussian(x, *param)
    diff = (yp - y)
    if dy is not None:
        diff /= dy
    return -np.dot(diff, diff)
dataset/ETHPy150Open glue-viz/glue/doc/gui_guide/emcee_plugin.py/lnprob/swapped operands of `-`