def METHOD_NAME(x, atol=1e-2, rtol=0.1):
    # We go in random direction using twice the original tolerance to be sure this
    # results in an error
    sqrt_n = np.sqrt(float(np.prod(x.shape)))
    tol = 2 * (np.linalg.norm(x) * rtol + atol * sqrt_n)
    noise = np.random.normal(size=x.shape)
    noise = tol * noise / np.linalg.norm(noise)
    return x + noise