def METHOD_NAME(x, atol=1e-2, rtol=0.1):
    # This noise affects just a single component
    sqrt_n = np.sqrt(float(np.prod(x.shape)))
    tol = 2 * (np.linalg.norm(x) * rtol + atol * sqrt_n)
    n = np.random.randint(np.prod(x.shape))
    noise = np.zeros_like(x)
    noise.reshape(-1)[n] = tol
    return x + noise