def gauss(x, sigma):
    """ Compute 1-D value of gaussian at position x relative to center."""
    return (np.exp(-np.power(x, 2) / (2 * np.power(sigma, 2))) /
            (sigma * np.sqrt(2 * np.pi)))